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The fundamental solution for small steady three- 
dimensional disturbances to a two-dimensional parallel 
shear flow 
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1. INTRODUCTION 
In addition to the enormous aerodynamical literature concerned with 
the disturbances to a uniform stream produced by the presence of obstacles, 
there is a very extensive literature on disturbances to irrotational non- 
uniform streams. ‘The smaller body of work on disturbances to rotational 


streams is concerned principally with two-dimensional disturbances to a 
parallel shear flow; these are fairly easy to treat because the vortex lines 
are straight; they remain parallel and do not get stretched at all during the 


motion. 

3y contrast, in flows where the disturbances are three-dimensional in 
character (even though the undisturbed stream may still be two-dimensional), 
stretching and rotation of the vortex lines must play an important part. 
In spite of the difficulties which this introduces, certain outstanding papers 
on the theory of these flows have appeared. 

Squire (1933) found the equation governing small three-dimensional 
disturbances to a two-dimensional parallel flow. If the velocity field is 


(V(y) +4, v, w), (1) 


H 





114 M. Ff. Lighthill 


then by neglecting squares of the disturbance velocities u, v, w and analysing 
them into Fourier components, he showed that disturbances varying 
with x and 2 like e"** must vary with time like e~®, where c is such that 
the equation 
(V —c){v" — (a? + B*)v} — V"v — (v/ia){v'* — 2(a? + B?)o” + (a? + B?)?v} =0 (2) 
(where v is kinematic viscosity and where primes signify differentiation 
with respect to y) is satisfied (together with boundary conditions) by a 
non-vanishing value of the y-velocity v. Values of c with negative imaginary 
part correspond to stable disturbances, values with positive imaginary part to 
unstable disturbances, and real values to neutral disturbances. Equation (2) 
shows that neutral disturbances with 8 = By ~ 0, « = %, v = vp can exist, 
if and only if neutral disturbances with 

B = 0, a = 4/(a)+ Bi), v = vo/(1 + Bo/%) 
exist; the latter are two-dimensional disturbances at a lower Reynolds 
number*., Squire deduced the classical theorem that two-dimensional 
disturbances become unstable at a lower Reynolds number than three- 
dimensional disturbances. 

Karman & ‘Tsien (1945) investigated small disturbances to a three- 
dimensional parallel stream (so that V in (1) is a function of 2 as well as 
of y) but limited themselves (as we shall do in this paper) to steady 
disturbances with viscosity neglected, evolving the theory of the ‘lifting 
line’ in such a non-uniform stream. ‘They found it convenient to work 
from the equation for the pressure p, namely 

divs") =f, (3) 

% 
both because it in this case is simpler than that for any of the velocity 
components, and because the boundary condition in the lifting-line problem 
is expressed most conveniently in terms of p. The results of the lifting-line 
theory are complicated, however, and only in the case when V is a function 
of y alone (the spanwise coordinate) do they become reasonably tractable. 

Another approach to the study of steady three-dimensional disturbances 
to a parallel shear flow has been to regard them as a ‘secondary flow’, to be 
analysed by the method initiated by Squire & Winter (1951). In this 
method there is no assumption that the disturbances are small; for example, 
the disturbance due to a sphere can be treated (Hawthorne & Martin 1955; 
Lighthill 1956). ‘There is, however, an assumption that the undisturbed 
stream is weakly sheared. ‘The ‘primary flow, (or first approximation to 
the real flow) is taken to coincide (either wholly, or at least as regards its 
streamline pattern) with that in which the undisturbed stream is uniform. 
The secondary flow is a perturbation on this, due to allowing a small shear 
in the undisturbed stream. ‘This shear is usually taken to be uniform, 
although such an assumption is not essential to the method; we shall call 
secondary flows calculated on this assumption ‘simple-shear secondary 
flows’. The secondary vorticity can be calculated either direct from 


* Equation (2) with 8 = 0 is the Orr-Sommerfeld equation. 





ee _cl eeenm 





eS 





1 


\ 
? 





EO 





SO 


Steady three-dimensional disturbances to a parallel shear flow 115 


Helmholtz’s equations (Hawthorne 1951, 1954), or by geometrical con- 
siderations (Preston 1954; Lighthill 1956) from the effect of the primary 
flow in stretching and rotating the vortex lines of the undisturbed shear flow. 
(From this point of view, the approximation adopted consists in neglecting 
the fact that it is really the exact flow, not the approximate ‘ primary flow’, 
which carries the vortex lines alongwith it.) ‘The calculation of the secondary 
velocity field from the associated vorticity field is more arduous, but can be 
achieved in special cases (Lighthill 1956, 1957 b). 

The work of this paper originated from a desire to clear up a difficulty 
arising in the application of the secondary-flow method to cases where the 
undisturbed stream is unbounded. The difficulty is that the secondary-flow 
disturbance due to the presence of an obstacle falls off more slowly with 
distance from it than does the primary-flow disturbance. For example, 
a ‘half-body’ extending from the origin to infinity downstream has a 
primary-flow disturbance of ‘source’ type, falling off like r~-* (where r is 
distance from the origin), but the corresponding secondary-flow disturbance 
falls off like rt. If the procedure of successive approximation were 
continued, and a tertiary-flow field calculated from the configuration of 
vortex lines as distorted by the sum of the primary and secondary flows, 
this would not even tend to zero as r > ©. 

As another example, a finite body makes a primary-flow disturbance of 
‘doublet’ type, falling off like r-*, but the corresponding secondary-flow 
disturbance falls off like 7~* (see Lighthill (1956), where also a not quite 
accurate attempt was made to calculate the tertiary flow, but the conclusion 
that it falls off like r-! was right). 

Now, there is no experimental evidence (e.g. from the flow visualizations 
reported by Livesey (1956)) that secondary flows are large far from the 
obstacle, and it seems clear that, as suggested in Lighthill (1956), the 
approximation sequence is not uniformly valid in this region, where 
(therefore) the first term gives no indication of the real behaviour of the 
disturbance. 

There is some analogy in this to the difficulties concerning flow about an 
obstacle at very low Reynolds numbers, often described as ‘ Stokes’s and 
Whitehead’s paradoxes’. In those problems a direct procedure of expansion 
in powers of Reynolds number gives nonsensical results, at any rate far from 
the obstacle (although for three-dimensional obstacles the first term in the 
expansion has value near the obstacle). The difficulty was resolved by 
Oseen, who showed that the correct form of the equations of motion at large 
distances is that obtained by neglecting squares of disturbances to the 
uniform stream but not discarding terms of lower order in the Reynolds 
number. 

Similarly in the present problem the true behaviour of the flow far 
from the body is found from the full equations with the squares of the 
disturbances neglected. Such an approach is obviously valid sufficiently 
far from the obstacle, and the only problem is how the results obtained 
should be related to secondary-flow and other results obtained by the 
methods described above. 

H2 
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The answer found below is that the small-disturbance solution is valid 
in a region, far from the obstacle, which overlaps the region, near the 
obstacle, where the small-shear solution (more accurately, the ‘ simple-shear 
secondary flow’) is valid. ‘This overlapping holds in the sense that the 
asymptotic behaviour of the small-disturbance solution as one goes near 
the obstacle is identical with the asymptotic behaviour of the small-shear 
solution as one goes away from the obstacle. ‘The overlap makes it possible 
to say how the flow behaves everywhere. 

These are good reasons for studying small steady three-dimensional 
disturbances to a parallel shear flow. Obviously, however, there are many 
others: to investigate, for example, the flow of sheared winds over hills. 
Again, one could follow up the work of Karman & ‘Tsien (1945) by studying 
further kinds of perturbation of parallel shear flows by thin shapes of 
aerodynamic interest. 

As a basic tool for any such investigation, we seek in this paper the 
‘fundamental solution’ of the equation. ‘This represents the flow due to 
a source (but a weak source, since our equation treats only small disturbances) 
in a parallel shear flow. ‘This solution will give as it stands the asymptotic 
behaviour of the disturbances due to a ‘ half-body ’, since these are equivalent 
far from the body to those produced byasource. Also, a simple diff-rentiation 
of the solution will give the disturbance due to a doublet, which wit represent 
the asymptotic behaviour of the disturbances due to a finite body. Finally, 
more complicated solutions of the equations can be built up (on familiar 
lines) once the fundamental solutions for a source and a doublet are known; 
work on these lines, however, is excluded from the present paper. 

We limit ourselves (like Squire 1933) to the case when the undisturbed 
flow is two-dimensional, its velocity field being (V(y),0,0). ‘The formal 
theory for the undisturbed velocity field (V(y, z),0,0) (used by Karman 
& Tsien 1945) can be derived, but so much extra effort is then needed to 
deduce any concrete results (as indeed they also found) that it is better to 
try to keep the theory intelligible by confining discussion to the simpler 
case. In this case the equation, after Fourier transforms with respect to x 
and z have been taken, becomes Squire’s equation (2) with c = 0 (because 
the disturbances are steady). ‘There are, of course, no continuous solutions 
of this equation which satisfy the appropriate boundary conditions (such 
solutions exist only for one value of c, which in practice is non-zero); on 
the other hand, to represent a source (say, at the origin), one needs a solution 
with a discontinuity of a certain specified kind at y = 0. ‘This solution is 
the Fourier transform of the required disturbance field, which must be 
deduced from it by Fourier’s inversion theorem. 

The Orr-Sommerfeld equation has in the past been used almost 
exclusively for studying harmonic disturbances. ‘The building up of 
anharmonic disturbance fields by Fourier synthesis from solutions of that 
equation, as in this paper and the author’s previous paper on boundary 
layers and upstream influence in supersonic flow (Lighthill 1953), is an 
approach that seems fruitful and may have a number of applications. 
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In the present paper the viscous terms in the equatior. are neglected. 
This is known to be permissible except where V(y) is near the ‘critical’ 
values 0 and c, which coincide in the case here studied. ‘Thus, shear layers 
in which the undisturbed velocity falls to zero anywhere are excluded in 
this paper as a result of the neglect of the viscous terms; such shear layers 
include boundary layers and jets. Boundary layers could be treated, 
provided that the viscous terms in the equation were retained in some 
‘inner viscous layer’ near the wall, as used by Lighthill (1953) to provide 
a corrected form of boundary condition for the non-viscous equation. 
However, in a preliminary survey this refinement was thought to be not 
worth making, particularly since laminar boundary layers when disturbed 
are so prone to large disturbances associated with separation. With jets, 
again, the assumption of parallel or even nearly parallel flow at the edge, 
where air is being entrained, is wrong. Accordingly, the theory as set out 
here is limited to a wake, or a mixing region between parallel streams, or 
any other layer in which V(y) is nowhere zero. 

Since steady inviscid flow is being considered the Karman—T’sien 
equation (3) could have been used. However, equation (2) for v (with 
v=c=(), although no simpler than equation (3) for p, has advantages 
in the present problem, in which v (but not p) turns out to be a function 
of vy and of s = 4/(x2+=37) alone, and in which the velocity components, 
rather than the pressure, are the main things one wants to know. 

Of the results of the theory, four main classes may be mentioned in 
advance. First, if the solution is expanded in ascending powers of r (distance 
from the source), the first term (of order r *) is the primary flow (the source 
flow itself) and the second (of order r~') is the small-disturbance approxi- 
mation to the simple-shear secondary flow (as already mentioned) and 
depends only on V(0) and V’(0), the velocity and shear at the origin. 
However, the next term (which might be called the small-disturbance 
form of the ‘tertiary flow’) depends on the velocity distribution throughout 
the shear layer. 

Another result of general interest is concerned with the behaviour of the 
solution for large r: a source in a shear layer produces in any region of 
uniform flow outside the shear layer a disturbance equivalent to a source of 
different strength in a different position. ‘The strength of the equivalent 
source can alternatively be predicted by simpler arguments (due to 
Mr M. B. Glauert), in which the shear iayer is regarded as a discontinuity 
between regions of uniform flow. ‘The displacement in effective position 
is of the order of the layer width, and has not yet been predicted by such 
simpler arguments. 

A third result is that in shear layers with small velocity spread (that is, 
a velocity minimum only a few per cent less than the velocity maximum) 
a specially simple form of solution is available; this can be given an elegant 
interpretation by means of images, again due to Mr M. B. Glauert. It can 
also be shown to be identical with the asymptotic form of the ‘ exact-profile 
secondary flow’, that is, of the secondary flow calculated by allowing the 
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vortex lines of the exact undisturbed-flow profile V() (not the corresponding 
simple-shear profile (0)+V'(O)v) to be stretched and rotated by the 
primary flow. 

Lastly, the theory is applied to the problem of the displacement of the 
dividing streamline in shear flows about axisymmetrical obstacles—a 
problem which the author has used in a number of recent papers, both as 
a convenient excuse for developing techniques for calculating shear flows, 
and as a point where theory could easily be referred to the results of 
experiments. In the last section of this paper it is shown what corrections 
need to be made, to values of this displacement as calculated from simple- 
shear secondary-flow theory, to take into account the improved picture of 
what the disturbances are like far from the obstacle. 


2. EQUATIONS OF MOTION 


The momentum equation for steady incompressible inviscid flow is 
1 

v.Vv+-—Vp =), (4) 
p 


where v is the velocity, p the density and p the pressure. If the departure 
from two-dimensional parallel flow is small, so that 

v = (V(y)+u, v, w), (5) 
where squares and products of u,v, w and their derivatives are to be 
neglected, equation (4) becomes 





| Ou ' 1 op 
if Vy — +e)’ "4 +-—= 0, 6 
(5, V0) + 5 5 (6) 
400 12 = 
Vy) + - + = 0, (7) 
3x p ay 
~~ oe ie 
Mik eae (8) 
OO p os 
We shall assume that all disturbances vanish far upstream (at x = — %), 
so that (7) and (8), with the pressure eliminated, give 
9 1p? a Se 
= {V(y)w} = — - —— dx = V(y)—. 9 
oy U (3 Ju p | : aya “ (¥) a ( ) 
Elimination of p from (6) and (7) gives, similarly, 
e ou | 072 el 
od Blaha be Behe. = oe Sel. 10 
oy I (3 Ox J I (Yas oy UV (y)z J ( ) 


These equations are to be satisfied together with the equation of 
continuity. Now, to obtain the ‘fundamental solution’, due to a source 
of strength m at the origin, from which other solutions can be built up, we 
write that equation as 


+ “ + = m6(x)5(y)d(z), (11) 
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where 6(x) is the delta function of Dirac, or ‘impulse function’; thus, 
the right-hand side of (11) vanishes except at the origin, but its integral 
over any volume including the origin is m. No similar term appears on the 
right-hand side of (6), (7) or (8) because the source is taken to be a source 
of mass but not of momentum. 

To obtain an equation for the single variable v, we first multiply (11) 
by V(y), giving 


Vy) = + Vy) = + V (9) = = mV (0)3(x)8(y)8(2) (12) 
(where the general theorem f(¥v)5(v) = f(0)6(y) has been used), and then 


we differentiate (12) with respect to y, using (9) and (10), to obtain 


A A ~ A 





SV Oe} + 55 VO) E+ Va) sy = mV (0)8(4)5'()8(2) 
. . . (13) 
or V(y)V20 — V'"(y)e = mV(0)8(x)3'(y)8(s). (14) 


Equation (14) for v is the basic equation of motion. After it has been 
solved, w and u may easily be deduced, from (9) in the form 


uv “ 


V(y)w = | V(5 y) “dy, (15) 
Os 

and from the equation of continuity (11). ‘The choice of lower limit in (15) 

arises from the condition that the disturbance vanish at y= + , but 

either limit gives the same answer. ‘This is because 


| Vive dy =0. (16) 


To prove (16) (which is a useful check on the accuracy of solutions, and has 
been applied as such to all those given in this paper), integrate (13) from 
y=-—-%0% to y= +0. This shows that the left-hand side of (16) is a 
harmonic function of the two variables x and z for all values of those 
variables; but it vanishes at infinity, and so must be zero everywhere. 
To obtain u, on the other hand, one must integrate the equation of 
continuity with respect to x from — © (not + 0) to x, since it is possible 
that fluid at x = + 0, whose vorticity may have been permanently re- 
distributed when it passed the source, may have non-zero disturbance 
velocity. 
It is possible that the solution of the basic equation (14), for given V(y), 
could be obtained by some sort of ‘analogue method’. Thus, the equation 
can be written 


div, V*(y)gr rad 7 T% To} = mV(0)d(x)s'(v)8(z), (17) 
which shows that v/I’(y) is proportional to the electrostatic potential of a 
dipole at the origin, with axis the y-axis, in a stratified medium with dielectric 
constant proportional to V*(y). Many similar analogies are, of course, 
possible, and some might, perhaps, be useful in giving intuitive information 
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about the solution, or an experimental procedure for computing it. ‘The 
author has not so far been able to make such use of them, and in the meantime 
offers an analytical method of solution. 

It may also be noted that the theory can be at once extended to probleins 
with density stratification, in which the undisturbed flow has its density p(y) 
as well as its velocity (V(v), 0,0) dependent on y; however, this extension 
is straightforward only if gravity is still negligible—thus, the variable inertia 
of the fluid, but not the variable weight, can easily be taken into account. 
The equation which results in this case, corresponding to (17), is 

div, V*( yy? grad —— $ 
~“ p(Q) V(y) 


so that the expression for v/l(v) may be derived from its expression for 


= mV(0)6(x)d'(¥v)d(z), (18) 


constant p by simply replacing V’(y) in the latter expression by 
Vy) Vv te(’)/p(0)}- 

Unfortunately this result is not of much value because in problems with 
density stratification the variable buoyancy of the fluid is usually important 
(or else the Mach number is not negligible). One may hope that the theory 
can later be extended to take the gravitational term into account, but this 
would require a new investigation. 

To solve equation (14) it is usually convenient to divide it by V(y), 
giving 





oe Oe. ree, os ee ( 
V°7 7 (y) 7 mo(x)o(z)< 6 ( yt V(0) d() es (1 )) 
where the form of the right-hand side follows from the general theorem 
" f(¥)8'(y) = F(O)8'(¥) —F'(O)8(y). 


One next writes v as a Fourier integral 
U(x, y, 2) = | | eltar+Pa(y, a, B) dadB. (20) 


Then, by substitution in (19), 


where all primes signify differentiation with respect to y. 

Apart from the delta functions on the right-hand side, representing the 
effect of the source, equation (21) is the same as Squire’s equation (2), 
with c and v put equal to 0. ‘The fact that « and f appear only in the 
combination 
at + fF = he (22) 


is again important. It follows from it that v itself depends only on y and on 


x? + 2% = 5°, (23) 
For if we substitute 


z= kcos8, B =ksin6, x = scosd, z= ssinA (24) 
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in (20), then, since dxdB = k dkd@, we obtain 


“of “Le 
¢=xi gft| o's 


“0 


= 27 | wkd (ks) dk. (25 
Thus v, a function of y and s alone, is the ‘Hankel transform’ (25) of vo, 
a function of y and alone. 


This v is to be obtained from (21), which states that for y > 0 and 
y < 0 the simple equation 
Un — Up (2 + T ] 218 (26) 
\ / 


is satisfied by vp, but that vw) and uv, are discontinuous at y = 0, in such 
a way that 
m V'(0) 


=e (1) a 


m 


vo( +0)—vo(—9) = 43° v,( +0)—~v( —0) = 





The problem has thus been reduced to terms not involving any generalized 
functions. 
SIMPLE SOLUTIONS FOR TWO PARTICULAR CASES 


There are two particular, although somewhat artificial, distributions V’(y) 
of the undisturbed velocity for which very simple solutions exist. First, if 


V(y) = U+Ay, (28) 
so that the oncoming flow is uniformly sheared, then (26) becomes 
v, —k*v, = 0, (29) 
and the solution of (29) which tends to zero as |y — © and satisfies (27) is 
m m V'(O)_, 
gq. — ——_g-Klyle poh 2 p—kly| 2() 
q 5 sen wae mas , 3 
0 Bo BNL BR V(0) is 
where sgny is +1 when y >0 and —1 when y <0. _ By the result 
(Watson 1944, § 13.2 (1)) 
| | e-klul J (ks) dk = a (31) 
Jo ? V (y? +5") 
and its derivative with respect to y, equations (30) and (25) give 
0 m V'(O) m 
(= == ce ee, (32) 
oy 4ar V(0) 4ar 


> 


where r = \/(y? +58") = 1/(x?+y?+ 27). That the coefficient of V’(0)/V(0) 
must be simply the potential of a source of strength m is independently 
obvious from equation (19) when the second term therein vanishes; 
similarly, the other term must be the y-derivative of the source potential 


(—m/4nr). 
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The first term in (32) is the ‘ primary’ flow due to the source, and the 
second term is the small-disturbance approximation to the simple-shear 
secondary flow (Lighthill 1957 a). ‘The full secondary flow in the simple- 
shear case (28) is the solution obtained by assuming not small source strength 
but small shear; it is a perturbation on the irrotational flow due to the 
source, derived by neglecting the square of the shear; at a large distance 
from the source, where the disturbance to the oncoming stream is also 
small, it takes the form shown in (32). What is interesting in (32) is that 
no higher powers of the shear (representing tertiary, quartary flows, etc.) 
appear, although in the present analysis these have not been neglected. 
It is, as we shall see, only for the particular, and artificial, undisturbed 
velocity distribution (28) that this is so. 

Moreover, such terms are not absent in w. By (32) and (15), with 


V = U+Ay, 
| Y a 0 A ms me 
77H) Cans + aes) 9 


which after an integration by parts gives 
ms A? mz 
* far ~ TUF Ay) Fer’ 
of which the last term is a tertiary-flow term. More seriously, this term 
makes the integration of the equation of continuity to obtain w impossible 
without allowing u to become logarithmically infinite at both v= + x 
and x= —«. ‘This is an extremely discouraging feature of the present 
solution. However, it is related to the circumstance that V(y) actually 
becomes zero for finite y, and we shall see that for more realistic velocity 





distributions the difficulty does not arise. 
‘There is one other V(y), not specially realistic, but at least not leading 
to unrealistic infinities, for which a simple solution is possible. This is 
V(y) = Ue”. (35) 
Again, the solution is most readily obtained from equation (19), which in 
this case becomes 


Vv—v = e + 1) ma(x)_(0)8(2), (36) 


whence é me~’" io 
od = (= +a) - 4 | (37) 
oY TY 


the term in square brackets being the fundamental solution of V*v = A 
Alternatively, (37) can be obtained by first solving (26) and then making 
use of the result (Watson 1944, § 13.47 (2) with » = 0, v = $) 


-% —i v(y*+s?) 
k sda es ‘ 


ol) VE + 48) | ¢ = : 38 
| —__.___— ¢ J (ks) dk T/C) (38) 


9 5) 
Ji (Rk? + A?) 
The modification of the primary- and secondary-flow terms in (37) by | 
the exponential factor inside the square bracket means that tertiary- flow 
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terms in A? and higher-order terms are present in v in this case, but that they 
all add up to a solution tending to zero rapidly as r+ oo. As to w, 
equation (15) with (37) gives simply 


C ( me~“" (39 
Ww —| — ——}, 3 
; oP 4ur Ss 


and for u the equation of continuity yields after some reduction 


“( me~*" mX/ e ’ a ae ie 40) 
u = BA dar) * Sahay +) | <<. — 


Thus there is an exponential falling-off as r— oo in wu also, except as 
x -> + 0 for values of y and z which are not large. This last result is 
realistic, and corresponds to the fact that fluid which has passed close to 
the source has been permanently affected, because its vortex lines have been 
deformed. In fact, 


; : mA/ e J env (a? + ut +2 





4m\cv V/ (x* + y* + 3°) 
mXr/ o Woes , : 
7 (= 3 \)2Ko(M (9° + 2")5, (+1) 


by Watson (1944), $6.22 (5). All the outflow due to the source has somehow 
been channelled into the region immediately downstream of it*. ‘That the 
total volume flow in the x-direction in this region is in fact m is easily verified 
from (41); we obtain 


| (lim uw) dydz =m | _ tK,(t) dt = m, (42) 
by Watson (1944), § 13.21 (8). The solution, however, is not worth pursuing 
further because of the impossibility that in a real flow the velocity distribution 
could be well represented by (35) over any large range of y. 


4. METHOD OF SOLUTION IN THE GENERAL CASE 
We pass to methods applicable for a general distribution of velocity V(y), 
subject only to the restriction that the limits V(o) and V(— 0) exist, and 
that V(y) is nowhere zero, even at these limits. (We shall also require V(y) 
to tend to these limits sufficiently fast for certain integrals to converge.) 


Thus, the parallel flow is supposed to become uniform as y - + « and as 
y — — ©, as in a wake or in a mixing region between two streams. 
The equation (26) as y > + tends to the simple form (29) and has 


solutions e-*” and e*”. Only the former tends to zero and so is of interest 
to us. We therefore introduce the notation v,(y,k) for the solution of (26) 
which is asymptotic to e-*” as y+ +. Similarly, v.(¥,k) signifies the 

* We shall see that this is due to the displacement of streamlines, as they pass 
the source, towards the region of lower undisturbed velocity. 
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solution which is asymptotic to e*” as y+ —o. ‘The necessary dis- 
continuous solution v) must take the form 
Uy = A(k)v,(y,k) (y > 9), B(k)vy,k) (vy < 9), (43) 


for some constants A and B. But, by (27), 











A(k)v,(0, k) — B(k)v2(0, k) = ae (44) 
; , V'(0 
A(k)v'(0, k) — B(R)2'(0, 8) = Ti (45) 
and the solution of these simultaneous equations is 
m (ies V'(0) 
ACR) = —~< v.(0,k) — ——~v,(0, k) > 
Ath) = Tap 00) — Hoy 7204) p 
fri : (46) 
i) = a ey 
Hk) = teW(h) |" (9, &) — Toy @6 , Dp> 
where W(k) = v,(y, k)vs(y, k) — vay, 2)v1(y, &) (47) 


is independent of y by a well-known property of second-order linear 
differential equations (and it is actually (47) for v = 0 that has been used 
to derive (46)). 

Equation (46) would be useless if the Wronskian W(k) vanished. 
However, this cannot happen. For it would mean that v, and vy were 
simply proportional to one another, so that a continuous (and continuously 
differentiable) solution vp existed tending to zero at y = + 0. Physically, 
a stationary wave could then be present without the need for any extraneous 
disturbance like the source here discussed. ‘This state of affairs is not at 
all plausible; but the possibility can be dismissed altogether on mathematical 
grdunds, since if such a v,» existed, it may be shown from the differential 
equation (26) that 


o> (f, Vv" \2 ) 
| } (5 Tv to) + Reo iy = 0, (48) 


so that vg must vanish identically (because otherwise the left-hand side 
would be positive). 

The character of the solution v,(y,) (and its sister solution v.(y, k)), 
on which all has been made to depend, will be sufficiently uncovered in 
what follows by three theories, one for large k, one for small k, and one 
for general k subject to the restriction that the total variation of V(y) is 
small compared with its absolute magnitude. In each case the corresponding 
properties of u(y, k), v(y,s), u and w are inferred from the above formulas 
and those of § 2. 


5. "THEORY FOR SMALL 7 
By properties of the Hankel transformation (25) we shall be able to 
deduce the behaviour of v for small r from the behaviour of vy for large k. 
To study the latter we solve equation (26) in the form 


y" 


a)” ery: eer ee ats 
%— kv = yp wo 
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by ‘variation of parameters’. Picking the solution v, asymptotic to e~*¥ 
as y > + 0%, we obtain 





tis 1 _ p-2k(q—-y) y" 
v4(y,k) = eS 1+ | : (q) 


{ } _— V(q) v1(9, k)eka dq >. (49) 


This integral equation is useful for many purposes, in particular for proving 
by Picard’s well-known method the existence and uniqueness of v,(y, A), 
and its analyticity as a function of k for #{k} > 0. It also tells us the 
expansion of v,(y,k) for large k, giving to a first approximation 
v,e*Y = 14+ O(k-") and to a second 


( ee l 
a = el 1+ = | Dy, + (js) rs (50) 


2k}, V(q) k? 

Similarly aut 1” Vv") 1\) 7 

ee Uo = cL < + — ———— 6 a > pe 
g = ebs 1 i | , Va) dq+O B) (51) 

and so the Wronskian of v, and v, is 

‘ ; ro I "(q) 1 ” 
(Rk) = 01 Vg—VgV, = 2R+ =~ dq+ Of 5 2 
W (Rk) = U1 Vo 9 V; i; | j (9) dq of; (5 ) 


its independence of y being a check on the work. Hence, by (46), 


( oo 1\)\-1 
A = al 2k+ | 79) dq+ o(;) r x 
77 | Geer q J 








47” | V(q) Rk 
(, 17° V%q) v0) /1\) 
Lea) Taye By +R) 
_ mf 1 V"(q) 1 V’(0) I bl 
"as" 2d, Vg) 2 k V0) (js) , (7) 


and similarly 
m | lf @ 1 V’(0) I : 
B oS 8x2) —] T Th | a dq meer V0) T O —— iy (54) 
whence by (43), (50) and (51) 
m oa 1 py V"(q) 1 V’(O) | 2 = 
= gare oe) ve (i) Ph 
The author has evaluated also the coefficient of k~* in the expansion of 
u(y, k) for large k, but the only property of this coefficient that tells us 
anything about the behaviour of wv for small r is that it is O(\y|) as |y| +0; 
the next term still is one in k-°. Hence, writing F(y,k) for the error term 
in (55), we have by (25) and (31) 


my [V(0) 1 (V%q),)\m 





o = Fs 1 Voy * 2), Vey tae? 





~ 4ar 
we ini S 
e ye é J o( ks) E (v, Rk) dk, (56) 
with Cily C, 
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for some C, and C,. The integral in (56) is therefore bounded as r > 0; 
it is less in modulus than 


0 e-k\u\ dk > dk 
e +Cy R 





Cily anise 


Yi), THk |, Waa 
which is bounded, because the first integral is O(log|y|-'). (Actually, the 
integral in (56) tends to a finite limit | . F(0,k) dk asr +0. However, this 


limit cannot be evaluated unless v)(0,k) is known for all k.) Part of the 
middle term in (56) is also bounded, and so finally 
my V'(O) m 

= far ~ V0) Far 

Equation (57) shows that the part of v which tends to infinity as r > 0 
consists simply of the primary flow and the simple-shear secondary flow. 
In other words, near the source the variations of velocity and of shear are 
small enough for the approximations of neglecting the square of the shear 
and regarding it as uniform to give the right answer. 

It follows from (57) and (15) that 


(- gs) +00); (58) 


this also is the sum of the primary flow and the simple-shear secondary 
flow, since the latter has no s-component if the squares of the disturbances 
are neglected. Finally, the equation of continuity gives 


0 m mV'(O) y . 

u -x(-g)-= V0) ata(1+ *) +0(), (59) 
the second term being obtained by integrating y/r? with respect to x from 
the lower limit — ©; again, this second term agrees with the small- 
disturbance approximation to the simple-shear secondary flow (Lighthill 
1957 a). 

The work of this section gives a check that fluid is truly emerging from 
the source at the rate m units of volume per unit time, as well as showing up 
the character of the simple-shear secondary flow solution as a valid second 
approximation in the neighbourhood of the source. 


+O(1) asr>0. (57) 


it | 
A) 
= a 


O8 


6. EVALUATION OF vy FOR SMALL k 
By a process inverse to that of §5, the behaviour of v for large r can be 
deduced from that of v, for small k, which is obtained as follows. 
A uniformly valid first approximation to the solution v,(y, k) for small k 





v, = VOY) ey, (60) 


To see this, note that when y is large V(y)—> V( co), so (60) tends to e-*Y 
as it should, but when y is not large e~*” = 1 (for small k), and so (60) is 
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approximately proportional to V(y), which in turn approximately satisfies (26) 
for small k. These facts suggest the transformation 








V(y) 
= - p— ky 61 
a KU Y,, ) 
ms) V ( 10) 1 ( ) 
which on substitution in (26) gives the equation 
d du 
—( Ye-2ku —*) = RY'e ky, 62) 
ral dy ! ee) 
where for shortness we have written V2(y) as Y(y). From (62) and the 
condition u, > 1 as y > © we obtain the integral equation 
mo dq oe ry 2k ] 5 
u(y)=14+k}| zal] Y'(Deheu,()) dl, (63) 
Jy ¥(Q) 24 


whence successive approximations to uw, for small k may be found, the first 
being uv, = 1+O(R), the second 
= ¥(«)—¥(q) | 
= t k Se ee — C 2 6 
u,=1+k | ; Yq) dq + O(R?), (64) 


and the third 

















ery Y(«)— Y(q) " dq 
u, = 1+k . =a ee | ; Y(q) * 
| Ly ( x - } ih (q) a (/)} di + O(R?). (65) 
Similarly, we put ae we) okvy (66) 
“Ta 
and find that 
af ete 
uy =1+k | V9) dq 
3” dq pe? {¥(— ©)-YOHV(Q)+ YO} . i 
_ Re +O(k). (67 
k | Ya) | . Y(I) dl+ O(k*). (67) 





Y¥(2)- VO), B/? W(2)- YOHYO)+ VO} ay, Ocaay 


(68) 











i Y(—o)-Y(y)_  F ’ fY(— 0)— Y(q)}H{Y(y)+ Y(q)} . 
_ Y(y) ¥(y) J -« V(q) . 

x dq+O(k*), (69) 
, ¥ ) , , 
_= WW ow (uy Uy — Ug U, — 2Ruy Ug) 
k{Y(00)+ Y(— o)} k 





(— 00) — Y(q)} dq +O(k3), (70) 
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where on collecting terms it was found that all those involving y explicitly 

disappear, which is a check on the result. 

Now, by (46) and (66), 
m V(0) 

—— Q,R ku.(O, R)} 

| We VW =o) 200 »k) + ku,(0, k)} 
m V(0) ‘ Y( - 0) = 
421V(k) V(— x)|* Y(0)  Y(0) ~ 


<}_—#)=Y@} ag +00") |, 71 








and similarly 


' m V (0) Y co nT eee ase al 
i= Ph) Vx | —k \ a + : | {Y(c«)— Y(q)} dq+ O(n) 
(72) 


results which with equation (70) give non-zero limits for A and Bas k > 0) 
and give also their derivatives with respect to k atk = 0. 
By (43), (61), (65), (70) and (71), for y > 0, 
: mV(y)e*" Y(-— a) 
“0 = “F2V(0) ¥(«)+¥(— @) 


E af Sega 




















Y(— o) 





{¥ (x) — V(@)}{¥(— ©) — ¥(q)} oer 
“ — k | Rats Wis pi ore) dg + O(k )} (73) 
and similarly for y () 

- mV (y)er" Y( a a 

“0 47°V(0) Y(«~)+Y(— x)” 











Now, equations (73) and (74) are to be used differently in different 
regions. For large |y|, they are best used as they stand, or at most after 
such simplification as is gained from replacing V(y) by V(0o) or V(— ), 
as the case may be, and suppressing the first integral inside each square 
bracket. However, for values of |y| which are not large, it is permissible, 
and valuable, to expand the e~*!¥' factor outside the bracket in each 
equation and take it inside. When this is done, both (73) and (74) can 
be rearranged into a common form in which the coefficient of k is seen 
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to be a continuously differentiable function of y even at y = 0. This form 
for |y| not large is 


mV(y) Y(0o)Y(— 0) | sen y h " dg 


+¥(— x) LV(wseny) *), ¥(@) * 


0° 40) Y(O)+ VC &) 





@ 








The term independent of k in (75) is still discontinuous at y = 0, but this 
will be found to produce no discontinuity in wv itself, which is influenced 
only by the term proportional to R. 


7. "THEORY FOR LARGE r 
We now show that the properties of 
v=27 | wo(v, Rk)RSo(Rs) dk (25 bis) 
for larger = \/(y” +s") can be deduced from the equations (73), (74) and (75) 
which give the behaviour of vo(v,k) for small R. 
The case when |y is itself large is especially easy and may as well be 
treated first. For large positive y, (73) shows that 


Mm, 








Vy ~ a3 ku(1—ke,) (76) 
as k + 0, where - V(«) 2V*(— a) eS (77) 
1 V(0) V%(0)+V%— 0)’ ais 
and 
+ F-@)- Fe q+ | (¥(20)— Y(g)}t¥(- &) MDS 4 
= J_. 7 Veo) » (a+ ¥(-a)¥@) 
(78) 


The ‘Borel’ theory of integrals with exponential factors whose exponent 
consists of the variable of integration multiplied by a large parameter 
(here y) tells us that as y > x 


ea A9 


a 
tz}, 


¢ 


Qj | 
m,e-*¥(1 —kce,)J (ks) dk = — (m, a + mM, Cy =a Nar (79) 


v,~ 
where (31) has been used. The behaviour of w follows at once from (15) 
with the lower limit replaced by + 0, in which V can be regarded as 
constant on both sides, giving 





0 0? \ ] 
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The equation of continuity then gives a similar result for u, and the three 
results may be written 


m f m 
(u,v, w) ~ grad( — ie) +¢, - grad( ~ z). : (81) 


Equation (81) shows that, for large positive y, the effect of the source 
of strength m at the origin, in the midst of the sheared flow, is the same 
as if a source of different strength m,, and a doublet of strength m,c, with 
its axis in the y-direction, were present at the origin, but immersed in a 
uniform flow with velocity V(0o). (These are the two leading terms in 
the expansion of the irrotational disturbance above the shear layer in terms 
of singularities at the origin.) 

Alternatively, we can say that the effect is the same as if a simple source 
of strength m, were present, not at the origin, but at the displaced position 
(0, —¢,0). 

Similarly, for large negative y, 








(u, v,w) ~ grad( — a) —Cy 7 grad( 7) (82) 
where _ V(- @) 2 «) os 
m= V0) V{a)+VX{— a)” se 

and 
__ [? ¥(%)-¥Qq)  (¥(~)— Vigyv(— ©) — Yq); 4 
o- | ag th eee oe 


Thus, for large negative y, the disturbance is that which would be produced 
in a uniform flow by a source of yet different strength m, at the origin, and 
a doublet of strength m,c, with its axis in the negative y-direction. Again, 
it can if preferred be regarded as equivalent to a source of strength mz at 
the displaced position (0, +c, 0). 

The strengths m, and m, of the equivalent sources can alternatively be 
deduced by an image approach due to Mr M. B. Glauert. ‘To obtain the 
asymptotic disturbance at large distances, we replace the shear layer by a 
concentrated vortex sheet. Now, suppose first that the source of strength m 
is ‘above’ the vortex sheet, in the stream with velocity V(«). By applying 
the condition of continuous pressure and flow direction across the vortex 
sheet, it may be shown that if squares of disturbances are neglected the flow 
in the upper region is equivalent to that due to the original source of strength 
m and a source of strength m; (7 for image) placed at its mirror-image in the 
vortex sheet; further, the flow in the lower region is that due to a source 
of different strength m, (¢ for transmitted) at the position of the original 
source. Here, 

7% — 00) — V%( 00) 2V(w)V(— 0) 


ni. = Na mM, = m (85) 


+= Pw) F VHa V*(0)+ V{— 00)" 








When the original source approaches close to the vortex sheet it almost 
coincides with its image, so effectively in the upper region we have a source 
of strength m+m,, and in the lower region one of strength m, Finally, 
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when the source is in the shear layer itself, at y = 0, the fluid produced by 
it fills a cylinder (stretching back from the source) of cross-sectional area 
m/V(0); this measures how much the rest of the flow is pushed out. Its 
effect is the same as if a source of strength mV()/V(0) were present just 
above the shear layer where the velocity is V(«). ‘Therefore, the effective 
source strengths m, and m, in the regions above and below the shear layer 
respectively should be 

(00) V(«) 


V 
m, = (m+m;) VO’ My = My Toy (86) 
which agree with (77) and (83). 

The effective displacements c, and c, of the source (equations (78) 
and (84)) can probably not be calculated by such a method, as they are of the 
order of the thickness of the shear layer itself. 

The limiting case when V(— «)/V(co) + 0 is of interest; then both 
m, and m, tend to 0. This is connected with the fact that m; > —m 
(the image of a source in a ‘free streamline’ separating a flow from a 
dead-air region being an equal sink), so that the source and its image 
cancel. ‘They leave as the leading term in the disturbance in the upper 
region of uniform flow, a doublet whose strength 


-() 


m 
m,C,—-> — V(0)V(~) | V(q) dq. (87) 
The case V(y) = Ue” studied in § 3 may be compared with this, since 
that distribution has in a sense V(— «)/V(oo) = 0. But for large r (except 
where y and 2 are not large and x is positive) the disturbance in that case 
was found to fall off like e~*’, so that not only the effective source but also 
the effective dipole and higher multipole strengths are zero; possibly the 
V( co) in the denominator in (87) might have led one to expect this. 
Note that the net volume flow through the part of the surface of a large 
sphere, with centre at the origin, which lies outside the shear layer, is 
V(a«)V(— w2){V(o)+V(— w)} 
3(M,+m,) = V0) Vac) + V{— co) m. (83) 
( ~ a 
This differs from m, being less than m for a mixing region and greater for 
a wake. ‘The balance is made up by flow across the part of the surface of 
the sphere immediately downstream of the source, where streamlines have 
been displaced. In a mixing region, for example, the displacement is 
towards the region of lower velocities, and so the velocity at any fixed point 
in space is increased. ‘This point will be returned to in §8. 
We begin now to investigate the flow for large r within the shear layer, 
that is, where |y| is not large and therefore equation (75) holds. ‘This work 
depends on the result that 





ro F(0) 3F"0) 3.5 F%(0) 
| ; F(k)kJ (ks) dk ~ — o ta-"" "oa 

ass — oo, if F(R) is a bounded analytic function of k forO <k < wo. The 
author has been unable to find a clear-cut reference for this result, though 
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it must presumably be well-known to people who frequently use the Hankel 
transform. Itcan be approached either by real- or complex-variable methods, 
the latter being simpler, though they assume rather more. 

If, in fact, F(k) is regular in a wedge |argk a, then we write 


Kol tks) — Ky tks) 


J (ks) me 


(90) 


(by Watson 1944, $3.7 (8) and $3.62 (5)), and divide the integral (89) into 
two parts, shifting the path of integration to argk +a for the part 
involving A,(—cks) and to argk x for the part involving A,(zks). 


‘The integral becomes 


1 wet? | pve? aa 
| F(k)RK,( iks) dk - — | F(k)RkK, (tks) dk 


= : F(it)itK (ts) dt 4 q | F(—it)(—it)Ky(ts) dt. (91) 


TJ La |) 


Since A,(ts) falls off exponentially along each path in (91), the integrals 
may be asymptotically expanded, by a method directly analogous to the 
‘Borel’ method used earlier in this section, as 


PH ies) dt +- —— >» 


n! Win n} 


a 


uv 


r rn +1 BYXn) ) -meN4 270) r\n +1 BY) 
l <: |’ ne | (—2) TO 


net( 4 — x) 


x | 11K, (ts) ar|. (92) 


The contour integrals in (92) are clearly identical (and their values are 
simple combinations of factorials; see Watson 1944, $13.21 (8)); hence 
the terms with n even in the sum vanish, and the rest take the form (89). 

We wish to apply (89) to the case F(k) = u(y, k), so the question of 
the analyticity as functions of k of the expressions (43) for v» arises. ‘The 
analycitity of v,(v,) and v,(¥, k) and their derivatives with respect to y for 
Rk (is easily proved by Picard’s method from the integral equation (49). 
Hence vy will be analytic provided that W(k) is non-zero. But the proof 
based on (48) that this is so can easily be extended to complex k with non- 
negative real part; in fact, if W(k) vanished, so that a continuously 
differentiable solution of (26) tending to zero as y + + © existed, then the 
integral (48) would still be zero if k® were replaced by its real part and all 
other squares replaced by squares of moduli, and the conclusion that vp 


2) 


is identically zero would follow as before provided that #{k?} > 0. Hence 


we may take « = !7 in the above discussion. 
“ 4 
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Equation (89), whose applicability to the case F(k) = u(y, k) has thus 
been established, now gives with (25) and (75) that 


mV(y) ¥( x )¥( — ©) ] (" dq ] 


2nV(0) ¥(~)+¥(—- o)L J, Vg) Y(~)+¥(— &) * 
( (2 YXoo)—¥%q), = ¥%X— o)— ¥%q) 
| Seong a- | - dq\ | 33 


| (5), (93) 
as § = (x?+ 37)? —. of for values of y which are not large. 

It is, perhaps, tiresome that the expression for v when r is large should 
take three different forms (79), (82) and (93), according as y is large and 
positive, large and negative, or not large at all, and it is natural to ask 
whether they can all be replaced by a single, uniformly valid asymptotic 
form. ‘This does not seem to be possible. By careful inspection of (79) 
and (93) one can show that both are included in the more general expression 


poo a V1 om f* ¥(@)-¥@) 1) of 
7 V( = ) mM, oy FM, Cy 53 iar — in | =a C q a =) 


(94) 
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which therefore is valid both when y is large and positive and when |y| is 
not large; but for large negative y (94) tends to 
Ma _ Mey (95) 
4rr®— 4 : 
whereas by (82) there should be an additional term 3m, cy y"/4zr° for 
large negative y. 

However, such an expression as (94) is just what one needs to integrate 
(15) (with lower limit + ©) to obtain an expression for w when s but not |y 
is large. ‘This gives simply 

ms amy oe + 5) (96) 
~~ V(y)VO)Y(c0) + Y(— «)} 488 ~~ \ st 7? 


whence, by (93) and the equation of continuity, 








; 2Y(«)Y(— o) mx mV'(y) Y(a)Y(-« 
“ = ViyyVONY(0)+Y(— ©) 428% 27 V(0) V(x) + Y(— ov) * 
‘UY dg I ro Vato) Fg) 
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LJ.Y@ Varroa, Wore 
“0 ¥?2(— c)— Ye )J1+-4x/s I - 
= a Eo te OF 

J» Y(-—o)Y(q) 2° s 
The leading terms of (96) and (97) show how the effective source strength 
varies inversely as the undisturbed flow velocity V(y) within the shear 
layer, tending to m, (see (77))as y > + and to mz, (see (83))as vy > — @. 
This is connected with the fact that the pressure is constant across the 
shear layer far from the source, and that the pressure varies like —pV(y)u 
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along any streamline. ‘The additional term in uw, not inversely proportional 


to V(y), is due to the displacement of the streamlines. 


8. ‘THEORY FOR SHEAR LAYERS WITH SMALL TOTAL VARIATION OF VELOCITY 


In §4to§7 the theory has been developed for a general shear layer, and 
asymptotic expressions for the velocity field found, both as r — 0 (in the 


form of terms of orders r~* and r plus an error of order 1), and as r 


~ L 


(in the form of terms of order r* and r* plus an error of order r~*). 
However, the detailed characteristics of the join between these solutions 
remain unknown, and it was to give some understanding of these that the 
solution to be described in this section was sought. ‘This solution is valid 
for all values of x, y and z, but only for a restricted class of shear layers, 
namely those within which the total variation of the undisturbed velocity V’(y) 
is small compared with V’(y) itself. If the minimum value of V(y) is (1 —e) 
times the maximum, we seek an approximate solution by consistently 


neglecting €?. 


As before, we begin by obtaining v,(v,k). ‘The integral equation (49) 


shows that, if the first power of « be neglected, v, is simply e-*”. 


‘This 


first approximation is valid uniformly, even for small k, since the fraction 
involving k within the integral is always less than (z—y); and indeed the 
solution (60), valid in the limit k — 0, differs from e-*” only by a term of 


order e. 
Applying equation (49) again, we now see that 


S| (Ee 2h(q—y) V"(q) ) 
ah "4 dq ¢ + O(¢ ). 
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Equation{(98) is similar to the earlier approximation (51), except that the 
exponential factor cannot now be suppressed because k is not being 
assumed large. However, the V(q) in the denominator can be replaced 


by a mean value V,, 
V(y)) with an error O(e?). 
Similarly, 
v1 VQ) ng 
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whence it follows easily that the Wronskian 2, z. 





W(k) = 2k+ | = dq + O(c) = 2k-+ O(e). 
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Next, by (46), 
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Hence, by (43), (98), (99) and (102), if <? is neglected, 
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Equation (103) has been written in the form which gives clearest 
confirmation of the result proved in §5, that wv differs from the sum of the 
primary flow and the simple-shear secondary flow by an expression bounded 
asy>0(Q. For it gives the Hankel transform (25), after the substitution 
2q—y! = 21 in the last integral, and use of the result (31), as 
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(105) 

The terms in curly brackets in (105) may be called the ‘tertiary flow’ 
contribution to v; they are the limit as r — 0 of what remains when the 
primary-flow and the simple-shear secondary-flow terms are taken away. 
The presence of the integral shows that this tertiary flow depends on the 
value of V"(y) not only at the source itself but throughout the shear layer. 
Expression (104) is not, however, the most suitable form of v for showing 

its behaviour away from the source. For this, we first integrate (103) by 
parts, substitute 2g = / and only then take its Hankel transform, namely 


’ -o @ ae _ 
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are +r 4 jx 2V,,, {(y — D2 + 528? 

(Alternatively, on integrating (104) by parts and changing the variable of 
integration we get (106).) 

Equation (106) has a simple and valuable geometrical interpretation. 

It says that the y-component of disturbance velocity v is that due to the 

source of strength m at the origin, together with a line distribution of sources 

along certain stretches of the y-axis; to be precise, the stretches from 2y 

x “0 

— | is the same 
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negative. ‘The source strength at the point (0,/,0) is V’(3/)sgn //2V,,, per 
unit length. 
The values of w and u can be interpreted in the same way. Thus, by (15), 
C ( 1 Pw, nes | 
w= Pal Vy) | V(q)v(x, g, 3) dq f 
€ 


¥ » V-Vy) _mg , 
| u(x, q, 2) dq+ | V te aay (107) 


L “ m 
if O(c?) be neglected, and when the integration is carried out we find that 


ms mf ° “0 V3) z ' 
e-pnta(|,-| ja oapeee% (108) 


which is the z-velocity due to the distribution of sources described above. 
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Also, on integrating 
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from — © to x to obtain u, we clearly get the x-velocity due to the same 
distribution of sources, except from one extra term in dv/éy resulting from 
the dependence on y of one of the limits of integration in (106). ‘This 
additional term on the right of (109) makes a contribution 
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to u, whence (by evaluating this integral) 
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Thus, apart from the last term in (111), the complete system of 
disturbance velocities can be identified with the flow due to the source 
distribution along the y-axis described above. 

This does not mean that it is irrotational. ‘The velocity at a point (x, y, 2) 
is that due, not to a fixed system of sources, but to one which itself varies 
with y. ‘This gives additional terms in the y-derivatives of uw and w which 
(together with the y- and z-derivatives of the extra term in uw) prevent any 
component of the vorticity from vanishing. An actual expression for the 
vorticity is obtained later (equation (117), in which by (120) ¢ must be 
taken as —m/4nrV,,). 

Mr M. B. Glauert has shown the author an ingenious derivation of the 
equivalent source distribution from his image theory mentioned in $7. 
He regards the shear layer as made up of a very large number of elementary 
vortex sheets. he velocity at a point (x,y,z) due to the source at the 
origin is then represented as the velocity field of the original source itself, 
together with its images in all vortex sheets which do not pass between 
the origin and (x,y,z). ‘This is permissible if O(e?) be neglected because 
(i) the ‘transmitted’ source strength m, by which the original source 
strength m needs to be replaced as a result of vortex sheets between the 
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origin and (x,y,z), is equal to m+O(e?) (as (85) shows) because the 
difference of velocities across all these vortex sheets is O(e), and (ii) the 
images of the image sources will have strength O(e?). 

It follows that there will be image sources from 2y to © and from — « 
to 0 if y > 0, and from — @ to 2y and from 0 to « if y < 0. Further, 
the source strength between (0,/,0) and (0,/+d/,0) due to the vortex 
sheet lying between (0, 3/,0) and (0,3/+3d/,0), across which the velocity 
change is V’'(4/)} dl, is seen from (85) to be {V’(3/)sgn//2V,,} dl. The 
source distribution is therefore exactly that obtained analytically. 

To explain the additional term in (111), note that Glauert’s vortex 
sheets are displaced by the source flow. ‘To a first approximation (neglecting 
O(e)) the y-component of displacement is 


1 rz my my(1 Tx r) (1 1 2) 


The disturbance velocity u will therefore have an additional component 
approximately 
ee. my(1 + x/r) Vy) = m V'(y) y(1+x/r) 
lr 4nV,,(y? +2?) f i eae 4n Vi yet?’ 
as in (111), with error O(e?). 

The equivalent source distributions appropriate to a point with y > 0 
and to one with y < 0 are illustrated for a mixing-region flow in figure 1. 
A general picture of the disturbance flows can be readily got from this figure. 
If y were actually greater than the greatest value it takes in the shear layer, 








(113) 


the upper part of the equivalent source distribution would disappear 
altogether. In this region, therefore, the flow is the irrotational flow due 
to a fixed system of sources. ‘Their total strength is easily verified to be m,, 
and their total dipole moment m,c, (see §7), to within an error O(e?). 
Similar results hold for values of y less than the least value it takes in the 
shear layer. Within the shear layers, the sources and sinks are so arranged 
as to produce a pronounced ‘downwash’ in the case of a mixing region, 
but this is somewhat mitigated in the case of a wake. 
Far downstream of the source, where x is large and positive but y and 2 
are not large, the disturbance velocities are dominated by the last term in (111) 
(interpreted in (113) as a displacement term); thus, 
yt , , 
ae ep ed eS) 5, 0, OF. (114) 
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We may use (114) to check conservation of volume flow. ‘The net volume 
flow, across that part of a large sphere (with centre the origin) which 
intersects the portion of shear layer downstream of the source where (114) 
holds, is 
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where (88) has been used with O(e?) neglected. ‘Together with the volume 
flow $m, across the part of the sphere above the shear layer, and 5m, across 
the part below, this makes up exactly the rate m at which fluid is being 
produced at the source. 

We here have exemplified from the case of small velocity spread how 
displacement of streamlines in the shear layer can remedy deficiencies 
between the rate of fluid output at the source and the rate of escape to regions 
far from the shear layer. In layers with large velocity spread the nature of 
the displacement is greatly altered (it includes secondary downwash and 
the like as well as the primary flow due to the source) but it must still produce 
this same effect. Of course, if terms of the order of the square of the source 
strength were included, the secondary trailing vorticity would be present, 
and then the streamlines far downstream would spiral about their mean 
positions; however, the conclusion would be unaltered. 
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Figure 1. The equivalent source distribution, representing the effect that a weak 
source in a mixing region with small velocity spread makes at the point marked 
with a circle (which has y > 0 for the left-hand of the two distributions shown, 
and vy < O for the right-hand one). The heavy plus sign represents the 
original source in each case, and the light plusses and minuses additional 
sources and sinks respectively. ‘These can be regarded (Glauert) as images 
of the original source in the component vortex sheets of the shear layer. 


9. RELATION BETWEEN THE PRESENT SOLUTION AND THE EXACT-PROFILE 
SECONDARY FLOW 


In the present paper the disturbance due to a weak source in a shear 
layer has been investigated by neglecting the square of the disturbance 
due to the presence of the source, no other approximation being made 
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in §2 to§7. In §8 the additional approximation of neglecting the square 
of € (the relative spread of velocity in the shear layer) was made. 

It is interesting to compare the results of the double approximation 
with those of another procedure, in which first the effect of a small non- 
uniformity in the oncoming stream on a basic irrotational flow (such as 
that due to the source) is considered by neglecting squares of e, and the 
result is then further approximated by assuming small disturbances. ‘This 
procedure was expounded briefly for the case of a general parallel flow in 
Lighthill (1956, p. 39), although the rest of that paper was devoted to the 
case of a uniformly sheared upstream flow. 

‘The term ‘secondary flow’ may appropriately be used to describe any 
flow field computed by taking an irrotational flow with uniform upstream 
velocity as the first approximation (or primary flow), and going to a second 
approximation by allowing a small spread in the upstream velocity but 
neglecting its square. ‘The secondary flow may be called a simple-shear 
secondary flow if the upstream flow is taken as uniformly sheared, and an 
exact-profile secondary flow if the true upstream velocity profile is used. 

As shown in Lighthill (1956), secondary flows can be calculated by 
finding how the vortex lines of the upstream parallel flow become stretched 
and rotated by the primary flow. ‘The secondary vorticity field having been 
found, some fairly arduous integrations are needed in general to determine 
the secondary velocities; but such integrations have been carried out for 
the simple-shear secondary flow about a sphere in Lighthill (1957 b). 

The secondary vorticity field for any upstream velocity profile can be 
expressed in terms of the ‘drift function’ ¢ for the primary flow. ‘This 
function is such that planes of fluid initially at right angles to the stream 
are distorted by the primary flow into shapes given by equations? = constant ; 
t represents the time at which a given fluid particle reaches a point, measured 
from when it would have reached x = 0 had the uniform upstream flow 
continued undisturbed. By considering how this distortion of planes of 
fluids affects the vortex lines of the secondary flow, which initially lie in 
such planes, one obtains the fairly simple equations (30) of Lighthill (1956) 
for the exact-profile secondary vorticity field. 

These can be further simplified in the region away from the source of 
disturbance, where disturbances are themselves small, since in this region 
the drift function takes the simple form 

im ets (116) 
where ¢ is a ‘ disturbance potential’, defined so that the full velocity potential 
of the primary flow is U(x+¢). The form of the secondary vorticity 
distribution in this region is easily deduced (Lighthill 1956, (32)) in terms 
of ¢; in the special case when the upstream velocity is V(y) (independent 
of z) it becomes 


curl(u, v, w) ~ d —V"(y) z= — [ 
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The disturbance velocity field associated with this disturbance vorticity 
field is 

) 

\ 


(u, v,w) = 4 —V"(y) | & as, V'(y)¢, 0;—grad ub, (118) 


where Vs = V"(y), (119) 
since the curl of (118) is (117) and its divergence vanishes (as was secured 
by subtracting gradw). In some cases it is necessary to add on to (118) 
an irrotational part representing the far field of the ‘secondary trailing 
vorticity’ (horseshoe vortices wrapped round the disturbance); this part 
is given in equation (85) of Lighthill (1957 a). 
When the disturbance is that due to a source, we may take 
m 


$=-77- (120) 





The asymptotic form (118) of the secondary flow for large r then needs 
no addition due to the far field of the second trailing vorticity, because the 
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latter is of smaller order, O(r-*), as r — co. Hence in particular 


=, Cus 
v~ V(y)b— = (121) 


as r > 0. 

Now, the solution (121) for the small-disturbance form of the exact- 
profile secondary flow, which is its asymptotic form far from the source, 
can be identified with the solution of $8, where the small disturbances due 
to a source are approximated for an upstream flow with small velocity spread. 
For if O(e?) be neglected, equation (19) (except at the origin) becomes 

V"(y) my 
V%n = i pe 7" (122) 
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and it is a simple matter to verify that the solution (106) does indeed satisfy 
(122). Hence, if we were to write v in the form (121), where ¢ is given 
by (120), then ¢% must satisfy 


A A 


= (Vp) = V"(y)$ +20") 5 V8 (123) 
= V"(y)b+ VO) 5, = FV") (124) 


and integrating with respect to y we obtain equation (119) for %. Similarly, 
the identity of w and wu in the two solutions may be checked. 

However, though this shows that the solution of §8 could be obtained 
by the alternative method leading to (118) and (119), the method adopted 
in this paper seems preferable, because it gives a lot of information about 
the solution to the more general problem in which é? is not neglected, and 
because the direct solution of the Poisson equation (119) would in any case 
involve the same Fourier and other analysis which had to be introduced 


above. 
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10. RELATION TO THEORIES OF THE DISPLACEMENT OF THE STAGNATION 
STREAMLINE 
One interesting property of shear flow is that when a Pitot tube is placed 
in such a flow with the open end of the tube pointing upstream the pressure 
measured in the tube is the stagnation pressure not on the streamline 
approaching the tube directly along its axis but on one displaced by a 
measurable amount in the direction of higher velocities (Young & Maas 
1936). ‘This may be regarded as evidence of the presence of ‘downwash’ 
(flow at right angles to the undisturbed streamlines, in the direction of 
decreasing velocity) on the streamlines approaching the tube along its axis, 
and so it becomes a problem to compute this downwash from secondary-flow 
or other theories and calculate the displacement from it. This is the subject 
of a companion paper (Lighthill 1957 b), but there is one aspect of it which 
must be mentioned here, as it can be properly treated only by the present 
theory. ‘This aspect is the behaviour of the downwash far ahead of the tube. 
The ‘downwash function’ D(s) is the value of (—v) on the axis of the 
tube (y = zs = 0) at x = —s, and is expected to be positive if V’(0) > 0. 
The most important region of downwash is near the tube; here D(s) is 
best found from the simple-shear secondary flow, or from a higher approxi- 
mation in a sequence starting from this. But the downwash so calculated 
falls off (see (57)) only as 

V'(0) m 
D(s) ~=>~ = 
V(O) 4x8 


at points far enough ahead of the tube both for the disturbances to be small 
and for the tube to be representable by a source of strength 
m = 1nd? V(0), 


(126) 
where d, is its external diameter. But the displacement calculated from (125) 
would be logarithmically infinite, since (125) is not integrable up to s = . 
If tertiary-flow terms were included, they would be even greater for large s, 
and still worse infinities would arise in the displacement. 

The difficulty can be resolved by using the complete solution of this 
paper in the region where s is large enough for the disturbance to be small 
and for the tube to be representable by a source of strength (126). In this 
region the downwash falls off only initially like (125) (see (57)), and later 
like s~?; more precisely, by (93), we have 

m Y(0)Y(—o) f{ (* Y%(0o)— Y%(q) 
Ins (¥(o0) + ¥(— a) |) Y(o¥q) “!7 
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ass—>oo. The difficulty of D(s) not being integrable up to s = oo then 
disappears. 

The problem is thus treated by means of two overlapping solutions. 
For small to moderate r, the simple-shear secondary-flow theory is used; 
for moderate to large r the exact-profile small-disturbance theory is used. 
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Expression (125) is valid in the overlap range of ‘moderate’ r; for it 
represents the asymptotic behaviour of the simple-shear secondary-flow 
theory as r -> 00 and also that of the exact-profile small-disturbance theory 
as r >0. 

To apply this to the calculation of displacement, let sy be a value of s 
in this overlap region. Then the displacement of the dividing streamline 
due to downwash D(s) for s < s)is calculated from the simple-shear secondary 
flow. For s > s, the disturbances are small, and the contribution to the 
displacement from downwash in this region is 


1 rx 
V(0) |. D(s) ds, (128) 
where D(s) is given by the theory of this paper. 
Now, we can write the contribution (128) as 
a i V0) m 


provided that 
V'(0) m 
V(0) 42 
and in the limit on the right the downwash function D(s) is taken to be that 
given by the small-disturbance theory right down to sy = 0. It is because 
this D(s) is asymptotic to (125) as s + 0 that no substantial difference is made 
by allowing sy to become smaller and vanish in (130), since sy has already 
been assumed sufficiently small for the approximation (125) to be good. 
We call s,, defined by (130), the upstream cut-off of the simple-shear 
secondary flow. ‘The displacement for s > sg, as (129) shows, is equal to 
what it would be if the downwash were equal to the simple-shear secondary- 
flow value (125) up to s =s, and zero for s > s,—in other words, if the 
secondary flow were cut off at s = s,,. 
A precise expression for the cut-off value s, can be obtained if the theory 
of §8 for the case of a small velocity spread be used. Then, by (104), 
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whence, by (130), 

1 na 
200) 2 {-—V"()+V"(—D)jlog/ dl. (133) 

Since the denominator 2V’(0) in (133) is the value of the integral if the 
log/ in it be suppressed, we can say that logs, is the average value of log / 
in the shear layer, if the ‘weight’ used in averaging is —V"(/)sgn/. ‘This 
makes s, itself a ‘geometric’ mean of /in the shearing layer. (Here, / means 
the difference in y-coordinate between a general point and the source.) 
Thus, s,. is of the order of the layer width, as indeed is easy to see 
geometrically from the ‘equivalent source distribution’ described in § 8. 
In more general problems it is reasonable to suppose that s, is still of the 
order of the layer width. 

We may consider also the implications of the present theory for problems 
like the displacement effect of a sphere, as treated by Hall (1956) and by 
Lighthill (1957 b) on the assumption that the inviscid flow is a good enough 
approximation (in other words, that separation and wake flow do not 
influence the displacement). ‘The upstream effect of the sphere is that 
of a doublet of strength 





logs, = 


m = 27a®V(0), (134) 


if the sphere has radius a and centre the origin. With this value of m, 
the flow in the region of small disturbances (due to the doublet) is given 
simply by taking the x-derivative of the velocity field of this paper. 

The upstream cut-off of the simple-shear secondary flow is not crucial 
in this case, since the simple-shear secondary downwash is itself integrable 
up to s= oo. However, the cut-off reduces the displacement below that 
obtained by this integration, and the reduction can be calculated, as follows. 

If D,(s) is the downwash function due to the source studied hitherto 


in this paper, and D,(s) = —D{(s) is that due to the doublet, then the 
displacement for s > sy is 
l L 
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where (105) has been used to approximate to D,(sy). Now, the first term 
in curly brackets is the simple-shear secondary-flow value, and so the 
difference, between the true displacement and that given by assuming that 
the simple-shear secondary flow persists to infinity, is given approximately 
(using the value (134) for m) by 

& (°° P()-V(-) 

Bn Jo 

This is presumably negative when V’(0) > 0, since the numerator of the 
integrand is then negative on the average. ‘Thus the correction (136) due 
to the cut-off is a reduction, as expected. 





dl. (136) 
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Actually, it was shown in Lighthill (1957 a) that as well as the doublet 
of strength (134) the upstream influence of the sphere includes also a rather 
smaller contribution due to the ‘secondary trailing vorticity’. ‘his is 
equivalent to the effect of a line of doublets with axis the y-axis stretched 
all along the positive x-axis. However, it may be shown that the displacement 
resulting from this line of doublets is estimated correctly even if the cut-off 
studied above is neglected, and hence the validity of the correction (136) 
remains unaffected. For the practical significance of this correction, see 
Lighthill (1957 b). 
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to its length 
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SUMMARY 

A flat strip, of infinitesimal thickness and infinite length, 
is located in a viscous incompressible fluid, and both the strip 
and fluid are at rest initially. ‘The strip is abruptly set into 
steady motion parallel to its length. An unsteady uni-directional 
flow of the fluid results, and there is a variable skin friction on 
the strip which must be overcome to maintain its velocity. In 
the early stages of motion, the skin friction is large, with a local 
behaviour which resembles that for a flat strip of infinite width. 
The skin friction near each edge of the strip can be more accurately 
represented by referring to the semi-infinite configuration that is 
realized on displacement of the other edge to infinity. At this 
stage, the results depend on separate consideration of the two 
limiting configurations, and the way to further improvements is 
not clearly delineated. ‘The object of this paper is to provide a 
formulation which contains the preceding information and allows 
a systematic evaluation of all additional refinements. ‘Thus, it is 
shown that the total skin friction on a strip of width 2a, moving 
with velocity V in a fluid whose coefficients of viscosity and 
kinematic viscosity are yu, v, takes the form 


D . a 2a ; 2 [ 1 i xf au ; 
aN rem 5 |, eal ee eee 


Here the first two terms stem from the infinite and semi-infinite 
strip distributions, and the integral, containing a complementary 
error function, furnishes all corrections of the lowest exponential 


order, in the initial stages of motion, when a/,/(vt) > 1. The 
calculation is based on a new integral equation which makes 
explicit the effect of an isolated edge, and by iteration provides 
the interaction effects between edges at a finite separation. 


1. INTRODUCTION 
A simple class of viscous incompressible fluid motions is that produced 
by the forced motion of a solid parallel to its infinite length. If the solid 
is instantaneously set in motion with constant velocity, and has a uniform 
aspect along its length, the resulting fluid motion is uni-directional, 
although unsteady in time. The fluid velocity is parallel to and uniform 
in the direction of motion of the solid, and hence described by a linear 
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partial differential equation with the time and transverse coordinates as 
independent variables*. Rayleigh (1911) analysed the simplest con- 
figuration, a flat strip of infinite width, for which the skin friction is 
uniform and has the magnitude pV /4/(zvt) per unit area (of each face). 
The strip of semi-infinite width and the more general wedge shape have 
been studied by Howarth (1950), Sowerby (1951), and Hasimoto (1951 a, b), 
and from these calculations the edge effect can be isolated. A characteristic 
dimension is introduced by the strip of finite width, say 2a, and the state 
of affairs depends on whether a/,/(vt) is large or small compared with 
unity. In the former case, the fluid motion is confined to a boundary 
layer with small thickness relative to the strip width, and the total skin 
friction is given approximately by supplementing the Rayleigh contribution 
with those of the isolated edges, whence 


4 
D= mal a +2| (a/\/(vt) > 1). (1) 


Vv (at) 

This result, obtained by Howarth (1950) and Hasimoto (1951 a, b) (Davies 
(1950) also considers the flat strip, presenting a formal series solution 
which does not lend itself to analytic approximation for small ¢), omits any 
interaction between the strip edges, and it is such effects that we propose 
to calculate. It turns out that their magnitude is rather small, in the domain 
where the skin friction may be analysed by the preceding scheme, for the 
corrections to (1) are of exponential order; thus 


D=-pgF Re sane ee (te) eit + Og (see a*/vt | 
\/ (zvt) Ln a 4 a J 
(a/4/(vt) > 1). (2) 
‘The developments of the skin friction for the strip and the (smooth) circular 
cylinder show a quite different form after the first two terms, since the 
latter is, according to Batchelor (1952), 
27a Trvt mvt 

Du pr Rar +a7— = + ——+ + 7 (a/\ (vt) > 1), (3) 
where a denotes the cylinder radius. For the later times, 1/(vt)/a > 1, 
Hasimoto (1954) and Batchelor have obtained other developments of the 
skin friction in the two cases; there is a smooth transition to the approximate 
expression (1) at a/4/(vt) = 1. Hasimoto (1954, 1955) has recently 
generalized (3) so as to apply to a smooth general cylinder, but the 
corresponding accuracy is not yet attained for a polygonal cross-section. 
The calculation to be described makes use of a double transform 
operation, one being a Laplace transform with respect to the time and 
the other a complex Fourier transform with respect to the (transverse) 
coordinate in the plane of the strip. After applying the Laplace transform 
to the differential equation for the fluid velocity, the resulting boundary 


— 


* An equivalent heat conduction problem relates to the variable temperature 
distribution in a uniform medium which is initially at zero temperature, and in 
which an internal boundary is maintained at a constant temperature. Here the 
heat flux at the boundary is analogous to the skin friction in the fluid motion problem. 
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value problem is formulated (§ 2) in terms of Green’s functions and integral 
equations. ‘lhe subsequent analysis is effectively carried out with complex 
Fourier transform and elementary function theoretic arguments, yielding 
finally a new integral equation (§ 3) which is well adapted to the determination 
of the skin friction (§ 4) in the initial phases of the motion. It is noteworthy 
that the inverse Laplace transform back to the time coordinate can be 
performed without approximation at any stage of iteration for the integral 
equation. 

Similar Fourier transform techniques can be employed to analyse the 
scattering of sound or light waves incident on a strip, and the results invite 
comparison with those in the problem at hand. In the time-harmonic 
two-dimensional situation, with primary plane waves, the scattering cross- 
section at normal incidence relates to an integrated strip distribution and 
is thus analogous to the total skin friction. For comparison purposes, 
consider the skin friction (3) and the leading terms in the cross-section o 
at wavelengths A short compared with the strip width 2a, 


4a/ xX \32 4na ot 2af2r\2 Sra ; 
Co = Aa = 7a (aes) cos “5 — 3) 5 a (5-5) acai (A a<s 1), (4) 





when the total wave function vanishes at the strip. The counterpart of 
the Rayleigh term in the skin friction is the geometrical term of the cross- 
section, and whereas the next contribution to the skin friction arises from 
isolated edge effects, these do not contribute to the cross-section. Succeeding 
terms of the developments, which account for interaction between the strip 
edges, are of exponential form in the skin friction and of triganometric 
form in the cross-section, reflecting the diffusion and propagation character 
of the problems. 

Carrier (1955) has also studied the boundary layer aspect of solutions to 
certain types of integral equation, where reference is made to a semi-infinite 
configuration for the purpose of characterizing the edge behaviour. ‘The 
solution to an inhomogeneous finite-range integral equation (cf. (21)) is 
given, neglecting interaction between the edges, and utilized in analysis 
of an oscillating plate in a viscous fluid. 


2. FORMULATION 
Let the strip be located in the (x, z)-plane, 
—-a<x<ta, —-o<8 + 00, 
and set into motion at time ¢ = 0 with steady velocity V along the z-direction. 
The viscous fluid begins a parallel motion, with velocity v.(x,y,t) which 


satisfies the equation 
- = = 
ov 02v, ou, a 
— =>-wras>+ aS b, 
oO (oe at) ( ) 


(v is the kinematic viscosity) and automatically complies with the equation 
of continuity, by virtue of translational invariance. ‘The same differential 
equation applies to the velocity ratio 
u(x, y, t) = v,/V, (6) 
K 2 
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which is further subject to the initial and boundary conditions 
p= ati=O, (7) 
%=1 aty=0 and |x| <a, | 
ob ~O as x? +? > o, J 
and the symmetry condition 
u(x, y, t) = B(x, —y, t). (9) 
If 
U(x, Vy, p) =p | ePhb(x, y, t) dt (10) 
0 


denotes a quantity proportional to the Laplace transform of ¢/(x, y, ¢), then, 
using (5), (6), (7), it follows that 


x” 07 
and 
%=1 aty=0 and |x| <a, | 
~ r (12) 
p>0 asx*+y?> 0, | 
where 
k? = p/v. (13) 


The boundary value problem stated in (11), (12), is advantageously 
reformulated with the help of solutions to the associated inhomogeneous 
differential equation 


{; nf a — )G (ay x,y’) = —8(x—x')d(y—y’) (Atk} > 0). (14) 


ox" oy 








A fundamental solution of this equation with logarithmic singularity at 
the ‘source’ point x = x’, y = y’, and a null value at infinity is the Green’s 
function 
VP en eee es (a. ee ee: 5 
Go(x, V3 x,y’) = (27) 1K (Rv/ {(x— x’)? + (vy —y’)*}), (15) 
where Ky denotes the zero order Hankel function of imaginary argument, 
namely 
, - *77(1)/ ¢,in/2 
K,y(¢) = $atH}(Ce"*). (16) 
In terms of the source function G5, a solution of the differential equation (11) 
which vanishes at infinity and is symmetric with respect to the plane of 
the strip can be represented by 


= Ley x. ae * ; ; 
Wx, y,P) = 5- | KolAV/{(w—a')? +y*Vg(x') de’ (17) 


a < a 
Placing y = 0, and writing 
b(x, 0, p) = f(x), (18) 


we have 


1 se 
| K(k\x — x’ |)g(x’) dx’, (19) 





f(x) 
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and from the condition (12), which gives f(x) = 1 for |x| < a, there follows 
an integral equation for the (vorticity distribution function) 


g(x) = -2(30) (20) 


on the strip. ‘The boundary value problem (11), (12), is thus linked with 





the solution of a one-dimensional integral equation, 


1 ea 


Lx K (Rix —x'|)g(x’) dx’ (|x| < a), (21) 


. a 
and a knowledge of g(x) determines % at all points off the strip by (17). 
Choosing next the Green’s function of the half space y > 0, 
Gi(x,93 «,9") = (2m) [Ky (v— x2 + (yy) - 
— K (kv {(x-— x’? +(y+y’P})], (22) 


which vanishes at y = 0, a representation for ¢ therein is, 


wy l a r Cf "\9 ) 
Wx, yp) = — oa | K(k {(x— x’)? +y?}) dx’ — 
. 7 OY / |z’\|<a 5 
I a r / ‘\2 9 4 , , 2) 
oe | K (Rv/{(x— x’? +y}) f(x’) dx’ (y>0), (23) 
where the boundary value (x,0,p) = f(x) =1 for |x| <a is utilized. 
A further differentiation in (23) with respect to y, and passage to the limit 


y = 0, yields 

















+ | K (Rv {((x— x’)? +" b)f(x’) ax’) 


r 


1 / ad . j . , i 
- -(i -k) (| = K,(R\x—x’|) dx’ + 


+ K, (Rix —x'|) f(x’) dx’), (24) 


| J |x! a 





for the K,-function is a solution of the homogeneous equation (11) when 
y’ =Oandy +0. Since f(x, y, p)=(x, —y, p), it follows that (eb/dy),, — 9 = 0 
for |x a, and the condition g(x) = 0 for |x| > a provides an integro- 
differential equation, 


ad? “ i 2 P , 
par 


rm | . K,(k\x —x'|) f(x’) ax’) = 0 (|x| >a), (25) 


for the transform ¢(x,0,p) = f(x), in the domain |x| > a. The solution 
to the boundary value problem (11), (12), can be based on this integro- 
differential equation, and the representation (23) for the even function of y, 
u(x, y, p), in the half-plane y > 0. It is evident that the two formulations 
given are of a complementary nature, with the primary role accorded to 
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the distribution g(x) on the strip in the former, and to the coplanar 
distribution f(«) off the strip in the latter; clearly, a knowledge of either 
function suffices to determine the other. 

The skin friction D(x) at the strip is 


Dx) = 4 (52) = VEL) (26) 


where Z~! denotes an inverse Laplace transform, namely (recall the 
additional factor p introduced in (10)), 


E-4 p(x, p)| = = | ‘i ee dp (c>0). (27) 


l 


“¢-i« 

In the total skin friction or drag, D, the integral of D(x) over the strip width 
is involved, and an additional factor of 2 is necessary to account for both 
faces, whence 


p= pk (| 


Although the drag is naturally related to g(x), the latter function need not 
be the primary quantity in an integral equation formulation of the boundary 
value problem; to calculate D in the early stages of motion, it will in fact 
be convenient to work from the integro-differential equation for f(v) and 


a 


2(x, p) as = val ; ~ oem [" 2(x, p) dx. (28) 


—-1x 


obtain g(x) secondarily. 


3. INTEGRAL EQUATIONS 

The basic integral equations (21) and (25) cannot be solved explicitly 
in closed form, and, as presently constituted, do not lend themselves to 
approximate solution, say, by iteration. Our objective is to obtain a new 
integral equation from (25) which is appropriate to the conditions prevailing 
at small values of t, and which also permits ready approximation to the skin 
friction. ‘The idea underlying this transformation is that, for large values 
of the parameter a/,/(vt), the distributions f(x) in the ranges . > a and 
x a are weakly coupled and exactly describable in the limit of no 
coupling, which corresponds to the geometrical configuration of a semi- 
infinite strip. 

A brief exploration of the latter problem is in order, and information 
obtained relative to the behaviour of f and g near the edge will carry over 
to the finite strip case. ‘The relation (19) provides a convenient starting 
point, and, if the origin is shifted to one edge, this relation becomes 


f(x) | Ko(R\x—x’|)g(x’) dx’ 


2a So (29) 


1 forO <x 2a. 


Whence, letting a —> ©, we obtain 
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as the form appropriate to the semi-infinite strip 0 < x < «©. Employing 
the integral representation 


1 ‘ x ie . x" : 24 ‘2 — y’| , 
Hersey». SS 
4 : \ (k? t 2) 
(31) 
and introducing 
f(2)= | e®f(x) dx, (32) 
o(f) = | e®" g(x) dx, (33) 
we then obtain the transform version of (30): 
- l &(S) , 
O+ i= Vere si 


Here « appears by virtue of an exponential attenuation factor introduced 
when calculating the transform of f = 1 in the range x > 0. The relation 
(34) has meaning if there is a common domain of regularity in the complex 
¢-plane for the functions occurring therein. According to (30), f() contains 
the factor e*” as x-> — «, and hence f(€) is regular in the half-plane 
JI \C\ > —Atk\; an exponential free behaviour for the strip distribution 
g(x), as x > ©, implies that g(¢) is regular in the half-plane /{Z} < 0. 
Thus, the common domain of regularity is a strip parallel to the real axis 
of the ¢-plane, defined by 0 > ./{Z} ~ Aik = J\—-ik\. Following 
the Wiener—Hopf procedure, (34) is rewritten in the form 

asin, ViC+iN- Vik +i) VR) HD) 

V(o+ 1k) f(¢) + a ee + 7G 


és e+1C e+1C 








j 
( 
and the right-hand side is regular in the (lower) half-plane 7 
integral function is thus defined throughout the ¢-plane, and, since each 
side vanishes at infinity in the respective half-planes (assuming only that 
f(x) and g(x) are integrable at the edge x = 0), this function is zero. Hence, 


ik e 
i= (eta) : 


A(t) = 2 (iky VE) 


J) = 


in the limit « = 0, while 





ux1 


(37) 


e+7C 





; > oo a ij? 
Sap Fe ak - aaa eee a a : ail 38 
2a a I J ~ :) au | 7 | 1 ur/(u- a 0h 


Aa — /(6—-ik) de = -— =| Vu +1) du 


€+1¢ T 








=2k+— | —— vV(u-1)du (x>0), (39) 
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with k real. From (38), it appears that 
1a du 
lim f(x - ~ <= 0) 
shay 9) 
which is the value of f for x > 0, and 


f(x) ~ e&*/(—kx)'? asx> — ow. (41) 
According to (39), 
g(x) ~ (kx)? asx 04 (42) 
and 
p(x) ~ 2k asx>o. (43) 
Since 
l 4 K,(kix—x’'|)(2k) dx’ (-« Pa eo) (44) 


the uniform distribution g(v) = 2k evidently refers to a strip of infinite 
width. ‘This synopsis of features of the semi-infinite strip configuration is 
sufficient for the present, and further references are deferred until later. 

Turning attention again to the strip of finite width, and choosing (24) 
as a starting-point for the analysis, we find that the concomitant transform 
relation becomes 


1a(0) = V(k2 +2) | FO ; (— 7) | (45) 





where 2(C) | e g(x) dx = g(-9), (46) 
and f() =i i(¢) + fol), (47) 
with fA) = | o*“f(x) dx, (48) 
7 pe , 
FAS) = | e f(x) dx. (49) 
Evidently g(¢) is regular in the finite part of the ¢-plane, and with the local 
behaviour g(x) ~ (k(a¥x))-!? as x +a, which is characteristic of an 
isolated edge (see (42)), it follows that 
a(L) ~ Fal FIL? as > o (MGSO). (50) 
The exponential decrease of f(x) as x — © implies that f,(¢) is regular in 
the (lower) half-plane .7 {Z' A\k', and a finite limit for f(x) as wx > a+0 
(see (40)) yields the asymptotic behaviour 
Fi) ~e®lE as i+  (F{L} < AIk}) (51) 


in the half- plane of regularity. By symmetry, f is an even function of x, 


so that - o£) 4A —), and thus f(0) is regular in the (upper) half-plane 


SJ {C3 —&‘k*‘, where 


fall) ~ e®4/f as |fj >a  (S{L} > —AIR}). (52) 


Accordingly, the transform relation (45) holds in the strip |-7{¢}) < A{R} 
and, after the substitutions 


f(0) = e*F Yd, 


f(0) = e"F (2), 








—— 
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takes the ni 


1 o(L) == 4/(Rk? 4 *) Me fa F (L) + e%4 F(0) + (e®¢ —e®4)/iZ]. (54) 
Mult: Si: in turn \ e'°", e- =", we get 
bee = P+ QF (+ Fd) + (2 1d), (55) 
ba(Qje ee = 4/(R2+ 2)fe 2 (0) + F,(0) + (1 —e*4) ig]. (56) 
On the factorization of the cai 


. (R? t (4) =“ 1(¢ t tk\(C -ik)', 


coakiat «ae ste (a ee ee ene] 
22(C S) (E+ ik) =v/(¢ -it) | FO: ee Fl) 4 —F I; 


The left-hand member is regular in the upper half-plane, and the first term 
of the right-hand member is regular in the lower half-plane. For each 
remaining term of the right-hand member, a decomposition may be effected, 


with the help of Cauchy’s integral theorem, 
1 f w(€) 
gol = ERR ied Cc 
W(z) = ia? é 2 dé, 


into a sum of functions regular in the upper and lower half-planes, 
respectively. Collecting all terms regular in the respective overlapping 
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Figure 1. 


half-planes, an integral function is defined throughout the ¢-plane, with 


the following representation in the lower half-plane: 











; edie 1 | ee 
I.F. = 4/(€—1ik)F,(¢) + a | 4/(€—tk)e"' FE) 2 + 
] (ea fe"5"_|\ d€é ne 
S Oni | v(E=ik)( = | a de (57) 
Here C_ is a path of integration conducted in the strip -4{C} AKik* (see 


figure 1), and ¢ is located below the path. ‘The terms in (57) tend to zero 


x% in the lower half-plane, as do the corresponding terms at 


when ¢ 
and thus the integral function vanishes, 


infinity in the upper half-plane, 
— 
“ 7 ee 
V/(f — tk) F i() + a |. % (€ —th)e* FX(E) z t 


aie 
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Similar considerations applied to (56) give 


cial ae ae ae eee ee 
V(b+tR) FAD) + a | VUE +tR)E Pm PE) 


L ff Pe oe ad ‘ 
ral, Vern (——) = 0, 69) 


o-6 








where C_, is another path conducted in the strip, this time below the point . 

The integrals on C_, C., are usefully recast by deforming the contours 
in the upper and lower half-planes, respectively, so as to run along the 
branch cuts drawn from ¢ = +7k. If these cuts are placed on the imaginary 
axis away from the origin (.7{k} > 0), then, taking account of the phases 
of the radicals at each side of a cut, it follows that 




















ela re e-tkau _ | du 
: - p12 ee 
vie Fo) “* " 1)( u \as —C = 
ei! . ~ du 
___.. Bal = “(ab p—2kau 
_ k |, VV (u— 1) F(iku)e a? 
- ens 42 ( ea/(u—1) du A 
om dy us tku—¢ 
eizi4 "ik | du 
= L3 2 =i F. of —2khau 7 
1 | _ ali (F aa iu) thu —€’ ey 
and 
_— eix/4 “© 4/(u—1) du 
ap a = ee [1 2 : + 
V(b +tk) F (6) 7 ; 4 u—s- thu + 
eriz/4 "a ~ 1 du 
+ —— k3? —1){ F\(—iku) — — )e?*«——_.. 
T wn ( 1 ia) thu+¢° 


(61) 
The integral equations (60), (61) are not independent, for reversal of the 
sign of ¢ in (60), and the use of relation F\(¢) = F.(—2), leads to (61). 
Let ¢ = tkv in (61); then 
1 ¢*/(u-—1) du 
ak), u u+ov " 











V(v+ 1) Fike) = 


leit du 
: el = ___._. |p -2kau 
. i V(u-1 )(F (tru) j ). ae 
; il V(vt]) |. 
ki kv , 
Lif = l di 
een / a Phe ee 2kau_ 9 
ce q \/(u 1) ( F,(ieu) iu)’ ee (62) 


which is an integral equation for the function F,(ékv). When a> oo, 
the integral drops out, and the resulting form of F, agrees with that found 
in the semi-infinite strip configuration, namely (36). A more symmetrical 
version of the integral equation is obtained by defining 

P(v) = : 
U 


—k F,(ikv), (63) 
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for then 
: tf du 
V(e+1)P(v) == + — | V/(u—1)P(u)e?kan eae (64) 
From (62) it also follows that 
= . 1 re , o1.,,,, dU , 
F,(0) = F,(0) = ae F a/(u—1)F (u)e*hau — , (65) 


The relations (64), (65) provide a suitable basis for calculating the skin 
friction in the early stages of the motion, as will now be verified. 


4, SKIN FRICTION 
To begin the skin friction calculation, an expression for g(0) is required. 

Using (54), (65), this is given by 

2 _ 4 “i du 

g(0) = 4ha+4k F,(0) = 4ha+2-—- | +/(u—1)P(uje?*™ yon (66) 

“i 1 

where P(u) satisfies the integral equation (64). An iterative procedure can 
be applied to the integral equation when #{k} > 0, and the leading terms 
of P(u) thus obtained are 


Plu) = 1 l ‘fe v—1\e-**er du - 
“) ~ uy/(ut+) * 74 wrt |, 2S vo utv (67) 


The corresponding approximation for g(0, p) is 


' ‘ 4 4 / ag du 
o(0, p) = gM, ed oes 
— Bh a 2 ae.” (; i i) u2 
4 Ni L / (u = 1)(v- 1) 2av(p|Y)(u+v) dudv 
” a | iN (eaiNe + 1) j os (68) 


after reverting to the Laplace transform variable, 

















p= vk’. (69) 
There remains only an evaluation of the inversion integral (28), 
in Po dp . 
D=75) e"aO0)> (¢ > 9), (70) 


and it will be noted that the conditions #{p} > 0, Aik} > 0 are indeed 
compatible with (69). ‘The resulting skin friction is 


2a nO u—1 au \ du 
| a er mere a 
D= Zul | - (zvt) | ; JG a i)erfe(- =) u 


2 x / (u—1)(v—1) wr a(u+v) dudv = 
“lA (Gainer) \ a) uo(u = “7 


Alb 








where 





~we~lvNnz for z > 1 


is the complementary error function. 
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When a/,/(vt) > 1, the principal contributions to the integrals arise 
from the neighbourhood of the lower limit, and thus the double integral 
is smaller than the single integral by a factor e*“/vt. This exponential 
decrease in magnitude evidently characterizes further contributions to the 
skin friction resulting from continued iteration of the integral equation. 
Successive terms acquire comparable magnitude when a/,/(vt) < 1, and, 
in the limit a/,/(vt) + 0, the full iterated series is required to compensate 
the time independent term in D. The development (71) is therefore useful 
only for the early stages of the motion, as anticipated. 

The first term in (71) arises from a uniform, local skin friction character- 
istic of an infinite strip, and the next term takes account of the non-uniformity 
introduced by the strip edges, here computed additively when the latter are 
infinitely remote. ‘lo verify this interpretation, observe that the excess 


skin friction on the semi-infinite strip 0 < x < 0, attributable to the edge, 
stems from the distribution (see (39)) 
zn (? du 
g(x) =— | V(u-le*"— (x > 0), (73) 
7), 
whose integrated value over the range 0 < x < © is unity. The second 


term in (68), and hence (71), thus bespeaks two such contributions, from 
remote edges. If the distribution (73) is integrated over a strip 
0 < x < 2a, there is an additional contribution which is evidently distinct 
from the third term of (68). ‘The necessary refinement of g(x) entails a 
direct consideration of the finite strip configuration, and it is of interest 
to see how this may be carried out in the transform formulation. 


In |x| < a, by (24), 


bale as See as [” itp 2 Sin a dt 
~ 28") = 77 ae —* ) | a. C 4/(R+2) 





“ox 


| eb efe-Ba F(0) +e FAOSarcE | ia 





R\ | ] - 
£8 I(t ny V(u—1)P(uje* how F,(—{). (75) 


tku—¢ 


A modification of the #) | is expedient, and results from deforming 
the contour into the upper or lower part of the ¢-plane. Thus, for |x| < a, 








ro 2 sin Ca dv 
8° ft VR 
= 3 | . (1 — e+ 9) — (1 — e+) + 
(fee SN pe i By 
- : | 2g re ge orp (76) 
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and 


ae . er dé 
| ent (a- FE (C) —s $ 


2 -% e-kla—a)u l du 
“5, = 0- ye) ye 


) > 00 k(a—ax)yu L 9) 
2 € dv ) du 
J a l P a 2kav 
— ———_—_—___—____——- < q — OU é on —— > Se a ae / / 
=a Vert) J, vi ee) u+v} /(u?—1)’ 7) 
while the corresponding integral with /*,(¢) calls only for a reversal in the 
sign of x in the above expressions. Substituting these expressions in (74), 














we obtain 


a d2 \ fp 2 @e k(a+axc)u 4 e k(a—.r)u 
g(x) = 2k —|— -kR : 
g(x) c+ ole (73 ) | | :' \ (w ; 1) x 


ye dv du 
: Lae 9 —1)P(o)e-2kav { 
" | vie ee ut+v | Toh | 








(1 -{* dv 
ae abe coe UY — D v)e akav | 
lu om J, vie ee u+v| au 
2k ;@ 
— 2k-4 na J Vv (u* — 1)(e k(a+ar)u —€ k(a—a “) P(u) du, (78) 


where the last form follows from a use of the integral equation for P(u). 
In the second form, there may be recognized contributions of the type (73), 
referred toeach stripedge. ‘lhe singularities of g(x) atx = + aare contained 
in factors (a+ x)~!*, which are also characteristic of an isolated edge. ‘The 
time varying function g(x, ¢) obtains after the inverse Laplace transform 
is applied to (78), and, with the first approximation P(u) = 1/uy/(u+1), 


the result is 
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After integration of (78), one finds 
oa 4 po 
| &(x) dx = 4ka+ - | s/(u—1)(1 —e- 2") x 
— i al i | 
: i dv \ du 
x{-4- ‘(vw — 1) P(v)e* #2” ——_>— (802 
Lu ro. lla lk u+of u iad 
4 (* = Rica a du ; 
= 4ka+ - | _Vv (u? — 1)(1 —e-***) P(u) —) (80 b) 


in contrast with (66). However, on substitution of the iterative expression 
for P(u), the same development is obtained with each of these representations ; 
it is noteworthy that (66) and (80 a) lead to the result (68) when P(w) is given 
by the two-term approximation (67), whereas a three-term approximation 
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is called for in (80b). This feature reflects a difference in basis of the 
calculation according as the off-strip distribution f(x) for \~| > a, or the 
strip distribution g(x) for |x| < a, is involved. ‘The semi-infinite distribu- 
tion is a natural first approximation for f(«), and each successive stage of 
approximation generates a distinct order-of-magnitude correction in both 
f(x) and the skin friction D. On the other hand, only a finite part of the 
semi-infinite distribution for g(x) is relevant, and successive approximations 
are partially coupled in magnitude, resulting in like order contributions to D. 


This work was supported by the office of Naval Research. 
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SUMMARY 


The motion of a viscous fluid contained between two rotating, 
circular cylinders whose axes are set slightly apart is considered. 
The equations of viscous motion are linearized by expanding the 


stream function in the form > %&, y”, where yisa parameter which 


depends on the distance between the cylinder axes. ‘The ensuing 
analysis appears to hold for all values of the fluid viscosity v, and 
in particular for small values of v. 

The asymptotic behaviour of the solutions for small v is 
examined, attention being mainly confined to the first order 
stream function #, and the corresponding component of vorticity 
¢,- Outside the boundary layers, where, for small v, 4, may be 


I 


expanded asymptotically as > ¢'}’v"?, the terms ¢") of the corres- 
ponding expansions for the vorticity are shown to be uniform 
throughout the fluid. It is noted that the asymptotic expansions 
of ¢s, for the region of the boundary layers and for the region outside 
the boundary layers may be combined in a single expansion which 
holds in both regions. ‘The leading terms of this expansion are 
calculated by boundary layer methods. 


1. INTRODUCTION 

This paper is concerned with a viscous motion whose governing equations 
may, in principle, be solved exactly. ‘The asymptotic behaviour of the 
solution when the fluid viscosity v is small is derived and discussed and a 
procedure formulated by which the higher order corrections to the boundary 
layer approximation may be determined. Inasmuch as the streamlines 
of the motion are closed the results bear particularly on closed flows. 

Now, the boundary layer approximation, as conceived by Prandtl, 
determines the stream function ¢ of a steady two-dimensional motion with 
an error which is, in general, O(UL{v/UL}) in the boundary layer and 
O(UL‘v/UL}"?) elsewhere, U and L being respectively a typical speed and 
a typical length of the motion. For certain configurations even this 
accuracy cannot be attained. If for example the fluid impinges on a sharp 
edge as in the uniform streaming pasta flat plate, the Prandtl approximation 





160 W. W. Wood 


fails to determine the stream function correctly near the edge. More 
serious difficulties attend the shedding of a boundary layer to form a wake. 
Such configurations aside, the problem remains of how best to determine 
the higher order corrections to the original approximation. 

One correction can be made at once; that is, the displacement of the 
streamlines outside the boundary layer caused by the slow movement of 
fluid inside the layer can be calculated. Ina recent paper, Kaplun (1954) 
has shown that by suitable choice of coordinates the solution of the boundary 
layer approximation can be made to hold throughout the whole region of 
flow and made to include, outside the boundary layer, the effects of dis- 
placement thickness. 

In order to formulate a procedure for determining the higher order 
corrections, it is helpful to have as a guide one or more exact solutions to 
the equations of motion, whose asymptotic behaviour for small v has been 
elucidated. One such solution, though for an unsteady motion, has been 
analysed by Lagerstrom & Cole (1955). Lagerstrom & Cole consider 
the motion set up by a uniformly expanding circular cylinder moving 
parallel to its axis in unbounded fluid. ‘The exact solutions obtained by 
Lagerstrom & Cole are of the form 

pr > ZB, (r, tv)" “eG *). 
where the coefficients -#, are significant in the boundary layer and 
elsewhere are transcendentally small, and they propose, from general 
considerations (which are not stated), that the stream function of steady 
motion with no wake may be generally represented in the form 
bpr~ > SF, (r)yv"? + > BZ (e, vv"? +0(v-™), 


i ) nm 0 
where the first term represents the truncated asymptotic series for % at a 
fixed point r. 

‘The example considered in the following concerns the motion of fluid 
between two rotating, circular cylinders whose axes are parallel and set 
slightly apart. When the separation of the cylinder axes is small compared 
with the difference in length of their radii the equations of motion may be 
solved by a perturbation method. Moreover the perturbation expansions 
appear to hold uniformly for all v, providing that the expansion parameter is 
sufficiently small. ‘The solution therefore determines the asymptotic 
behaviour of the motion when v is small. 


2. LINEARIZATION OF THE EQUATIONS OF MOTION 
We begin by defining the perturbation expansions and deriving equations 
of motion for the first order terms. 
Let the radii of the two cylinders be a, b and the distance between their 
axes ae. ‘lake polar coordinates r, @ in any cross-section, with origin at the 


axis of the inner cylinder and reference line @ = 0 along the radius to the axis 
of the outer cylinder (see figure 1). _ Introduce further the non-dimensional 
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coordinates p, ¢ defined in terms of 7, 6 by 





Z=6¢ 7 —, 
ry | (1) 
' = a", C = pe™, ) 
with 
y = —2e[(b/a)? —1—«* + +/{(62/a? — 1 —&?)? —4e*}} 1. (2) 
The circular sections of the inner and outer cylinders are then the coordinate 
lines p = 1 and p = f, where 
; (b/a)+e-y ss 
en Raat (3) 
1 —(b/a)y —ey 
} 
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Figure 1. 
Che remaining coordinate lines p = constant and ¢ = constant constitute 
the two orthogonal pencils of circles generated by the circular cylinder sections 
} The coordinate mesh is clearly the same as that of the two-dimensional 
bipolar system in which the cylinder sections are coordinate lines. ‘The 
coordinate variables p, ¢ are however different, having been modified so as to 

\ have the advantage in this problem of degenerating to polar variables r, 6 
when the cylinders are coaxial (« = y = 0). 

The rotary motion between cylinders is supposed two-dimensional and 

, may therefore be represented in terms of a non-dimensional stream function 
, the velocity components in the coordinate directions p increasing, ? 
increasing, being respectively 

) vJ Ou Oys 

ul, = 9,— => ul, = —q,\J =, (4) 
Z Pp Op op 
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where q, is the peripheral speed of the inner cylinder. ‘The function J is 
in effect the Jacobian of the transformation (1) and is defined by 
(1 +2ypcos¢+ yp)? 
(d-pPp 
On eliminating the pressure, the equations of motion reduce to 
1 (Ys, 6) = 1 2F C= —JV%a, (6) 
p Op, d) R 
where ¢ is the axial component of vorticity, R is the Reynolds number, 
defined here by agq,/v, and, in both equations, 
o aoe 2 ™ 
op” pap * pap 
(In using q, as a typical speed, it has been tacitly assumed that the inner 
cylinder is not at rest. If it were at rest, the typical speed could with only 
trivial changes be taken to be the peripheral speed q,. of the outer cylinder.) 
The boundary conditions at the inner and outer cylinders are respectively, 





5 


JT 
— 


= 


0 ] 
(1, = 0), — (1, eT > a. eS 
(9) ap I #) Vv id(L, g); 
(3) 
(8, b) = constant, ap WB, 6) = — av (Bo): . 


Further the solution to (6) must be such that the pressure is single-valued. 
These equations and boundary conditions suffice to determine %. From 
them it is clear that 4 depends on the geometrical parameters f, y, the ratio 
of the peripheral speeds g,/g., and the Reynolds number R. 

It is now proposed to consider the motion when the separation ae of the 
cylinder axes is small compared with the difference in length b—a of their 
radii (that is, y<1). In the particular instance where the cylinders are 
coaxial the motion is axi-symmetric and independent of the Reynolds number. 
The stream function is then 





bo(p, B, 92/41) = - LA p" = B log p, (9) 
where 
(qo —] 8? — B(qs/q,) 
A = Pldaidi) — If 6 P— Pde! 44) 4) : (10) 
p?-1 p?-1 
and the vorticity is 
Co(B, G2/q1) = 2A. (11) 


In general, when the cylinders are nearly coaxial, it is assumed that the stream 
function % and the vorticity component ¢ may be expanded in the forms 


p= bo(p, f, 42/41) “4 2 nlp, d, B, 92/415 R)y", | 


n 


(12) 


va 


(1 =e = 2A(f, 2/91) + 2 Ci (ps ds B, G2/ 91> R)y". 


n=1 
It is the first order motion described by the first order terms with which this 
paper is principally concerned. ‘The higher order terms are discussed briefly 
in the final section. 
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Concerning the mode of expansion, note that the coordinate description 
of the bounding surfaces does not involve the expansion parameter y. ‘This 
is important if the expansion is to be used for large Reynolds numbers, and is 
the reason for introducing ‘bipolar’ coordinates. If the coordinate des- 
cription of the bounding surfaces involved y then the boundary conditions 
would have to be applied at surfaces y = 0 different from the cylinders. For 
this purpose, the velocity field at each boundary would have to be expanded 
in powers of ‘distance’ from the surface y =0 to the boundary, and it is 
unlikely that such series would converge when R is large. Note also that 
when y 0, the motion described by %(p, g2/g,, 8) is not axi-symmetric. 
In this motion the circles p = constant and, in particular, the cylinders are 
streamlines, the speeds at the inner and outer cylinders being respectively 
avd, y)} q2V J (B, d, Y)5- 

When the expansions of #, ¢ and the analogous expansion of J are sub- 
stituted into the equations of motion (6) and the boundary conditions (8), and 
the leading terms isolated, we get 


B\ 0, | , 
(4 5 =) ob = R V Gp (13 a) 
C, = — V*(b,— Ap’ cos 4), (13 b) 
with the boundary conditions on ¢,: 
da(1,4) = 0, 5, Yall, $) = 2cos4, 
(14) 


ws,(8, 6) = constant, a (8, 6) = 2 qo/q, Bos d. 
op J 


Further, from the corresponding exact condition, the solution of (13) must be 
such that the first order perturbation in pressure is single-valued. 

The equations governing the first order motion contain the Reynolds 
number in the same way as the exact equations but have the merit of being 
linear and hence more tractable. In that the convective velocity in the 
vorticity equation (13) is known, this equation is mathematically similar to 
Oseen’s equation and the generalisation oi it discussed by Zeilon (1927) and 
Burgers (1921). Inthe present case, there is, however, reason to believe that 
the perturbation to % described by the equations (13) is uniformly valid to 
O(y) even at large Reynolds numbers. A similar linearization of the 
equations of motion was discussed by Proudman (1956) who investigated the 
flow between two concentric spheres which rotate with slightly different 
angular velocities. In the simpler configuration considered here the 
linearized equations may be solved exactly. ‘This we proceed to do. 

3. "THE FIRST ORDER SOLUTION 

It is evident from (13) and (14), that %, and ¢, vary with ¢ according to the 

simple forms 
1 = Ap? cosd-+2A(f (pet, 


22 g(p)e*}. 


II 


$y 











164 W. W. Wood 


The way in which %, and ¢, vary with p is then determined from the differ- 
ential equations 





a RBi-+ =f 1 
all ie gay, 
' 5. ie | (16) 
if f iis J —- —o 
oe 
with the boundary conditions 
f(1) = —3A, f'Q) = 1-34, ] 
| 
. 17 
fis) ee 14 B, f'(B) = . —3A B, | ( ) 


where primes denote differentiation with respect to p. 

Provided that the vorticity of the motion described by ys, is not zero 
(A +0), g(p) satisfies a Bessel’s equation of order wp = 1/(BRi+1) with 
independent variable z = \/(—ARi)p. ‘Thus 


2(p) = —2CJ,(z) -2DH\(z), (18) 


where J,,(z), H!(z) denote respectively the Bessel function and Hankel 
function of order », and C, D are constants to be determined. For definite- 
ness we take arg pu Sa and |argz| = 47. A general integral of (16) may 
now be seen to be 


> 


, ii 
f(p) = Aip+ = + CI(p) + D1,(p), (19) 


where 


fi Lye, ) 
I(p) =p} J,A2)dp— - | pJ,(z) dp, | 
Jy 
1 (20) 
m ™ 
Ip) = p| H}(z) dp— : | p?H, (2) dp. | 
The constants A,, B,, C, D may be evaluated from the boundary conditions 
(17). Whence 
A, = }—A—3D{I,(1) + 1,(1)}, 
B,=- “ 1DiT,(1) —15(1)}, 


AC (2 ie -(B- 5M) (21) 


AD (2 3— +\1(8)4+(62-1) 2118), 


W here 


— 
A (1 BL (PE) L,(8)fe(1)} 4 
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Thus, when the vorticity of the unperturbed motion is non-zero (A + 0), 
the solution for ¢, is 


> 


I 
ab, = Ap® cosd + QR: | 4 per “1 + CL,(p) + DIse) |e’ . (23) 
Pp } 


where J,, J,, 4,, B,, C, D are defined in (20) to (22); whilst the solution for 
1» 49) 44, Dy 


Gy 18 


¢, = —4&{[C/,(2)+ DH,(2) le}. (24) 


‘The exceptional case in which the motion described by yp is irrotational 
A = 0)arises whengq,/g, = 1/8. Inthis case g(p) has the form 
42/91 5\f 


2(p) =a 2Cp" -~2Dp : (25) 


The function f(p) may be obtained in the same way as when A 40, the 


‘Bessel’ functions ,(z), H(z) being replaced by p", p~". Thus 
en 
f(p) = A; p4 Qs + CL,(p) + D12(p), (26) 


where J/,, /, are now 


The constants A,, B,, C, D are again related to J, and J, as in (21). The 
solution for ys, is therefore formally the same as in (23), the leading term now 
vanishing and the functions /,(p), J,(p) being defined asin(27). ‘The solution 


for C,, however, changes to 
bi — —42#\[Cp" + Dp ‘Je’ \ (28) 


The significance of these solutions is that they express in closed form the 
dependence of a viscous motion on Reynolds number. _ In what follows they 
are used to derive the asymptotic behaviour of the motion when the Reynolds 
number is large. Certain properties of this motion may of course be anti- 
cipated. ‘Thus we may expect viscous action to be negligible save in thin 
layers at the cylinders, and in these layers the Prandtl approximation should 
hold. Further, the vorticity at infinite Reynolds number should be uniform, 
save at the cylinders (Batchelor 1956). Our main interest in deriving the 
asymptotic behaviour of ¢,, however, lies in the higher order corrections to 
the Prandtl approximation, or, equally, in the higher order terms of the 


asymptotic representations of 4. 
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4, BEHAVIOUR OF THE SOLUTION TO THE FIRST ORDER VORTICITY 
EQUATION FOR LARGE R 

As a preliminary to deriving the behaviour of the first order stream 
function y, for large R, it is helpful to consider first the behaviour of the 
vorticity component ¢,._ It is sufficient, for this purpose, to confine attention 
to those properties that follow directly from the vorticity equation (13 a), 
when the way in which the vorticity component ¢, varies periodically with ¢ 
is known. In view of the similarity of the approximate vorticity equation to 
the exact vorticity equation expressed in coordinates for which the stream- 
lines are coordinate lines, these properties may be of interest in suggesting 
general properties of the vorticity distribution for motions with closed stream- 
lines. 

Before proceeding further, it should be noted that the method of expan- 
ding the stream function as a series in y cannot be expected to hold for large R 
if the outer cylinder is at rest. ‘The reason for this may be seen by referring 
tothe boundary layer approximation. Ifthe series were valid and a boundary 
layer existed on the outer cylinder, the speed of the inviscid flow at the outer 
cylinder would be expected to be O(yq,). The thickness of the boundary 
layer there would then be O(ay"1”, RR”). ‘Thus for a given value of y and 
a Reynolds number much greater than 1/y (a necessary condition for the 
boundary layer to exist), the flow quantities of the boundary layer would 
not be analytic functions of y. For instance, the velocity gradient cu,/op 
at the outer cylinder would behave like 

q, y®?R'? x (a function of ¢, independent of y and R). 
A similar objection would apply if the inner cylinder were at rest. From 
here on it will therefore be assumed that neither cylinder is at rest. 

It will also be supposed that the cylinders rotate in the same sense. 
With this provision, the speed of the unperturbed flow does not vanish at 
any point of the fluid and some simplification thereby results in the behaviour 
of the vorticity at large overall Reynolds numbers. 

The general solution to the first order vorticity equation for which ¢, 
has the appropriate dependence of ¢ is given by (24) or, if A = 0, by(28). In 
both equations C, D are for the present to be regarded as undetermined 
constants. ‘The behaviour of ¢, for large Reynolds numbers is thus governed 
(C, D apart) by the symptotic behaviour of J,(z) and H}(z), or, if A = 0, by 
the asymptotic behaviour of p” and p™. 

When the cylinder speeds are such that A 4 0 and B 4 0, the behaviour 
of J,(z), H/(z) when both order and argument are large is required. With 
some minor modification due to substituting for the variables » and z the 
more readily interpreted variables p, A, B and R, the asymptotic representa- 
tions obtained by Debye may be used. An auxiliary variable o is introduced 
which is related to the local speed 


B 
n0l0) = 4 Ap+ 5) 


P /O\12 
by emi (=) dp. (29) 
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The asymptotic expansions of J,,(z), H/(z) for the values which z/ assume 


in this problem may then be written (see Watson 1944, also Olver 1954) 
J (3) ~ ky evF P(R), 


30 
H}(z) ~ k,e~ “(Rie Py — R! 3 | (39) 


where k,, k, are constants independent of p, whose values need not be known 
here, and P(R"*) is an asymptotic series in R' whose coefficients are regular 
functions of p, A, B. 

If the cylinder speeds are such that in the unperturbed motion A = 0 or 
B = 0, the behaviour of the vorticity component depends on the behaviour of 
prvBie dD. p-ve 1 (see (28)) or of J,(4/(— Ri)p), Hi(4/(— Ri)p). The asymp- 
totic representations of these functions are analytically simpler than the 
asymptotic representations of J,(z), H\(z) when A 40, B 40 but are not 
otherwise importantly different. In fact, when A =0 the asymptotic 
behaviour of pv“, p-¥"* is correctly represented by the respective 
expressions for J,,(z), H(z) set out in (30) providing that log p is substituted 
foro. Likewise, when B = 0, the asymptotic behaviour of J,(4/(— Rz)p), 
Hi(4/(— Ri)p) is correctly represented by the above expressions for J,,(z), 
H}(z) providing that p is substituted for o. In either case, the function 
adopted for o(p) is related to the local speed of the unperturbed flow as in (29). 

Thus, when the cylinders rotate in the same sense, the behaviour of the 
first order vorticity component at large Reynolds number may in every case 
be represented in the same asymptotic form : namely, 


é. Par Bix, P(R! 2)ev( Rilo bie P(— R'2)e-v ‘Rie * (31) 


where x,, Kg are constants independent of position, P(R"*) denotes an asymp- 
totic series in R? ,and o(p) is related to the speed of the unperturbed motion 
as in (29). 

It follows immediately from the way in which a is related to the speed O 
of the unperturbed motion, that o(p) is monotonic increasing; the compli- 
cation that might have arisen if the speed of the unperturbed motion had 
vanished at some circle p = constant between the cylinders having been 
avoided by supposing that both cylinders rotate in the same sense. ‘The 
vorticity component ¢, thus resolves into two parts, both of which oscillate 
rapidly with changes in p, the amplitude of the first part decaying exponen- 
tially with increasing distance from the outer cylinder and the amplitude of 
the second part decaying exponentially with increasing distance from the 
inner cylinder. ‘This ‘exponential’ property of the two parts that contribute 
to ¢, will be seen to play a significant role in determining the behaviour of the 
motion at large Reynolds number. 


5, BEHAVIOUR OF THE FIRST ORDER STREAM FUNCTION FOR LARGE R 

We now consider the behaviour for large R of the first order stream 
function ,. It is to be anticipated that the behaviour of %, outside the 
boundary layers at each cylinder will differ from its behaviour inside these 
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layers. Accordingly the asymptotic behaviour of 7, is derived firstly for 
points in the fluid whose distance from the nearest boundary is very much 
greater than a/R~'? and secondly for points whose distance from the nearest 
boundary is O(a/R-"?). 

It is convenient to resolve the solution for y%, into the components 


wp = 2A‘(CI, t DI,)e' : 


) 
ub, = Ap* cos ¢ + 24> (4 1pt =e 
The component #, gives rise to the vorticity distribution ¢, and will be referred 
to as the strongly rotational component, whilst +, makes no contribution to ¢, 
and will be referred to as the weakly rotational component. ‘The strongly 
rotational component may further be divided into two parts, namely 


tp = 2ACI, e*} dpe = 2A{DI, e*}, 


fs 
each corresponding to one of the two parts of ¢, distinguished above. 

For the case where the cylinder speeds are such that A + 0, the strongly 
rotational terms ¢sp1, Ye give rise to the respective vorticity distributions 
2A\CI(z)e®', 2AIDH'(z)e®'. The functions J, and J, which define the 
way in which¢p,,%p2 vary with p, are expressed in terms of simple integrals 
of the corresponding functions J,(z), H'(z), and their asymptotic behaviour 


is therefore related to that of J,(z), H}(s). The detailed definitions of J,, J, 
were given in (20). Whence also 
7°) ] O 


I= | J,(z) dot | p,(z) dp, 
1 Poet 


Cc 


; Pp 1 rp She, 
I,= | H\(z)dp+ | p?H\(z) dp. 
SB p* J 
On substituting for /,(s) and H(z) the asymptotic expansions (30) we find 
that 
I, ~ «,4a,R 4+ > a, Rom +e" Rio 


» 


n 


I, ~ Ko b, R24 > b, Rom? pe- RO, 
are (33) 


lm wy V(io'a, + > at Rothe Rio 


i ~~ Ko} — (1)o’b, - p. b* R-nl2 \e-s i me 


n=S 


where the coefficients a,, a*, b,, b7, are independent of R. As was to be 
expected, in view of the corresponding behaviour of /,(z), both /, and J; are 
smaller at any point p than at the outer cylinder (p = £) bya factor of O(e~*' ni 
where d, = 0(8)—o(p)>0. Similarly, 7, and J; are both smaller at any point 
p than at the inner cylinder (p = 1) by a factor of O(e-¥°”), where 


d, = o(p)—o(1) 0. ‘The asymptotic representations of dp1, %p. and also 
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~ 


ws, now follow from (32) 
C, D. Hence, 


i es aca /(2 LS a, Rome fe-vetniors—aipn +i | (34 a) 
y R1 R! 2° \ igs road a) 5 . c 


and the expressions (21), (22) for the constants 4, B, 
i 


p Ba. : a Ri : 
j M12 - oO n/2 J(Ri ( i i 2 \ 
Une 5 ~ Ris” | - qini T x | 2 R i . (34 b) 
A. aad 
b,~ — —(p*- ] )( 6" — p*) cos d+ 
p 





Balers | (Bs (Yb fo f(A] 
re Tp — te 7 N \ Ye p dh NNGS) 


+S tiwems NR), (34c) 
where «,, 8,, are independent of R, and /,, m, are independent of Rand p. 
It is evident from (34), that the constants A, B, C, D are such that the 
strongly rotational components ¢/p 1, %pg and the weakly rotational component 
ys, are each O(R-"?) at the bounding cylinders. ‘The fact that the components 
bri Ue are transcendentally small save in thin layers of thickness O(aR "”) 
at the inner and outer cylinders follows automatically from the ‘ exponential ’ 
behaviour of these components. 

In the special case where A = 0 the analysis may be retraced with only 
small changes, caused by replacing J,(z), H(z) by p", p-“; and, on putting 
log p for o, the results (34) are obtained unimpaired. ‘The asymptotic 
representations (43) therefore hold for all relevant values of A. 

The asymptotic representation of the first order stream function y, 
appropriate to points whose distance from either boundary 1s much greater 
than aR-!? is evidently that given by (34) for %,. For points whose distance 
from either cylinder is O(aR-!?), the strongly rotational component is signi- 
ficant. Consider the flow near the inner cylinder. ‘The asymptotic 
representation of %, for a point p = 1+€R*, where € is independent of 
Reynolds number, clearly includes %, and yp. but not yp; Whence, from 
(34 b) and (34 c), we find after some rearrangement that for fixed €, 


Z Bs P Bq a 
— FR « ais Jouibienede &€ a — on ETS mu 
~ Re L( ae) qi! ) | 


+ > AP, e IR2+ FY AO, &“YR"2, (35) 
where P,, and O,, denote polynomials in € of orders n and n—1 sar seen 
In the same way the asymptotic representation of %, fora point p = B- -_ *, 
where 7 is indepe ndent of R, near the outer cylinder includes , and &,, but 
not ypo. Thus, from (34 a) and (34 c) we find that for fixed 7 
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where S,, and 7,, denote polynomials in y (different from before) of orders n 
and n—1 respectively. 

The main features of the motion described by these expansions are two- 
fold. ‘The first is that the stream function 4, may be expressed both outside 
and, with suitably rescaled coordinates, inside the boundary layers as an 
asymptotic series in R'?, ‘The zero order stream function yy is independent 
of R, so that similar expansions hold for %)+y,. ‘The second feature is 
that the vorticity ¢, associated with y, is zero, to the accuracy of the 
asymptotic representation for y,, in the region outside the boundary 
layers. ‘The vorticity ¢) is constant, so that the vorticity ¢) ++, associated 
with yw» +y%, is again, to the accuracy of the asymptotic representation, 
uniform throughout the region outside the boundary layers. ‘The leading 
terms of the expansions (34 c), (35), (36) for %, are those that would be obtained 
by the use of the boundary layer approximation. ‘The leading term of y,, 
in particular, represents the first order perturbation to the inviscid motion 
(that is, the flow at infinite Reynolds number). As expected, the speed of the 
inviscid flow is non-zero at the bounding cylinders being of amount 2fyq, cos ¢ 
at the inner cylinder and 2yq, cos ¢ at the outer cylinder. 

The behaviour of the stream function ys. + ys, outside the boundary layers 
is typical of more general closed motions. It was to be expected that yy + yy, 
would be expressible as an asymptotic series of the form > #” R-"? where the 
coefficients i" are independent of R, and further (Batchelor 1956) that the 
inviscid motion would have uniform vorticity. It now appears that the 
vorticity associated with the higher order terms of this expansion is also 
uniform. With the general assumption that the stream function behave at 
large Reynolds number R like >” R-"?, where the coefficients ¥” are 
independent of R, this property is true of any two-dimensional motion whose 
streamlines are closed and lie in a domain free of shear layers. ‘The proof 
of this result entails only a trivial extension by induction of the argument due 
to Batchelor, and will not be given here (though it may perhaps be remarked 
that in the iteration, the fact that the vorticity associated with ¥ is uniform 
follows from the requirement that the pressure distribution associated with 
i") should be single-valued). 

In the present problem the vorticity €)+y¢, comprises a uniform 
distribution 2A and a distribution which outside the boundary layers 
becomes transcendentally small (i.e. o(R~”), for arbitrarily large n) when the 
Reynolds number is large. ‘This latter distribution, it will be recalled, 
divides into a part which decays ‘exponentially’ rapidly with increasing 
distance from the inner boundary and a part which decays ‘ exponentially’ 
rapidly with distance from the outer boundary. Further, each of these 
parts represents the vorticit: distribution, apart only from a constant term 
and a term which is transcendentally small, in the appropriate boundary 
layer. By regarding the vorticity distribution as a uniform field with two 
‘exponentially decaying’ fields superposed, we gain a picture of the distri- 
bution which comprehends at once both the behaviour of the vorticity in 
the boundary layers and the result that the vorticity is uniform, apart from a 
transcendentally small term, outside the boundary layers. 
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The concept of vorticity decaying so rapidly with increasing distance from 
the generating surfaces as to become transcendentally small outside a thin 
layer is familiar, and, in the boundary layer approximation, fundamental. 
The possible occurrence of a component in the vorticity distribution which 
does not vary rapidly with distance is less well recognised. ‘The occurrence 
of such a component is, however, a priori to be expected. For any solution 
of the vorticity equation for which the diffusion term vanishes identically 
may supply a component which does not change rapidly when the Reynolds 
number is large. In the case of a two-dimensional motion with closed 
streamlines whose stream function can be represented as above in the form 
> 4 R? it appears that the only possible slowly varying component is one 
for which the diffusion term vanishes. ‘This is also true of certain axi- 
symmetric motions with closed streamlines (Batchelor 1956). 

Concerning the mode of representation of the stream function yf, it is 
interesting to note that the asymptotic expansions (34 c), (35), (36) may be 
embraced in a single representation which holds uniformly throughout the 
fluid. Such a representation is achieved by writing for %, its asymptotic 
series (34 c) and for &p1, pe their ‘ boundary layer’ representations and then 
combining all the terms. ‘Thus, we put 

P hiro 2p 
(p* — 1)(B? — p*) cos 6 + Ri 782 —1) ve 
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a > ANT, eie—v (ide Bay)n\ R- 24 o(R™ 2), (37) 
n=2 


In the boundary layers, the weakly rotational component may be re-expressed 
in re-scaled coordinates, and (37) then transforms into (35), (36). Outside 
the boundary layers, the strongly rotational component is transcendentally 
small, and (37) then reduces to (34c). A similar representation holds for the 
stream function %)+ys,. Both representations correspond with the division 
of the vorticity field into a uniform component and two ‘exponentially 
decaying’ components. 


6. "THE HIGHER ORDER APPROXIMATIONS OF BOUNDARY LAYER THEORY 


In the following, an expansion for ¢, which holds uniformly throughout 
the flow region is calculated by the methods of boundary layer theory. For 
simplicity attention is confined to the case in which A = 0. The key to the 
procedure to be used is the knowledge that ys, can be expressed as in (37). 
Lagerstrom & Cole have proposed that a stream function may generally be 
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expressed (flows with separating boundary layers excepted) as a sum of 


components of the form > %,(2,,2.)R-"? and } &,(%,, %)R?, where 
21, S) are coordinates such that , vanishes on the solid boundaries, and Risa 
Reynolds number. ‘The second component is further supposed to be 
transcendentally small in the region outside the boundary layers. Once 
it has been established, or can be regarded as established, that the stream 
function may be expressed in this form it is not difficult to devise a procedure 
for calculating the higher order corrections to the Prandtl approximation. 
The aim in what follows is to illustrate this procedure for a motion which is 
closed. 
We begin, then, by assuming that for large R, 


py = 2 JI ,(p, p)f sale > BE, ¢)R-"? T 


n=( n=1 


a: > BAN, bt w+ o(f od (38) 
t=] 
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the terms R,,,(€,¢) being supposed transcendentally small for large é 
(i.e. o(€-"), for any 7) and the terms R,,.(y, 6) being supposed transcendent- 
ally small tor large 7. ‘The leading summation represents the (truncated) 
asymptotic series for ys, appropriate to a fixed point in the fluid. 

As was remarked in the preceding section, the vorticity associated with 
the terms 7, must be uniform throughout the motion. ‘Thus 


V2s, =5,, (39) 


n 


where the 6,, are constants independent of position. ‘lo obtain equations for 
R,,(€,¢) we consider the boundary layer at the inner cylinder. ‘The 
contribution of the terms -#,,. to %, may here be neglected. On expressing 
the variables that occur in the equations of motion (13) in terms of €, ¢ and 
then expanding them in powers of R-!?, we get 











———— a = Q) ‘ 
0&4 a 0&7 db 
4G 93 2? ABZ 2p 2€ A3g \ 
Ay 1% _ _2OAy LP UHAy (40) 
0&4 o” 0&0 of og a3 C&0d 0 5 00d : 


where « = 1. ‘The equations for #,,.(n, 6) may be obtained in a similar way 
and, in fact, the first two equations are correctly given by (40) if Ry, Ro, € 
are replaced by Rj», Ros, — 7 respectively and « is put equal to 8. It may be 
noted that these equations for R,,,, R,,. contain no terms from the weakly 
rotational component of %,. "This would not be generally true. In general, 
the weakly rotational component will contribute to the convection velocity 
in the vorticity equation (though, of course, it will not in general contribute to 


the vorticity derivatives). 
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The boundary conditions for the coefficients -7,, of the weakly rotational 
component are 


JF (1,¢4)=0,  -%,(8,¢) = constant, 





> (41) 
J (1, ¢)+Z,,(0,¢) = 9, %,(8,4)+Z,,(0,¢) = constant (nm > 1), } 
whilst the boundary conditions for the coefficients #,,;, 2, are 
oR (0, 4) OF (1, 4) 
vt ig. fe ie it a 2 fo Ss a 0 
aE i cosd (n ) | 
= @ (n> 1), | 
oR 9 0, ¢ oF B, | 
PAnl0, #) - oF n1(B #) =2cosd (n=0), | 
on op : 
| (42) 
= @ (n > 1), | 
‘ .. OA Ae. 
lim 2,,(€,¢) = lim ils ¢) = 0, 
G-> f—> O¢ 
. 08, n, 
lim 2,(7,¢) = lim nal #) 0. | 
n> @ "> F 07) 


J 


In addition, each of the %,, #,, and #,, and the associated pressure distri- 
butions must be periodic functions of ¢. ‘The conditions on #,,,, #,,9 at 
the outer edges of the boundary layers result from supposing that these terms 
are transcendentally small outside the boundary layers. ‘The boundary 
conditions for the -%, derive from the requirement that the normal component 
of velocity at the boundary associated with these terms should annul the 
normal component of velocity at the boundary associated with the terms 
#4, A. The boundary conditions for the #,,,, 7, arise from adjusting 
the tangential component of velocity associated with the terms -7,, to the 
tangential velocity of the boundary. 

The iterative solution of these equations is complicated by the fact that 
the equations that govern -7, ((39) and (41)) are incomplete. ‘That is, if 
#,,, #9 are supposed known, -%, is undetermined to the extent of two 

and represents the (uniform) 
say, arises because the fluid 


constants. One of these constants is 4, 
vorticity associated with .%,. The other, I’, 
encircles a closed boundary, and may be defined as the circulation round the 
inner cylinder of the motion associated with -7,,. 

In the present problem it is a priori obvious from the symmetry of the 
configuration that 6, and I’, vanish. (‘This may be seen by considering the 
motion in which the sign of y is reversed.) ‘The indeterminancy in -/,, is 
thus resolved immediately. ‘To resolve the indeterminancy without 
recourse to a symmetry argument, it is necessary to pass to the equations 
governing #,,.44, B,.1. On solving for #,.,, a condition on the 
tangential velocity distribution 0-7 ,,(1,4)/cp emerges. Similarly, on solving 
for #,,,, 9a condition emerges on the tangential velocity 0.4 ,(8,¢)/ep. From 


these two conditions, we may then deduce that 5,, = 0, [’,, = 0, and so on. 
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The equations for 4,,, Z,,,, 4,2. may otherwise be solved without further 
difficulty. For the leading terms we get 
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when these terms are combined together as in (38) they yield an expression 
for , which holds to O(R~) uniformly throughout the fluid. 


7. ‘THE TERMS OF HIGHER ORDER IN y 
For the sake of completeness it is worth noting that the higher order terms 
of the perturbation of the motion may be solved in essentially the same way as 
the first order terms. 
When the equations governing the higher order terms €,, %,, of the 
expansions for the vorticity component ¢ and stream function ¢& are examined, 
it is seen that each varies with ¢ according to the respective forms 











- Ea > En, n—2s ex -, fy, = > 4. n—2s el(n : —, (44) 
8 a 8s o 
The equations governing g,, ., f,,, are then 
Py RBsi +s? 5 
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y) 1 ; 2 
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where G,, ,, /,,, are defined by 
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> Fi. n—2s(pye*s* = 4p cos V7xs,_, + 2p?(2 + cos 24) V7,» + 


+ 4p? cos A6V7us,,_5 + p*V7us,,_4. 
Further the boundary conditions are 
fn(1) = 9, fro(8) = constant, f,,(8)=9 (s40), 
fil) =(=1)", fiB) = ~(-B)" (ila) (8 = £m) $47) 
fill) =0, fi(B) = —(-B)-(B2-1) (si<n), | 
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together with the condition that each term of the expression for the pressure 
is singled-valued. 

These equations for the component functions f,, ,(p) of the higher order 
perturbations y,, are similar to the equations for the component term f(p) of 
the first order stream function y,, and may be solved in a similar manner. 
The only features which call for special treatment are the interaction terms 
F,, (p) and G,, (p) and the components of the vorticity that do not vary with 
¢. ‘The interaction terms in the equations for f,, , depend on the components 
f,.; for r = 1, 2...n—1, and are determined iteratively. ‘To determine the 
components of the vorticity that do not vary with 4, it is necessary to use the 
pressure condition. It is found on examining the momentum equations 
that the necessary and sufficient condition for the pressure to be single- 
valued is that 


P fee 
Sno a p | PGing dp, 
whence ae eee P . ; 
£,,9 = constant+R| - | p’'G, o(p’) dp’dp. (48) 
JEP SS 


The corresponding term of the stream function, namely /,, 9, may now be 
determined in the same way as the other terms f, ,, s 40. 

It is not proposed to investigate these equations further. Enough has 
been set down to show that the solution may be completed in essentially the 
same way as for the first order terms. (The homogeneous form of (45 a) now 
has the general solution 2C/,*(s*)+2DH)}«(z*), where n* = 4/(— RAsi,) 
3* = ,/(RBsi+1) and 4 +0). Moreover the solution may be expected to 
behave at large Reynolds numbers in much the same way as did the first 
order terms. In particular it may be tentatively inferred that for large R, 
the vorticity may be expressed as the sum of a constant term > {5,9 y”", 


n=0 
and an infinity of terms which vary exponentially rapidly with R. It may 
also be tentatively inferred that the stream functions ¢,, may at large Reynolds 
numbers be expressed as in (38). If this proves true, the expansions of #, ¢ 
in powers of y may be expected to converge uniformly with respect to R, 
when y is sufficiently small. 


The author is grateful to Dr I. Proudman and Dr H. C. Levey for many 
helpful criticisms of this paper. Acknowledgment is also made to the Chief 
Scientist, Department of Supply, for permission to publish. 
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The supersonic motion of an aerofoil through a 
temperature front 


By J. W. CRAGGS 


King’s College, Newcastle-upon-Tyne 


(Received 28 July 1957) 


SUMMARY 

A solution is given for the sound waves generated when an 
infinite wedge passes at supersonic velocity through the plane 
interface between different media. ‘lhe results are then applied 
to the particular case when the two media consist of the same 
perfect gas at different temperatures. It is observed that the 
flow can be deduced for any symmetrical aerofoil of infinite span 
by superposition, and the case of a double wedge is considered 


in detail. 


INTRODUCTION 

When an aerofoil passes through the surface separating two media, 
as for example into a cloud, or through a temperature front, certain transient 
effects on the lift and drag are to be expected during the transition from one 
regime of steady flow to another. ‘The transient effects will be examined 
by the method used in an earlier paper (Craggs 1956) on the reflection of 
sound waves. 

The method used is first to examine the motion due to an infinite wedge 
moving with constant velocity. ‘The fluid motion may then be expected 
to possess dynamic similarity, the velocity distribution being the same at 
all times after the wedge has broken the surface, except for a scale factor 
depending on the depth of penetration. ‘The solution for the wedge can 
then be extended to any symmetric aerofoil by the principle of superposition, 
within the limitations imposed by the linearized theory. It is assumed 
throughout that the fluid motion is sufficiently small for linearized theory 
to be adequate, and that the motion is strictly two-dimensional. 


1. THE EQUATION OF SOUND WAVES WITH DYNAMIC SIMILARITY 
Use polar coordinates (7,4) based on the point where the apex of the 
wedge cuts the surface of separation of the two media. Let ¢ be the time 
measured from the instant at which the apex breaks the surface. Let c be 
the local velocity of sound and s the condensation, measured by the 
proportional increase in density due to fluid motion. ‘Then the fluid 


velocity q is given by 


4 ~ CVS, (1) 
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and the equation of motion is 
c8V%s = d*s/0#*. (2) 
Make the assumption of dynamic similarity and use A, 8 as independent 
variables, then (2) becomes 
A?(1 — A?/c?)d2s/0A? + A(1 — 2A?/c?)es/0A + 0s /00? = 0. (3) 
Equation (3) is the same as the equation arising in the linearized theory 
of steady three-dimensional flow when the cone-field method is used 
(Ward 1955, p. 136) and the treatment here is similar. 
For A > c equation (3) is hyperbolic and can be reduced to the canonical 
form 


os os 
> Ape — 0 4 
op” og? ( ) 
by the substitution secu = A/c. (5) 


The general solution is then of the form 
$= f(u—6)+g(u +8) (6) 
where f and g are arbitrary functions, and particular cases can be solved 
by the method of characteristics. 
Equation (3) is elliptic for A <c and is reduced by the substitution 


A/c = sech(— v) (7) 


to the form 25 ABs 


w= s+ir, (9) 

then from (38) w= w(v+ie) (10) 
and Os CT os : OT 
ay 00? 0@ 01 


2. ‘THE SOLUTION IN THE REGIONS WHERE THE EQUATION IS HYPERBOLIC 

For simplicity, assume that the wedge moves in the direction normal 
to the surface separating the two media. Let the motion be in the direction 
6= 47, and let p,, pg be the densities in the regions 0 > 6 > —7a and 
0 <6 <7, respectively. Use the suffixes 1 and 2 consistently to refer 
to the two media. Consider a plane with polar coordinates A, @ (figure 1). 
Then in the regions 0 < 6 < —z7, A>c, and 0 >@> —-7z, A> ¢, the 
equations are of hyperbolic type, the characteristics being the tangents to 
the sonic circles, A = cs, A = ¢,, respectively. 

Consider ‘first the motion of the wedge with velocity V > c, before 
the apex reaches the point r= 0. ‘Then the motion is the familiar steady 


F.M. M 
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motion of a wedge with supersonic velocity and the solution is well known. 
There is a weak shock wave at the Mach angle ¢, = sin“‘(c,/V) to the 
direction of motion, attached at the apex of the wedge. In front of the 
shock wave the velocity of the fluid is zero and the condensation is also 
zero. Behind the shock wave the condensation is 

$1] = 0, = M, Ssec dy, (11) 
where M, is the Mach number, M, = V/c,, and 6 is the semi-angle of the 


wedge. 











Figure 1. 


To discuss the motion after the wedge has broken the surface, measure f¢ 
from the instant at which the apex reaches the surface. ‘Then the condition 


of steady motion gives the motion as t>0+, A+ o. ‘The boundary 


conditions on the equations as A > © are therefore 

$1; = 9}, —$7 <O0< —}r+q, | 

— | —~Int+¢, <0<0, (12) 
and s, = 0, 0<9@< fn, J 
with similar conditions for $7 < @ < 3(37). 


Use of characteristics theory shows that there is a weak shock wave, 


of trace 
A = c,sec(d, — 4) (13) 


and of strength o,. ‘This is represented by the line HJ in figure 1. 

It is convenient to assume that c,sec¢, > ¢,, so that the point H, 
where the shock wave cuts # = 0 lies outside of the circle A = c¢,. (Only 
slight alteration of the argument is needed when this assumption is false, 
and for the application contemplated, in which the difference in properties 
between the two media is small, the assumption certainly holds.) Further 
simple deductions from the theory of characteristics are that there are 
weak shock waves HK, HL corresponding to the reflection and refraction 























The supersonic motion of an aerofoil through a temperature front 179 
of the incident wave H./ at the interface, and a new shock wave FG attached 
at the apex of the wedge. ‘The condensation behind FG is 

Sy = 0g = M,Ssec 4), (14) 


where ¢, = sin-\(c,/V). ‘The conditions behind HK, HL can be deduced 
from the shock wave equations and the continuity of pressure 


Pret Sy = P2352, (15) 

and of normal velocity cj 08,/00 = c50s,/08, (16) 
over the interface 6 = 0. ‘The values are 

5, = 20, tan d,/{tan d, + k(m? sec? d, — 1)"?! = 6} (17) 

and o, = mko,, (18) 


respectively, where m = ¢,/C, k = p,/po. 
The condensations are now known throughout the regions A > ¢, 
\ > ¢y except for the triangular region ABE, and in this region the theory 
of characteristics shows only that 
5, = 0,+f,(9—p), (19) 
the function g of equation (6) being here constant and equal to oy. 


3. REGIONS WHERE THE EQUATIONS ARE ELLIPTIC 


To discuss the region 0 < @ < 3a, A < cy in which the equation for 
the condensation is elliptic, it is convenient to introduce a transformation 
Co = €g +g = sech 2(v, +70) (20) 


which maps the region on the upper half of the € plane. ‘The closed boundary 
of the elliptic region maps into the axis 7, = 0, and the points O, A, B, 
L, G, C of figure 1 correspond to the points 
{, = 0, m?/(2—m?), 1, m? sec? ¢,/(2 — m sec” d), sec? 6.(2 — sec? d.), —1, (21) 
respectively. 

Now consider the section AB of the interface. On the lower side of AB 
s, is governed by a hyperbolic equation, and from (19) 


ra § Os 2 1/2 °C S 
i a — — A — 1 ) me ; 
06 a Ff @ 








Op. 1 
using (5). The continuity conditions (15) and (16) therefore give 
\(1 2 V2 OS> } OSo C T9 P (1 A2\ 126 To 79 
MLA a = RS = — RS HS Se ao —. LL 
( ) oA ot OV 4 oA ( ) 


using (7). ‘This relation holds on the line AB, and on this line 


2 ft ((2—m2)E,-1\12 
a dw, {y_ Te) rF(L.) (23) 


df, | mk 


on the line, where F(,) is real on the line, and by the principle of analytic 
continuation equation (23) holds throughout the region 7, > 0. 


M2 
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The function F(¢,) is to be determined subject to the following 
conditions. 

(i) Singularities of F are to be expected only at the points (21) and at 
at any zero of the function in brackets in (23). 

(ii) For finite pressures, no singularity of dw,/d¢, may be of higher 
order than a simple pole. 

(iii) The point at infinity in the ¢, plane is an ordinary point, since it 
corresponds to an ordinary point, 6 = }7, A = c, in the (A, @) plane. 

(iv) From the conditions on the quarter circle BGC, 


ds, 
a," 
for 7, = 0, —-0 < & < —1,1<& < o. 


(v) The normal velocity on the surface of the wedge is constant, so 
d7,/00 = ds,/eA = 0. On the corresponding section of the real axis yn = 0, 
—1 < &, < 0, the condition 0s,/0¢, = 0 holds. 

(vi) A careful examination of the conditions in the immediate 
neighbourhood of the point A of figure 1, with use of the continuity 
conditions across the interface, shows that F(¢,) can contain only half- 
integer powers of the quantity (2—m?)¢, — m*. 

A sufficiently general form for w, is then given by 
dw, [| mk(1— Cy)? +1{(2—m®)¢,— m?}'? 

- | m*k?(1 — €,) + (2 —m?)¢, — m? |. 























with a B - Cc (25) 
2 (2—m* sec? d,)6.—m? sec? d, (2—sec” do)C, — sec” dy ie 

and O, = D - E (26) 
<"  (2—m* sec? g,)C,— msec? d,  (2—sec® do), — sec” dy" | 

For the solution in the region A<c,, —da 86< 0, a similar 


transformation 
¢, = sech 2(v, +20) 
is used, and (24), together with the continuity conditions across the 


interface OA, leads to 


dw, kil —(2m? - 1)¢, 2? -(1-¢,)¥? 
mk 7 — <= ay — — 7 x 
eG k?}1 —(2m? —1)¢, —(1— 4) 


[aoe + (2m? = (27) 








m*k(1 + ¢,)"? 














with Pp Lb mC (28) 
1 (2—sec? ¢,)¢,—sec? 6, (2m? — sec? d,)¢, — sec” dy ioe 
ind D mE 
: O (29) 


> ; 5 ai > ; ; : . 
<1 (2—sec? 4,)6,—sec?d, — (2m? —sec? d)C, — sec? dy 
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The constants B, C, D, E are easily determined by an examination of 
the singularities at the points corresponding to G, L, K of figure 1. At G, 
for example, where ¢, = sec? ,/(2—sec? 4.), there is a discontinuity of 
magnitude o, in s,, whereas at the corresponding point in the ¢, plane, 
f, = sec? d,/(2m? — sec? d,), there is no singularity. Applying the theory of 
residues in each plane, one obtains 
C+~2E(sec? d.— m?)'* cos dy 





ah = mk tan $, — (sec? dy — m?)!” 
and C—~2mkEsin ¢, = 0, 
whence E = V2mkoy sec $2 (30) 
and C = 2m*k?o, tan do. (31) 


A similar argument based on the singularity in the ¢, plane at 
f, = m* sec? ¢,/(2—m? sec? d,), corresponding to H, where there is a 
discontinuity o, in s,, and on the singularity at the corresponding point K 
in the ¢, plane where s, changes by (c,— 4), gives 
B= —2m'k’o, tan 4, (32) 
and D = —~2m'k*o, sec 44. 33) 
The functions dw,/d¢,, dw,/d¢, are now determinate, and the solution 
of the problem is complete. The pressure at any point can be obtained 
by integration of (24) and (27) in terms of elliptic functions, but since in 
most cases only the pressure on the face of the wedge will be of interest, 
it is easier to evaluate the integrals numerically as required. 


4, 'THE SPECIAL CASE OF A TEMPERATURE FRONT 

When the two media concerned are both composed of the same gas, 
but at different temperatures, some simplification is possible. 

The gas laws then give k = p,/p, = 7,/T,, ¢?/c3 = m? = T,/T,, and 
mk = 1. It will generally be sufficient to use a first approximation, writing 
the absolute temperature ratio as T,/7, = 1+.« and retaining only linear 
terms in «. Then k= 1+ 6, m* = 1—€« and sec” 4, = sec? $,(1+ tan? ¢)). 
Equations (25), (26) and (30) to (33) then lead to 
4M, dC, sin d, 








2 - ] 
a (C, cos 24, — 1)? (34) 
and 
a wi v2M, de{(1 : $a) + sin? (520s 26, —1)} . (35) 
7 cos* ¢,(¢, cos 24, — 1)? 
and equation (24) gives 
a ds 2asind (1—«)!? 











Mode dz (acos2h6+1)(1—a)!2 (2a)! cos? (1 + ~ cos 24)? 
tan? ¢ 
(2x)"2(1 —«)"9(1 + cos 24) 





+ 


where a= —€ = A?/(2c?—A?) 
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and the sufhxes have been dropped, since the use of M/, for M, introduces 
errors only of smaller order than the terms retained. ‘The corresponding 
result in the other medium is, from (27), 
a ds 2asind . (1—«)!? 
Mie dx (xcos26+1)%(1—«)' | (2x)! cos’d(1 + acos 26)? 
tan” oe 
~ (2x)! 2(1 —«)} a + «cos 2¢) = —-i) 8 
Integration of (36) gives 
Ss—G, sindcos2¢d 2 sin 6(1 — «?)!? 2 sind cos 2d : 


“Mée ~— sin? 26s zsin? 24(1+acos2¢) = zsin®2¢ 




















us gf 08 2¢ 1 (1—«)! ee 
«sin ’\acos26+1)* baa ii cos? (1 + cos 24) 





l ; 1/2) sin? d 
~ 2ncos?d sai cos d(T i 4cosh 
bP a f 1 —a(1+2cos* $)| (39) 
~ 2rcos? d | l+acos26 |’ 


and a similar result follows from (37). 


5. APPLICATION TO GENERAL SYMMETRIC AEROFOILS 
In the linearized theory, superposition of results is allowable, and the 
motion of a general symmetric aerofoil can be obtained by regarding it 

as being made up of a number of wedges. 
* Define a function G(A), — 0 <A < o, as follows. For 0 <A < o, 
the value of 5G(A) is equal to the condensation, sj, in the above solution 
on the line @ = }z at the point A. For — «0 <A <0, the value of dG(A) 


is the condensation at 6 = — $7, X. 
Now consider a symmetric aerofoil defined by 
—X(y)<x« < X(y), —@<y <9, 


where x, y are cartesian coordinates moving with the aerofoil and based 
on its leading edge, and write 

o(v) = —dx/dy 
for the slope of the surface, assumed everywhere small. 

The condensation on the aerofoil before the leading edge breaks the 
surface of separation of the media is the same as in steady motion in medium 1. 
When the aerofoil is part way through the interface the condensation due 
to that part which has passed into the second medium is given by the 
foregoing theory, but the flow due to the part which has not yet left the 

| medium 1 is still the steady flow appropriate to medium 1. Consider the 
conditions at a time ¢, less than a/V, after the leading edge breaks the surface, 
then the foregoing arguments give, on the aerofoil, 


Vit+y\ 
$= WOo(—*)- | (0) — $0") x 
¢ / —Vt 


72 ae ex gr! 
«oP V) aye > ¥> V0), 0 


Vct + cy* 
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and 


Vt+y\ 0 V2t+ Vy—Vy* 
— z{ ——— }— tT —u(y*)'G —— = = 
6 = n6(“2)— f"  6o-sone( ara) a 





— Msec 4, | ; {b(0)—b(y*)} dy* (—Vt>y> —a), (41) 


where V is, as before, the (supersonic) velocity of the aerofoil. When the 
aerofoil has completely passed the interface only the expression (40) is 
required. 


Y 


6. NUMERICAL RESULTS FOR A DOUBLE WEDGE PASSING THROUGH A 
TEMPERATURE FRONT 
As an example, consider the motion through a temperature front of 
a double wedge, of angle 5 and length 20. 




















“1-0 "1-0 
r H | H dx A H | H dn 
1:0 to 1°5 0-1341 0-00000 0-05 0-3978 058093 
0:95 0:4896 0-01559 0-10 0:3775 0:60031 
0-90 0:5867 0:04250 0-15 0-3567 0-61866 
0-85 0:6328 0-:07299 0:20 0-3359 0-63598 
0:80 0:6557 0:10520 0:25 0:3153 0:65226 
0:75 0:6654 0:13823 0-30 0:2948 0-66751 
0-70 0:6636 0-17145 0-35 0:2744 0:68174 
0-65 0:6570 0):20447 0-40 02532 0-69493 
0-60 0-6461 0):22704 0:45 0-2415 0-70730 
0-55 0:6322 0:26900 0:50 0:2139 0:71868 
0-50 0:6167 0-30022 0-55 0-1939 0:72888 
0-45 0-5960 0-33054 0-60 0-1779 0-73818 
0-40 0:5822 0-36000 0-65 0-1535 0-74646 
0°35 0:5623 0:38861 0-70 0-1335 0-75363 
0-30 0:5424 041623 0:75 0-1124 0-75978 
0-25 0:5223 0:44284 0:80 0-0935 0-76493 
0-20 0-5017 0-46844 0-85 0-0738 0-76911 
0-15 —():4810 0-49301 0:90 0:0539 0:77230 
0-10 0-4602 0-51654 0-95 0-0334 0:77449 
0-05 0-4399 0-53904 1-00 to —x 0-0000 0-77532 | 
0-00 0-4188 0-56051 | 








Table 1. 


It is convenient to tabulate a function H(A) defined by 
G A = 
H(A) = 2, 
€O71 
Then (39) leads to the results in table 1. 
The drag en the aerofoil then follows from (40) and (41). In table 2 
a correction factor e is given such that the total drag is 


pci a, bd{2—ee}. (42) 
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The correction factor is tabulated against V’t/b, that is the number of half 
lengths by which the leading edge has passed the interface. 





Vt/b| 4 1 2 3 4 5 6 | 
e |0-2808 0-2752 0-3103 0-5476 0-0817 0-0209 0-1712 02682 | 
| 
| 





Table 2. 


When I’t/b exceeds 6, the transient effects no longer influence the drag 
and the drag is constant and equal to that for steady motion in the new 
medium. 
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SUMMARY 

A mechanism for the generation of surface waves by a parallel 
shear flow U(y) is developed on the basis of the inviscid Orr 
Sommerfeld equation. It is found that the rate at which energy 
is transferred to a wave of speed c is proportional to the profile 
curvature — U"(y) at that elevation where U=c. The result is 
applied to the generation of deep-water gravity waves by wind. 
An approximate solution to the boundary value problem is 
developed for a logarithmic profile and the corresponding spectral 
distribution of the energy transfer coefficient calculated as a 
function of wave speed. ‘The minimum wind speed for the 
initiation of gravity waves against /aminar dissipation in water 
having negligible mean motion is found to be roughly 100 cm/sec. 
A spectral mean value of the sheltering coefficient, as defined by 
Munk, is found to be in order-of-magnitude agreement with total 
wave drag measurements of Van Dorn. It is concluded that the 
model yields results in qualitative agreement with observation, but 
truly quantitative comparisons would require a more accurate 
solution of the boundary value problem and more precise data on 
wind profiles than are presently available. ‘The results also may 
have application to the flutter of membranes and panels. 


1. INTRODUCTION 

The primary aim of the following analysis is the prediction of the energy 
transfer from a prescribed two-dimensional parallel shear flow in an inviscid 
incompressible fluid to a surface wave of prescribed wavelength and wave 
speed. The principal result to be obtained is that this transfer is propor- 
tional to the curvature of the velocity profile at that point in the profile where 
the mean air speed is equal to the wave speed. ‘The formulation, as developed 
in §3—5, is valid for surface waves of a rather general type and may be of 
some importance in the flutter of thin panels or membranes, but we shall 
focus our attention primarily on two-dimensional deep-water gravity waves. 

We shall not attempt to develop the mechanism by which the surface 
wave energy is dissipated except in the initial formation of gravity waves 
on a previously smooth surface, and even in this case it will be necessary to 
introduce assamptions which would appear more questionable than those 
implied in the original model for energy transfer from the air to the surface 
wave, 
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Previous models of wave formation have been critically reviewed by 
Ursell (1956). ‘The mathematical (as opposed to empirical) models may be 
briefly characterized as follows. 

(1) The Kelvin-Helmholtz model (Lamb 1945, § 268) assumes a tangential 
discontinuity between uniform flows in air and water and predicts a 
critical wind speed of 650 cm/sec for the initiation of surface waves. 
Jeffreys’ (1925) sheltering model (Lamb 1945, §348) assumes a 
periodic component of wind stress in phase with the wave slope, 
which is supposed ot originate in the separation of the air flow over 
the wave crests and is described by a ‘ sheltering coefficient’. No 
method of calculating this coefficient is given, however, and the form 
of the results would be the same for other mechanisms. 

(3) The laminar flow models of Wuest (1949) and Lock (1954) assume 
laminar flow in both air and water and yield the critical wind speed 
as in the corresponding analyses of stability for channel and boundary- 
layer flows. 


—_~ 
bo 
— 


(4) Eckart’s (1953b) stochastic model assumes an ensemble of gusts 
having prescribed durations in time and space. The details do not 
permit brief summary, but it should be noted that the assumed wind 
structure differs widely from that considered here. It may approxi- 
mate a storm more closely, but it would be far more difficult to 
reproduce in a controlled experiment. 


The model to be developed here improves on the Kelvin-Helmholtz 
model by allowing for distributed, rather than concentrated vorticity*, 
resembles Jeffreys’ model in predicting an equivalent value for his sheltering 
coefficient, and is similar to the investigations of Wuest and Lock inasmuch 
as it leads to a boundary-layer stability problem. We define this model by 
the following assumptions. 

(a) The air. We assume the air to be inviscid and incompressible and, 
in the absence of the wave motion, to have the prescribed mean shear 
flow U(y). ‘The subsequent disturbances in velocity and pressure 
associated with the wave motion will be assumed two-dimensional 
and sufficiently small to justify the linearization of the equations of 
motion. ‘Turbulent fluctuations, albeit decisive in maintaining the 
mean shear flow, are neglected in the perturbation equations. The 
end result of these assumptions is the inviscid Orr-Sommerfeld 
equation (3.4). 

(b) The water. We assume the motion of the water to be inviscid, 
incompressible, and irrotational, and the slope of the displaced 
surface to be small. Following Jeffreys (1925; Lamb 1945, § 348), 
the dissipation in the water will be calculated in the second approxi- 
mation as a perturbation on the inviscid flow, assuming laminar 


* The assumption of a concentrated vortex sheet in the Kelvin-Helmholtz model 
was criticized originally by Rayleigh (1880; also 1945, § 366). 
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motion and a free surface, but this result will appear only in the 
calculation of minimum wind speed (§6)*. We shall neglect any 
mean motion of the water that may be induced by the mean air flow. 
This neglect could be justified if either the actual mean velocity 
were small compared with the wave speed (so that it could only 
produce negligible inertial forces) or were laminar and confined to 
a layer thin compared with the wavelength. 

(c) Wave speed. We assume that the access of inertia associated with 
the disturbance in the air has a negligible effect on the magnitude 
of the surface wave speed, only that component of the aerodynamic 
force that is in phase with the wave slope being regarded as important. 
‘This assumption, which is introduced only to simplify the algebraic 
operations by virtue of an expansion in the specific gravity of the 
air (and which is identical with that made in Jeffreys’ sheltering 
theory), rules out the Kelvin-Helmholtz mechanism of wave 
formation a priori. 

Referring to the assumption of an inviscid fluid, an appropriate Reynolds 

number for travelling waves in a slightly viscous fluid is 


R = |U—c\/kv = |U—c|A/2av, (1.1) 


where U, c, A, v, and k denote the speed of the undisturbed flow, wave 
speed (phase velocity), wavelength, kinematic viscosity, and wave number. 
R may be evaluated for either air or water, but it will attain its smallest values 
in the air, both because of the higher kinematic viscosity and because U —c 
may tend to zero. In the examples treated subsequently (§6, §7) R will 
have values of the order of 10% at the outer edge of the boundary layer, 
though it tends to zero at the critical point where U = c. ‘This implies 
that viscous forces in the air are likely to be important only in the neighbour- 
hood of U = c, where their omission leads to a singularity in the equations 
of motion ($3). Previous studies of channel and boundary-layer flows 
(Lin 1955, ch. 5) suggest that this singularity provides an adequate repre- 
sentation of viscous effects (on growing waves) for sufficiently large average 
values of R. We specifically remark, on the other hand, that (in contrast 
to the Poiseuille or boundary-layer stability problems) the perturbation 
inertia forces for sufficiently large values of ck/v may be expected to dominate 
the perturbation friction forces in the immediate neighbourhood of the 
boundary by virtue of the surface wave there. This implies, in particular, 
that the viscous drag forces of the air are negligible compared with the 
normal pressure in their effect on the surface wave (cf. preceding footnote). 

* It might be thought that the air, behaving as in the boundary layer over an 
oscillating plane (cf. Lamb 1945, § 351) and having the higher kinematic viscosity, 
could contribute more damping than the water acting as a free surface (Lamb 1945, 
§ 348). The ratio of the damping coefficients for small, laminar motions is found 
to be (u,/4u,,)[c4/2gv,(c-U)]!, where ¢ denotes wave speed, v kinematic viscosity, 
ft true viscosity, and c-U an appropriate average through the region of unsteady 
viscous influence. This ratio is about 0-1 in the calculation of the minimum wind 


speed (§ 6). 

















188 John W. Miles 


The foregoing arguments appear sufficient to justify an inviscid model 
for the determination of a first approximation to the disturbed motion of 
the air and the consequent energy transfer to the wave, but they throw no 
light on the perturbation Reynolds stresses associated with the interaction 
between turbulent fluctuations in the original and perturbed flows (see 
Appendix). It is implicit in our model that these also are negligible com- 
pared with the perturbation inertia forces except in the neighbourhood of 
U=c. It would be difficult to provide an adequate a priori justification 
of this hypothesis, but it appears reasonable for gravity waves in consequence 
of their relatively long wavelength. 

It perhaps may be regarded as obvious that, on the one hand, models 
such as we consider here represent a gross simplification of the real problem 
of generation of surface waves by wind and cannot be expected to yield 
quantitative results except under very special circumstances ; but that, 
on the other hand, our present very limited knowledge of turbulent flows 
makes substantial theoretical progress unlikely without appeal to simplified 
models. ‘The ranges of validity of such models must be established primarily 
by comparison with experimental data; this is scarcely available in sufh- 
ciently detailed and reliable form to warrant firm conclusions at the present 
time, but it seems unlikely that any single model (for which mathematical 
analysis is feasible) will prove adequate for all circumstances*. 


2. EQUATION OF MOTION FOR SURFACE WAVES 
We consider first the equation governing the propagation of a small, but 
otherwise rather general type of surface wave that gives rise to an aero- 
dynamic pressure p, acting on the surface. Let m be the effective mass per 
unit area and L a linear operator such that Ly gives the stress resisting a 
deformation 7 of the surface ; then the equation of motion reads 


Ly+mny = —Ppua- (2.1) 
We assume the surface wave to have the form 
als, 8) = af, (2.2) 


where a denotes the amplitude, & the wave number, ¢ the phase velocity, 
x a coordinate measured in the direction of propagation, and ¢ the time; 
here and subsequently, following the usual convention, the real parts of 
complex quantities are implied in the final interpretation. We also assume 
the aerodynamic pressure to be represented by 


Pa = («+ i8)p,,U? kn, ot 


* ** How does wind acting on water give rise to waves ? This is a question which 
has never been satisfactorily answered in detail. This is partly because we do not 
know the exact constitution of an ocean wave, nor of the wind, nor the exact way in 
which the one acts upon the other. Notwithstanding our ignorance of the details, 
the general nature of the operation has been fairly well made out.’”’ (Durand 1896). 
Ursell, in his (1956) survey of the problem, opens with the statement that ‘‘ wind 
blowing over a water surface generates waves in the water by a physical process 
which cannot be regarded as known’ and concludes that “ the present state of our 
knowledge is profoundly unsatisfactory ’’. 
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where p, denotes the density of the air, U, an as yet arbitrary reference 
speed for the air, ky (in magnitude) the local slope of the wave, and «+78 
a dimensionless pressure coefficient; in general, « and f are functions 
of both c and k which depend on the solution to the aerodynamic boundary 
value problem (§ 3). Substitution of (2.2) and (2.3) in (2.1) yields 
Ly —mk*c*n = —(%+78)p, UZ kn, (2.4) 
which is an eigenvalue equation relating c and k for the assumed wave motion. 
The operator L may be eliminated from (2.4) by referring it to the free 
surface wave speed c,, in the absence of the aerodynamic pressure. In the 
latter event (2.4) reduces to 
In = mk?c*,n, (2.5) 


which implies that mk*c%, is the eigenvalue of the operator L. Substituting 
(2.5) in (2.4), we may place the result in the form 


c? = c2+s(a+78)Uj, (2.6) 
where s denotes the relative mass parameter 
S = p,/mk. (2.7) 


We emphasize that (2.6) does not represent an explicit solution for ¢c inasmuch 
as # and f exhibit an implicit dependence thereon. 

We turn now to the special case of gravity waves on deep water of density 
p,- The effective mass then is p,,/k, while the operator L is simply p,,g, 
and we have (cf. Lamb 1945, § 228, § 229) 


~ = gin, S = pp... (2.8 a, b) 
We remark that in the Kelvin-Helmholtz model (Lamb 1945, § 268) 
4+18 =—1 if U, is chosen as U—c, where U is the (constant) air speed; 


in this approximation the aerodynamic pressure is directly proportional 

but 90° out of phase with, the wave slope and represents an access of 
inertia. In Jeffreys’ model (Lamb 1945, § 348), on the other hand, only the 
component of aerodynamic pressure in phase with the wave slope is 
considered, and f corresponds to his sheltering coefficient if U; = U—ce. 

We shall attempt a solution of the eigenvalue equation (2.6) only for those 
values of c and k that render the magnitude of the third term small compared 
with c7,, on which hypothesis we may take c = ¢c,, as a zeroth approximation, 
evaluate «+73 with ¢ = c,,, and take 


c = c,,{1+4s(a+iB)(U,/c,)7] (2.9) 


as a first approximation. If we define the negative damping ratio ¢ as the 
energy rate of growth per radian*, the first approximation to C, (a denoting 
air) is given by (now Seni the subscript w from c¢, as is permissible in 
this approximation) 

C, = 2K {c}/Bfc} = sB(U,/eF. (2.10) 
Our primary problem, then, is oe determination of f. 


‘ The corresponding logarithmic decrement is —7z{, while the rate of growth of 
n per unit time is }fkc. 
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‘The approximations (2.9) and (2.10) correspond to what we may designate 
as weak coupling of the two fluids. In the Kelvin-Helmholtz model, on 
the other hand, we have strong coupling, and at the transition from stability 
to instability [U = U,,, s < 1, see Lamb 1945, § 268 (4)] 

s( Unite? mae i, (2. I 1) 


where c,, (which depends on surface tension, as well as gravity) is the 
minimum value of c,. We infer from this that the approximation (2.9) 
may be regarded as valid only if 


ja + 2B\(Oy/Un)? < (¢/Cm)?. (2.12) 


‘The fact that waves do begin to form at wind speeds very much less than the 
Kelvin-Helmholtz value U,, implies that the inequality (2.12) may be 
assumed in describing incipient wave formation. Of course, (2.6) could 
be solved without further approximation, but the boundary value problem 
to be solved (§3) would be far more difficult in consequence of the fact 
that the dependence of c on k therein could not be specified a priori. 


3. EQUATIONS OF MOTION FOR THE AIR 
The equations of motion governing a small perturbation of a two- 
dimensional shear flow U(y) in an inviscid incompressible fluid of density 
Pq are (see Appendix) 


Pa(u,+ Uu,+vU,) = —p,, (3.1 a) 
Pa(%+ Uv,) ae (3.1 b) 
u,+v, = 0, (3.1 c) 


where u and wv denote the x- and y-components of the perturbation velocity, 
p the perturbation pressure, and subscripts partial differentiation. Intro- 
ducing a stream function according to 


u= —y,, v=4,, (3.2, b) 
and assuming (by virtue of linearity) ¢ and p to exhibit the same dependence 
on x and ¢ as y in (2.2), we obtain 

pal(U— chy — U,#) =p, (3.3a) 
Pa k?(U —c)b = py. (3.3 b) 

Elimination of p then yields 
(U—c)y,,,, —[k(U —c) + Uy] = 0. (3.4) 


Equation (3.4), which was studied originally by Rayleigh (1880; also 1945, 
§ 366(7)), is the inviscid form of the well known Orr-Sommerfeld equation 
(Lin 1955, (1.3.15)). We observe that the omission of the viscous terms 
leads to a (regular) singularity at U =. 

‘The subsequent analysis will be carried out in terms of the dimensionless 


variables €, w, and ¢, defined by 


€ = hy, U-—c= U,w(6), wb = U, d(E)7(x,t), (3.5, b,c) 
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where U, is the arbitrary reference velocity introduced in (2.3). Introducing 
these in (3.4), we obtain 
$d’ —[1+(w"/w)]d = 0. (3.6) 
The boundary conditions to be imposed on ¢ are dictated by the require- 
ments that the interface (originally y = y») shall remain a streamline and 
that the disturbance shall die out at infinity. Noting that the horizontal 
velocity relative to the wave is approximately U—c, we obtain for the first 
condition 
Y,/((U—c) =ikn at y=Vot+n = Vo (3.7) 
We remark that the approximation u —c = U—c assumes only that the slope 
of the wave is small (amplitude small compared with wavelength, that is, 
ka <1). If y were interpreted as the distance above a fixed plane it also 
would be necessary (in order to linearize the boundary condition) to assume 
the amplitude to be small compared with the length of any region in the 
neighbouring flow over which local variations must be considered. This 
would be a rather severe restriction in consequence of the large velocity 
gradients near the boundary, but it may be avoided by assuming y to specify 
the streamline that would have been a distance y above the undisturbed 
interface and applying the boundary condition at y= yy. Perhaps the 
most satisfactory proof of this statement is to replace the independent 
variables x and y in (3.1 a, b,c) by x and ‘’, where '¥’ is the total (as opposed 
to perturbation) stream function (von Mises transformation), write U = U(¥) 
rather than U(y), linearize with respect to the velocity perturbations about 
U = U(¥), and apply the boundary condition (3.7) at ‘Y = ‘Yo, where 
U, = U(¥>). The form of the equations is simpler in terms of y, however, 
and we shall retain it with the implicit understanding that it specifies a given 
streamline, rather than distance from a fixed plane. We also note that the 
boundary condition may be applied at any streamline, y = yg, for which 
ky) < 1; thus, U, need not vanish, and it is not necessary to prescribe the 
shear profile all the way to the surface. In effect, we assume that a small 
transition layer, within which 0 < U < U), moves with the surface wave*. 
Introducing the dimensionless variables of (3.5), the condition (3.7) 
becomes 
de = We, (3.8 a) 
where the zero subscripts imply € = &). ‘The requirement that the dis- 
turbance die out at infinity yields simply 
d>0 asE> om. (3.8 b) 


The perturbation pressure is given by (3.3 a) or, in terms of the dimension- 
less variables, by 


p = p,U; R(wd' —w'd)n. (3.9) 


* The assumptions kyy < 1 and U-c > v alone, without further linearization 
or the neglect ef the dominant viscous forces, imply that p and v may be assumed 
constant across the transition layer within the approximations of Prandtl’s boundary- 
layer theory, as applied to the disturbed flow. 
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Comparing (3.9) with (2.3) at € = &j me invoking (3.8 a) we obtain 

4+iB = Wo(d)—w)- (3.10) 
This completes the formulation of the sien boundary value problem, 
which now is stated by (3.6), (3.8 a, b), and (3.10). 


4. IMPLICIT SOLUTION 
We next obtain an implicit integral expression for «+78, both in order 
to illuminate the decisive role of the singularity at U = c (w = 0) and as the 
first step in an approximate determination of 8. Multiplying both sides 
of (3.6) by ¢*, the complex conjugate of ¢, integrating from € = &) to € = o, 
integrating #’¢* by parts, and imposing the boundary conditions (3.8 a, b), 
we obtain 


| (d'2+[1 + (w"/ee)] IB} dé = [6819 = —ao di (4.1) 


Now, from (3.10), since w is approximately real, f is given by the imaginary 
part of w,» ¢,, and the only contribution of the integral in (4.1) to this imaginary 
part must come from the singularity at w = 0. It follows that 
‘a ) 
B= —F. | |bi(e0"/w) dé}. (4.2) 
J & J 
‘The path of integration in (4.2) must be indented either over or under 
the singularity at € = €, where w(€.) = 0, and on this choice depends the 
sign of 8. ‘This rather delicate and evidently decisive question has been 
discussed in some detail by Lin (1955, § 4.3 and ch. 8), and we merely state 
the conclusion that the path must be indented under the singularity* 
Applying the calculus of residues then yields 


§ = —7\d, 2(z,. it w, = V, (4.. 


¢ 


s*) 
— 


a 


where the subscript ¢ implies evaluation at € = 
The result (4.3) implies that, 7 the absence of Micipatie e forces, a motion 

of the type (2.2) will be stable or unstable according as the curvature of the 

wind profile (U”) at that elevation where the wind speed is equal to the 

wave speed is positive or negative, respectively. We infer from this result 

that: 
(1) only those waves having speeds in that range of the wind profile for 
which — U” is large may be expected to grow; the lower limit for c 
may be imposed by the existence of a sub-layer of linear profile or 
by the interaction of the waves with the wind profile, while the 
upper limit will be rather less than the wind speed outside of the 
boundary layer ; 

(2) in the initial phases of wave formation, those waves having speeds 
well down in the profile (large — U”) may be expected to predominate ; 


‘This assumes ¢ >,»is real. In the next approximation .f {c} 0, so that the 
singularity lies slightly aon the real axis (assuming w/z, < 0), and the path of 
integration in (4.2) passes under the singularity without the necessity of indentation. 








| 
| 
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(3) experimental measurements of aerodynamic forces on stationary 
wave models (Stanton et alia 1932; Motzfield 1937; Thijsse 1951) 
may not yield significant values of such parameters as Jeffreys’ 
sheltering coefficient, since the point at which U = ¢ then occurs 
right at the boundary. 

We remark that the result (4.3) is directly related to ‘Taylor’s (1915) 
theorem that the momentum transfer from mean to disturbed flow in a 
shear profile is proportional to —U"; alternatively, the result may be 
interpreted in terms of the vortex theory of instability developed by Lin 
(1955, § 4.4), following earlier ideas of von Karman (see also Rayleigh 1945, 
vol. 2, p- 39] i 

We now proceed to develop an integral expression for ¢, which may be 
used for an approximate determination of 6 in conjunction with (4.3). 

3 


We start from (3.6) in the modified form 


(wd' —w'd) = we. (4.4) 

Integrating both sides of this equation between € = €, and & ©, imposing 
the boundary condition (3.8 b) at € ©, and noting that 

lim wd’ =0 (4.5) 


by virtue of the fact that 4’ can have only a logarithmic singularity at the 
singular point of the differential equation*, we obtain 


] . 
d, ; | wd dé. (4.6) 


" 
ie 


‘The principal attribute of (4.6), aside from the smoothing effect intrinsically 
associated with integration, is that it is locally insensitive to errors in (the 
approximation to) ¢ in the neighbourhood of the singular point w = 0. 
We conclude this section by remarking that the approximations 
introduced up to this point may be regarded as inherent in our model. 
Additional approximations will be made in the following sections in order 
to simplify the mathematical analysis, but these are, at least in principle, 
subject to improvement without modification of the model. 
5. APPROXIMATE SOLUTION FOR f 

Che exact solution of the differential equation (3.6) in terms of known 
functions for typical velocity profiles does not appear to be poss'ble. ‘The 
usual (Heisenberg) expansion in the parameter k? for the problem of thin 
(in wavelengths) boundary layers (Lin 1955, §5.5) is not well suited to the 
relatively thick turbulent boundary-layer profiles that are of interest in the 
present problem. 


It may he proved that 


b = dell +(w? wi(E—E,)log(E—E,) + OF —E,)] 


pal 


in the neighbourhood of & €, (cf. Rayleigh 1945, § 369; Lin 1955, (8.1.8)). 


F.M. 
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We shall attempt here only an approximate determination of 8 through 
(4.6), using the approximation 


$= af. (5.1) 
It is evident that (5.1) satisfies the boundary conditions (3.8 a, b), while the 
differential equation (3.6) is satisfied if either |w”) > |w (small €) or 
w’ < w| (large €). ‘These conditions are not, of course, sufficient to 


guarantee that the approximation (5.1) is uniformly valid for all €; in 
particular, it fails completely at the singular point €.. We anticipate, 
nevertheless, that, in conjunction with (4.3) and (4.6), it will lead to an 
approximation to the magnitude of ¢, which is adequate for at least the 
qualitative description of the spectral distribution of the parameter f. 
We shall be content with this goal, remarking that the difficulties associated 
with the accurate determination of wind profiles scarcely justify anything 
more elaborate at this time. 
Substituting (5.1) in (4.6) and that, in turn, in (4.3), we obtain, 


We shall evaluate the integral in (5.2) for a logarithmic profile, which, 
for turbulent flow over water, has the support of both theory (see Coles 
(1956) for a recent survey) and experiment (Roll 1948; see also Hay 1955; 
Ellison 1956). ‘The form that most directly fits the description of §3 
above is 


wn 


U(y) = U,+ U, log(y/yo), (5.3 a) 
where U, is the lower bound of the logarithmic profile and U, the reference 
velocity, which is now defined as the coefficient of the natural logarithm in 
the velocity profile. ‘The majority of experimental results are presented 
in the form (except that common logarithms are apt to be used) 


U(y) = (Ux /«)log(v/z9), (5.3 b) 


where Ux = (79/p,)"? is Prandtl’s shearing stress velocity, t) the shear stress 
at the surface, « Karman’s universal turbulence constant, and 2» an effective 
roughness parameter, in the notation of Roll (1948); the definitions of U, 
and « are those commonly adopted, but z) has been used in somewhat 
different senses by different writers. We remark that the reference pressure 
p, UF now is proportional to the shearing stress 7. 

The distribution of (5.3) evidently holds only for a limited range of U. 
We designate the wind speed at the outer edge of this range as U,, ; for our 
purposes, U 10U, = 25U, is a rough but adequate approximation for 
the wind speed at two metres above the water*. ‘The determination of 
the lower limit of validity, Uj, is both more important and more uncertain, 
reflecting the fact that, in reality, there is not a sharp transition from 


*'The reference position of two metres is that adopted by Roll; at 10 metres 
U 30U, might be a better estimate. 
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logarithmic profile to sub-layer (laminar or otherwise). Estimates of both 
U, and z, will be discussed in $7. 
Substituting (5.3.4) or (5.3 b) in (3.5 b) and (5.2) and noting that, by 
definition, 
c= l ‘+ U, log( Ve Vo) sie U, log(y, Sy); (5.4, b) 


we obtain w = log(é/E.), (5.5) 


B= al | e~* log?(E/E..) ae | (5.6a) 
= réil | ee" log? udu | (5.6 b) 
¥ J 


The integral in (5.6b) may be evaluated for small or moderate values 
of the argument €, by writing the range of integration (1, 0) as the 
difference (0, 0)-(0,1), identifying the infinite integral as the Laplace 
transform of log?u [Erdélyi 1954, §4.6(13)], and evaluating the finite 
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Figure 1. The function B (ky,), as defined by (5.6) and (5.7) 
integral as a power series, following the expansion of the exponential. 
An asymptotic development for large €, may be obtained by the usual device 
of repeated integration by parts or by expanding log” in powers of (u—1). 
The end results of these operations are 


P ; P i= 'er i? — 
B(E.) = wé,4 7? + log? (yé.)+2 > at a (5.7 a) 
( 
BE.) ~ 4€e-2°*(1 — 6€71 + 31€5*...), 5.7 b) 


where log y = 0:5772 is Euler’s constant. We find that 8 has a maximum 
value of 10-2 at €, = 0-017 and that (5.7 a) is adequate for numerical calcul- 
ation over the range of practical interest (€,< 1). ‘The asymptotic series 
(5.7 b), on the other hand, is accurate only for very large €, and is of interest 
primarily ag an indication of the exceedingly sharp cut-off of the damping 
ratio for very large values of ky,. 

The numerical values of £ are plotted in figure 1. 
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6. DAMPING RATIO vs WAVE SPEED 
We now proceed to express f and @, as functions of the dimensionless 
wave speed c/U,. ‘The corresponding wave number for gravity waves 
may be determined from (2.8 a), whence 


ke = p/e*. (6.1) 
Solving (5.4b) for y, and multiplying by (6.1) then yields 
€, = hy, = 1(U,/ep~e™, (6.2) 
where we have introduced the dimensionless, wind-profile parameter 
OQ <unsth (6.3) 


"he function ky,/Q is plotted in figure 2. The corresponding values of p> 
as calculated from (5.7 a), and ¢,/s, as given by (2.10), namely, 
C,/s = (U,/c)PB[Q(U,/c)e!™}, (6.4) 


are plotted in figures 3 and 4 for values of Q that appear to be typical (cf. § 7). 
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Figure 3. The parameter f, as determined 
by (5.7a) and (6.2); cf. figures 
1 and 2. 











c/U, 
Figure 2. The variation of ky, with wave 
speed ; see (6.2) and (6.3) 
The damping ratio associated with viscous action in the water, based on 
the assumptions stated in $1 (free surface with laminar flow), is given by 
(Lamb 1945, § 348) 


C,, = —4v,, Ric. (6.5) 
Eliminating k through (6.1) and superimposing ¢, and ¢,,, we obtain 
C = BsU* oc? —4y,, gc. (6.6) 


Equating this expression to zero yields 


U, = (4, g s)! 3l(¢, l 1) P] Ais (6.7 a) 
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for the wind speed at which the energy available from the air stream is just 
sufficient to maintain a small amplitude gravity wave of speed ¢ against 
laminar dissipation in the water. If we take the c.g.s. numerical values 
v, = 10-*, g = 980, and s = 1:2 x 10-%, (6.7) reduces to 





















































U, = 32[(c/U,)B]-* cm/sec for ¢ = 0. (6.7 b) 
This result is plotted in figure 5 using the values of f given in figure 3. 
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Figure 4. ¢,s, as given by (6.4). 
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c/U, 
Figure 5. The minimum reference speed U, for the initiation of gravity waves, as 
given by (6.7 b); the wind speed at two metres is roughly 10U,. 


It appears from the results of figure 5 that the minimum wind speed 
(at two metres) for the initiation of gravity waves is of the order of magnitude 
of 100 cm sec and might be as small as 80 cm/sec; the corresponding wave 
speeds would lie somewhere between 40-50. cm/sec. ‘These figures are close 
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to the respective measurements of Jeffreys (1924) and Scott Russell (1844), 
but it is important not to lose sight of the approximations on which they 
are based. The equivalent value of Jeffreys’ sheltering coefficient, as 
given by 


t(Ua—e), (6.8) 


7. 'THE WIND PROFILE PARAMETERS 
‘The available measurements on wind profiles over water (see Roll 1948; 
Neumann 1948; Hay 1955; also other papers cited by Ursell 1956 and 
Ellison 1956) all tend to confirm the logarithmic law posed in (5.3), but 
they provide little information on the laminar sub-layer and conflicting data 
on the effective roughness parameter. We consider here three possibilities. 


(a) . lerodynamically smooth flow 

If the departures from a plane surface do not penetrate the laminar 
sub-layer and the neglect of surface currents is justified the velocity profile 
should resemble that for flow near a smooth wall. Nikuradse’s data yields 
(Prandtl 1952, p. 126) 


U(y)/Us = 5-75 logyo( UVey/v,) +550, (7.1) 
corresponding to x = 0-40 and 2) = v,/9-06U,% in (5.3b). We then have 

%o 1-7 x 10°-*U;', Q = 2-7U,3, (7.24, b) 
in c.g.s. units (g = 980, v, = 0-154). 


Roll has made measurements for wind blowing over smooth ponds 
which yield values of 3) which are not only appreciably smaller than those 


predicted by (7.2 a) for Ux 10cm/sec but also decrease with Uy. The 
scatter in these data for small Uy is considerable, but the trend is quite 
definite. ‘lhe most plausible explanation would appear to be in the exis- 


tence of a surface current in the water, but further investigation of this point 
is desirable*. 

If the logarithmic profile of (7.1) is matched in ordinate and slope to 
a linear profile it is found that Uy = 7-8U% at Ux Vvo/v, = 1/« (cf. Roll 1948). 
\n examination of pipe data (Prandtl 1952, p. 128) reveals that the logarith- 
mic and linear profiles are established only for U/U, >15 and U/U, < 5, 
respectively, but the logarithmic profile does not appreciably overestimate 

U"/U’ as a function of U for l SU,. 


(b) Moderate wind over wavy water 

Roll found that for values of Uy between about 10 and 30 cm/sec and 
heights up to two metres above a disturbed surface his measurements 
(which exhibited a much smaller scatter than those for Uy, 10 cm/sec) 


* Prandtl (1952, p. 130), referring to this data, states that “‘ for values of U 


below about 10 cm/sec the flow in the lowest 2 metres seems to be laminar’’. ‘This 


* 


interpretation is difficult to accept, since his own “ universal velocity function ”’ 


thid, p. 128) implies a value of y from 1-3 mm for transition at U, 10 cm/sec. 
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could be fitted to (7.1) when y was measured from the wave crests. In terms 
of the distance from the mean surface, he found 


eer 5-75 log, o[Us(y + @)/v,,] + 1:35, (7.3) 
where a was the amplitude of the surface waves. ‘This implies that 
Zo 3x 10°E-*, Q 11-6U-%. (7.4 a, b) 


It should be recalled that we have defined y as a streamline parameter, so 
that the velocity profile of (7.3) is not directly significant and serves only to 
illustrate the magnitude of possible deviations from (7.1). We shall see 
(figure 6) that the resulting differences in @, are negligible. 

The results (7.3) and (7.4) are characteristic of aerodynamically smooth 
flow, so that it appears reasonable again to assume that (7.3) is valid down to 
about U, = 8U,. 


(c) Fully dew lope d rough jiow 

It has been argued (cf. Ursell’s (1956) and Ellison’s (1956) remarks) 
from dimensional considerations that, since (i) the apparent roughness 
of the water surface depends on the waves produced by the wind and (ii) the 
roughness parameter in pipe flows is found to be independent of viscosity 


and to depend only on the nature of the roughness for U, 2/v, > 3 (rough 
flow), zy) should be proportional to U{/g for sufficiently rough flow. We 


remark here that this conclusion appears to assume that capillary ripples 
are unimportant in determining the apparent roughness, although this 
would appear doubtful; nevertheless, Hay (1955) has made measurements 
for values of U, between about 25 and 55cm/sec that tend to confirm this 
hypothesis and lead to the approximate results (which easily might be off 
by as much as 30°,,; cf. Ellisons’ remarks) 

% = 8x 10°U3, QO = 125% 10". (7.5 a, b) 
The corresponding value of U,is probably best estimated as thataty, = 3025, 
namely U, = 8:5U,. 

The discrepancy between (7.4) and (7.5) emphasizes the order of 
magnitude of the uncertainties with which we deal; for example, at 
U, = 30cm/sec, where the data of Roll and Hay overlap, the values of 29 
obtained from(7.4 a)and(7.5 a)are 2 x 10-*cm and 7 x 10-* cm, respectively. 
In any event, our model cannot be expected to have more than qualitative 
significance for rough flow. 

The foregoing results are illustrated in figures 6 and 7, where we have 
plotted ¢, vs A (= 2z/k) for the combinations of data indicated in table 1 
(the values of © used in the actual calculations were rounded off to the 
parenthetic values in order to reduce the numerical work; the rounding 
errors are smaller than the other uncertainties). "Uhe differences in 3, for 
given Uy are seen to be relatively unimportant. We again emphasize 
that the results for the higher wind speeds must be regarded as especially 
dubious in consequence of the linearization. The value of —{,, given by 


>w 


(6.5) also is plotted as a reference base, but the actual damping from the 
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water almost certainly would be very much larger for values of Uy as high 


as 30cm/sec. The values of c/U, = 2, 3, 4, and 5, corresponding to various, 
possible lower limits Uy = c are marked off; on the basis of the preceding 


7 | | : 
l ¥ } 0 | l 200 


| 
| Item | ; 2 Eqs 

| (cm/sec) | (cm) | (cm/sec) | | | 
| a; | 5 | 34x00? | 100 | 22x10 | @.1),07.2) | 

| | | Ox | 

ae 10 | L7 X10 | 290 | 2-:7x10 | GTDAG2y | 
| | | | (3x10-3) | 
| be 10 | 73x10 250 | 1:2x10-2 | (7.3), (7.4) | 
| | | | (10-*) | 
| by | 30 | 24x10 g50 =| 4:3x10-8 | (7.3), (7.4) | 
| (3 x 10-*) | 
| te 7 30 7210-2 | 600 | 1:2x10-2 (7.5) | 
| | | (10-2) | 


Table 1. The parenthetic values for {2 were used in calculating figure 6. 
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Figure 6. Damping vs wavelength for the data of table 1. 


remarks, the logarithmic profile probably does not overestimate ¢, fot 


{ c > 4U,( = 10U,) but certainly would do so for values of c as low as 2U, 


fal’. 
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Figure 7. Linear plot of cases a; and dj; cf. table 1 and figure 6. 


8. MEAN TANGENTIAL WAVE DRAG 
A rough check on the mean value of 6 (the minimum wind speed estimates 
of § 6 provide a check only on values of f in the neighbourhood of its maxi- 
mum) may be obtained by calculating the (space or time) mean tangential 
drag on the waves, namely, 


Tt, = KR p,(ikn)*}. (8.1) 
Munk (1955), using data from Van Dorn’s (1953) experiments, has obtained 
the expression 

et, = ¢ plore? U.. (8.2) 

where c’ = 0-68 x 10-6 if U,, is measured at 10 metres above the surface. 

Let S(k, #) be the power spectral density of the surface* in the notation 

of Eckart (1953 a), @ being the angle between wavefront normal and wind 

direction (so that the effective component of the wind is U cos @); the power 

spectral density of the wave slope then is k?S(k, 4), and substituting p 
in (8.1) from (2.3) yields 


i pal ta | B cos?@k2,Sk dkdé. (8.3) 
The spectrum for a fully developed sea (Neumann, as cited by Munk 
(1955)) may be approximated by 

R2S(R,O)k dkdO = 4(27)-1 20 10% 2" &"d(c/U,) dO 0| < 34%), 


0 (6 4,), 


(8.4) 


2 


where o? is the mean square slope, which may be expressed in the form 
o = Wgv,y FU. (3.5) 


The surface displacement now is assumed to be composed of a distribution of 


sine waves having random phases. 
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with 6 = 1-1 x 10-4 (Cox & Munk 1954). Munk (1955) also introduces a 
sheltering coefficient f,, (s in his notation), defined such that 

c’ = By 5f(%), (3.6) 
where, for our model (in which U cos 6 would appear as the effective wind 
speed), f(4,) is simply the mean value of cos?6 over (— 54, 54). Substituting 
(8.4) and (8.5) in (8.3) and comparing the result with (8.2) and (8.6), we 


tind 
, 4 l b's s -U ~/U, een Cc ae 
Py ve 5 (; ) | y Be~ 20 a a(z-). (3.7) 


Using the values of 8 given by figure 3 with Q = 10-? (the value most 
likely to be representative of a fully developed sea), assuming 8 = 0 for 
c < 3U,, andchoosing U 12U,, (at 10 metres above the surface), we obtain 
>4, = 20x 10°", which compares with the value 10? inferred by Munk 
from observed rates of wave growth. ‘The corresponding value of c’, 
using 6 = 1-1 x 10-4, is 2-1 x 10-® cos?@; with Munk’s value of 6), namely 
130°, we have c’=1-5 x 10-8, which compares with the experimental value 
of 0-68 x 10°-&. We emphasize, however, that our estimate of 8,, might 
be otf by a factor of two either way in consequence of our estimates of Q 
and ¢p. 

We conclude that the calculated value of c’ is in agreement with the 
experimental value in at least order of magnitude, from which we infer 
that the mean value of 8, weighted as in (8.7), also is of the right order of 


magnitude. 


9. CONCLUSIONS 

We conclude that the model of an inviscid shear flow in air over water 
having zero mean motion provides an explanation of the generation of 
gravity waves which is in qualitative agreement with observations of mini- 
mum wind speed and mean tangential drag. Quantitative comparisons 
would require, on the theoretical side, a more accurate solution of the 
boundary value problem, and, on the experimental side, much more precise 
data on wind profiles. We also emphasize that the model is intrinsically 
limited by the neglect of turbulent interaction between surface wave and 
wind profile and the linearization of the equations of motion. Finally, we 
note that the results ot § 2-5, regarding energy transfer from shear flows 
to surface waves, are applicable to such problems as the flutter of membranes 


ind panels. 


! am grateful to Lester Lees of the California Institute of ‘Technology, 
Nicholas Rott of Cornell University, Cark Eckart, Walter Munk and Charles 
Cox of the Scripps Institute of Oceanography, and my colleague Kurt 


Forster, for helpful discussion and criticism during the course of the above 


research. 
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APPENDIX 
The Euler and continuity equations for a viscous incompressible fluid 
can be written 


A 


CU; Ou; Op;; 
yO 1 Shy (Ala 
( t ( Xv; p ¢ Xj 

Ou 

~ = (), (Alb) 


where x; denotes a Cartesian coordinate, u; a velocity component, p;; a 
stress tensor component, p the fluid density, and the usual summation 
convention is implied. We resolve the dependent variables according to 

u;= U;+u;+u;, bij = Pip +i, +Diis (A2a, b) 
where U,+u; and P,;4 p;; represent a solution to (Ala, b) having the two- 
dimensional (x, and x.) mean values (with respect to either x, or t) U; 
and P;,; plus turbulent fluctuations uw; and p;,, and u; and p;; represent a small 
perturbation with respect to this solution. Substituting (A2a, b) in (Ala), 
neglecting terms of second order in the perturbation flow, and invoking 
the requirement that the unperturbed flow satisfy (Ala), we obtain 


) 





Cu Cu; »(U,+u: 1 cp’ 
= +(U;+u,)— 4 Pa sk Pi (A3 a) 
ct ‘ OX; OX; p OX; 

Cu; 


Taking mean values with respect to x3, we may place the results in the form 








CU: : our — ol f low=;w; “Wk 

> “13, le dy, Pi Mud (A4a) 
Our 
ot SO | (A4b) 
OX 


where we have introduced the perturbation Reynolds stress (invoking the 
equations of continuity for both u; and u; in its derivation) 


7 


Vr; = p(u; u; +u-u;) (A5 a) 


a) 





= plu;(u; —u’)+u‘(u; —)u‘], (A5b) 
with (A5 b) following from (A5a) by virtue of u; = 0. 

Equations (A4a,b) are the equations of motion governing the two- 
dimensional perturbation of a two-dimensional shear flow if teo-dimensional 
is understood to imply mean value with respect to the transverse coordinate. 
They differ from the equations of motion for a non-turbulent perturbation 
flow only in the presence of the dl which represent the interaction between 
the fluctuations in the original and perturbation flows and, as might have 
been anticipated, are the first order perturbations of the usual Reynolds 





Stresses pu. ue 
Equations (3.1a,b,c) follow from (A4a,b) if we take x, = 4, x, = y, 
U, = U(y), U. = U, = 0, u; = 4, , = V, p;; 6,;,p, and r;; = 0. 
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NOTE ADDED AT PROOF STAGE 


A stochastic analysis recently presented by Phillips (1957) develops a 
model which appears more realistic than that proposed by Eckart (1936 b) 
and arrives at results in reasonable agreement with experiment. It seems 
that the mechanism proposed by Phillips is complementary to that 
presented here, but further investigation is required to establish the 


appropriate connections. 
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SUMMARY 

The turbulence produced by a multiplicity of small air jets 
has been investigated and comparisons are made with other 
turbulent flows. ‘The eddy Reynolds number is large (R; = 250). 

The energy spectrum was measured, as well as the skewness 
and kurtosis of du/dx. The effect of the finite leneth of the hot 
wire is considered and the corrected results indicate that the 
spectrum follows the law of Kolmogoroff in an intermediate range 
and appears to fall off with the (—6)-power of the frequency in 
the viscous range. ‘his range is limited at the upper end of the 
frequency range by electrical noise. Special precautions reduced 
this noise to the level of thermal agitation in the hot wire. 

The ultimate limit to spectral analysis is imposed by the 
molecular agitation of the gas. ‘This limit is evaluated and 
compared with the spectrum of turbulence. It appears that the 
spectrum of o%u/dx? merges with the spectrum of molecular 
agitation without a distinctive separation. 


1. INTRODUCTION 

The fine structure of turbulent flows has been the object of numerous 
investigations but complete understanding has not yet been reached. In 
this paper, we describe certain recent results that may contribute to the 
clarification of the situation. 

The theory of turbulence becomes greatly simplified if one assumes 
incompressibility, homogeneity and isotropy of the velocity field. Under 
such conditions, one can define a one-dimensional spectrum F(k) which 
is a function of a wave-number k. Let us define k,, as the wave-number 
such that half the kinetic energy of the motion is contributed by the region 
k > k,, and let us define k, as the wave-number such that half of the viscous 
energy dissipation occurs in the region k > k,,. 

Two important theoretical predictions about F(k) are as follows 
(Batchelor 1953, ch. 6). According to Kolmogorofl, the spectrum should 
fall as R~*? in some subrange of wave-numbers in the range k, < k < k,, 
provided that k, <k,. According to Heisenberg, the spectrum should 
fall as k~“ in the region k > k,. ‘These two statements can be formulated 
in a single expression, and the theory shows that F(k) should cross a limit 
given by cv?k in the vicinity of k = k,, where c is a universal constant and v 
the kinematic viscosity. Physically, this means that spectral components 
such that F < cv*k are strongly influenced by viscosity. 
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The theory of Heisenberg implies that the third derivative of the velocity 
does not exist, or, more exactly, that the mean-square of this quantity, 
given by the integral of k° F(R), does not converge. 

From Kolmogoroft’s theory it can also be deduced that the parameter S$ 
(skewness of du/ex) should be a constant, provided hk, < k,,. 

‘The experimental studies of the fine structure of turbulence have been 
carried out in a single gas (air), virtually at a single temperature and pressure. 
The velocity fluctuations are sensed by hot-wire anemometers and the effect 
of the finite length of the wire is often disregarded. ‘The hot-wire anemo- 
meter gives an electric signal proportional to the velocity fluctuation plus 
a certain amount of electronic noise, so that some kind of filter is desirable 
to reject unnecessary noise. . 

In general, the turbulent fluid passes by the hot-wire with a mean velocity. 
U which is much larger than the velocity fluctuation and, with the hypothesis 
made by ‘Taylor, the time derivative of the signal is taken as a fair approxi- 
mation to the space derivative. ‘This allows the measurement of quantities 
such as: 

== (y*)'2, A = u' ((du/dx)*)-??, 
where < ) denotes a time average, wu is the velocity fluctuation in the 
direction of the x-axis and A is the dissipation length parameter. Ifk, < k,, 
it can be shown that A = k,,k-, and the eddy Reynolds number 


R,=w'dalv 
becomes proportional to (k,/k,)?*. The verification of the theoretical laws 
thus requires a large value of R,, say R, > 100. 


Grid flows 

A large variety of investigations have been carried out in the turbulent 
wake of a plane grid obstructing a parallel flow. ‘The requirements of 
homogeneity and isotropy are reasonably well satisfied here and for low 
velocities incompressibility is a reasonable approximation. R, usually 
varies from 20 to 100. In practice, it becomes difficult to observe the fine 
details of the motion when R, > 60. This is due to the fact that w’/U < 0-02, 
resulting in a weak electric signal from the hot-wire. At some frequency w, 
the spectrum of the signal merges into the noise of the instruments, 
consequently the smallest observable details of the motion are comparable 
to the length /, = U/27w,. With a given grid size and a given hot-wire, 
R, increases as U1? whereas /y increases at least as U. ‘Therefore, as R; 
increases, the resolution of the fine structure deteriorates. ‘To avoid this 
effect, one should use lower viscosities or larger wind tunnels. Because 
of these limitations, it was not possible either to really support or to disprove 


the theoretical spectral laws with the work on grid flow. 

It has been found the it S varies very little and that quantities such as 
{(du/ox)?) and <(0*u/ox?)?) are related to each other in the way required 
by the laa wl Pie Measurements of the quantity ¢(0%u/dx3)? ) 
have been reported by Batchelor & ‘Townsend (1949), ‘Townsend (1951), 
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and Stewart & ‘Townsend (1951); Liepmann, Laufer & Liepmann (1951) 
carried out similar investigations of the third derivative but report more 
cautiously that ‘‘it appears to be finite’. 

‘lownsend (1951) proposed a simple model for the fine structure of 
turbulence based on stretching of vortex sheets or vortex lines. It leads to 
a spectrum falling as an exponential function of k? when & is large. 
Measurements in grid flows in the range 15 < R, < 30are in good agreement 
with the model, but discrepancies seem to appear as R, increases. 


Shear flows 

The structure of turbulent shear flows has been studied by Laufer 
(1950, 1953), Klebanoff & Diehl (1951), Corrsin & Uberoi (1951) and others. 
The departure from isotropy and homogeneity are serious but it affects 
principally the spectral components with wave-numbers near k,. Ina flow 
with R, = 300, Laufer (1950) found that the fine structure is nearly isotropic, 
that a portion of the spectrum falls as k-®* and that the high frequency end 
of the spectrum falls as k~’. (He obtained this without applying a correction 
for the finite length effect of the hot-wire.) ‘This last result implies that 
((0%u/ex*)?)) is infinite or at any rate that it is determined by values of k 
too large to be measured. 

In an attempt to clarify this situation, we designed a simple device to 
produce turbulence with large R; and small /,. The object of this paper 
is to describe the procedure and to comment on the results. 


2. EXPERIMENTAL INSTALLATION 
The * Porcupine’ 
A plywood box 60 x 60 x 120 cm has five perforated sides, as shown 
on figure 1. The sixth side is faired into a duct leading to a 20 h.p. blower. 


Figure 1. The ‘ Porcupine’. The mixing of 80 small jets 
produces a strong turbulence in the region marked A, B, C. 


Each of the 80 openings is fitted with a round tube, 10 cm long and 2:5 cm 
inside diameter. In operation the air enters the box through these tubes, 
forming 80 jets which merge into a turbulent flow exhausted into the duct 
by the blower. ‘The Reynolds number based on velocity and diameter 
of a single jet is about 35000. ‘lhe mean velocity of the flow in the 
parallel section is 450 cm/sec and the turbulence levels measured at the 
points A, B, C on the box axis are respectively 30°,, 12% and 5%. 
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At point A, the turbulence is so strong that the hot-wire signal is not 
proportional to the velocity fluctuation, because of the intrinsic non-linearity 
of the instrument. At points B and C this distortion becomes less serious. 
The resolution length /y is larger at B than at C by about 60%. and, as a 
compromise between the distortion and the loss of resolution, the author 
decided to perform all measurements at point B. 

No systematic variations of U and wv’ over the plane normal to the mean 
flow were found, at distances up to 15 cm from the box axis. Both quantities 
varied by 5 to 15°,, between one run and another and it appears that the 
merging of the jets produces different flow patterns, as observed elsewhere 
by Corrsin (1944). 


Hot-wire and preamplifiers 

Pt wires (diameter 1-25 microns, length 1 mm) heated up to twice the 
cold resistance were used. Whenever it became essential to minimize the 
electric perturbations (60 cycles ripple or contact noise) each wire was 
heated with a separate 22 volts battery and a fixed wire-wound series resistor ; 
all other metering devices were completely disconnected. 
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HOT WIRE TRANSFORMER CIRCUIT PREAMPLIFIER 


Figure 2. Input circuit for low electrical noise. 


‘l'wo low noise preamplifiers were used, each with a Weston 417 A triode. 
(he noise of each preamplifier is equivalent to that of a 600 ohms resistor 
applied to the input. With the amplifiers in parallel, a single channel is 
formed with a noise of 300 ohms. For measurements in the high frequency 
range, lower noise figures were achieved with the transformer input circuit 
of figure 2. The relatively small coupling condenser renders the output 
proportional to the time derivative of the signal, up to 20 kc/s. Above this 
frequency, the wire resistance and other impedances become significant 
and the overall response deteriorates. With a transformer ratio of 10 to 1, 
the input noise at 15 kc/s becomes equivalent to the thermal noise of a 
50 ohms resistor at room temperature. 

With an unheated wire of 60 ohms the total equivalent noise-generating 
resistance is 110 ohms. With the same wire heated to 120 ohms, that is, 
to approximately 600° K, we can expect a noise of about 300 ohms. ‘This 
means that the amplitude of the noise should increase by about 260” ,,, when 
the heating current is turned on. (‘This increase is noticeable on figure 5, 
at k 300, in the range where the electronic noise dominates the turbulent 


signal.) 
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We verified that the battery did not make an additional contribution 
to the noise. For this purpose we replaced the hot-wire with a 200 ohms 
wire-wound resistor, and observed that the passage of the heating current 
did not affect the noise of the wire-wound resistor. 


Thermal lag compensation 

The thermal inertia of the hot-wire appreciably attenuates the signal 
above 600 c/s. With the transformer input of figure 2, the coupling condenser 
provides the desired compensation from 600 c/s to 20 kc/s. It would be 
inadvisable to modify this circuit and obtain a constant gain between 3 and 
600 c/s. Indeed, the low frequency components of the signal would be 
so large that they would saturate the preamplifier. 

With a capacitive coupling into the preamplifier, it was necessary to 
limit the response of the preamplifier to the range of 200 c/s to 20 ke/s, 
in order to avoid this saturation by low frequency components. When the 
transformer was not in use, the thermal compensation was provided by an 
operational amplifier with an appropriate feed-back loop. However, the 
presence of a large turbulent signal made it difficult to adjust the time 
constant of the circuit by the usual square wave method. 

For measurements of w’ or of the low frequency portion of the spectrum, 
the preamplifier was by-passed completely. 


Differentiation and filtering 

An operational amplifier was used in the manner of analogue computing. 
The response of each differentiating unit increased with w up to 20 kc/s 
and decreased with w! from 40 ke/s upward. ‘This provided the appropriate 
filtering of the high frequency components of the electronic noise. 


Wave analysers 

Between 3 and 7500 c/s a General Radio wave analyser with a band-width 
proportional to the tuned frequency was used, and, from 100 c/s to 15 ke/s, 
a Hewlett Packard analyser built around a modulator and four amplifiers 
tuned for 20 kc/s. In order to avoid the formation of modulation images, 
it was necessary to insert ahead of this wave analyser a low-pass filter with 
cut-off at 15 ke/s. From 7 ke/s to 60 kc/s a Sierra analyser with constant 
band-width was used. 

The output of each analyser was squared, averaged, and reduced to 
a constant band-width, and the responses were matched at 100 c/s and 
10 kc/s. Good agreement was found in the two regions where analysers 
overlap. 

When the spectrum of the signal falls very rapidly, an analyser may be 
unable to separate the components properly. Indeed a tuned amplifier 
driven with a frequency lower than the resonance frequency has a gain 
proportional to the frequency. With four tuned amplifiers and a signal 
whose spectrum falls as w *, the output of the analyser is seriously in error. 
This difficulty can be avoided if the signal is differentiated once or twice 


F.M. O 








210 R. Betchov 


before reaching the analyser, since each differentiation multiplies the 


9 
spectrum by w*. 


3. ‘THE ENERGY SPECTRUM (WITHOUT CORRECTION FOR WIRE 


LENGTH) 
With a wave analyser tuned on frequency w the corresponding wave- 
number is k = 27w/U and the one-dimensional spectrum F(h) is equal to 


the mean-square of the analyser output multiplied by a constant such that 
(u2) = | F(R) dk. 
“OU 
In figure 3 we show F(A) as a function of k, with linear scales. A check with 
a variable low-pass filter indicated that half of the mean-square of w is 
contributed by components below k = 0-1cm™! and this confirms the 
existence of a maximum of F. It also suggests an energy scale of the order 
of 10cm, which is comparable with the size of the box. A check with 
two crossed wires indicated that w’/v’ = 1-3, where v’ refers to a velocity 
component normal to the mean flow. ‘This represents an appreciable 


deviation from isotropy. 
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Figure 3. Energy spectrum vs wave-number with linear scales. 


In figure 4 we show the spectrum with logarithmic scales. It follows 
the law of Kolmogorotf in the range 0-15 <k < 15. In the range 
70 <—k — 200 it appears to follow the law of Heisenberg. For k > 200 
the curve rises, but this is due to the electronic noise. ‘These results have 
yet to be corrected for wire length effect; this will be done in the next 
section. 

As previously discussed, it was suspected that the power law with 
which the spectrum falls off could be affected by the analyser itself. As 
a controlling experiment we measured the spectrum of 0?u/0t?, which 
is proportional to k'F. ‘Vhe results (see figure 5) show a corresponding 
portion of the curve falling as k-°, in agreement with the measurements 
of /(k) (for a wire of finite length). 

Figure 5 shows the spectrum of the noise, measured with a cold wire. 
At high frequencies it remains slightly below the signal obtained with the 


heated wire, as expected. % 
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Figure 4. Energy spectrum vs wave-number with logarithmic scales. 
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Figure 5. The spectrum 


of 02u/dx? (not corrected for wire-length effect). 
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In figure 6, the spectrum of figure 4 is compared with other results 
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Figure 6. Comparison of several spectra. 


We indicate by a dotted line the limit cv?k with the arbitrary choice 
c=5. Below this line, the spectrum is controlled by viscosity (see § 1). 
The grid flows nos. 1 and 2 (Liepmann, Laufer & Liepmann 1951; 
Klebanoff & Diehl 1951) appear to have energy scales too small for a 
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well-developed range of Kolmogoroff’s law. ‘The shear flows nos. 5, 6, 7 
(Laufer 1950, 1953 ; Klebanoff & Diehl 1951) and the author’s turbulent flow 
no. 4 have such ranges. ‘lhe curves 4, 5 and 6 seem to come together near 
k=20cm™. As for atmospheric turbulence, no. 3 (Crane & Chilton 
1956), the measurements have not been extended to the fine scale. 
(A pioneering contribution by Gédecke (1935) in the range 0-1 <k < 10 
unfortunately does not yield the spectrum.) It would be interesting and 
not too difficult to study the fine structure of atmospheric turbulence ; 
it might reveal a viscous cut-off somewhere beyond k = 10. 

Beyond the viscous limit, all the spectra fall rapidly enough to ensure 
the existence of ((du/dx)?) and <(0?u/Ox?)*). As for ¢((0%u/dx)*), it appears 
to depend upon spectral components beyond the range of present 
measurements. 


4. CORRECTION FOR WIRE LENGTH 

The finite length of the hot-wire introduces a particular type of error 
in the measurement of spectra. A complete theory of this effect has been 
given by Kovasznay & Uberoi (1953), with the assumption that all portions 
of the wire are equally sensitive. In the case of the present investigation, 
the ‘cold ends’ are about 0-05 mm wide each and the wire length is 1 mm, 
and, according to Betchov (1948), this means that a correction for finite 
length is necessary for k > 20. 

In figure 4 a dotted line shows the approximate values of the corrected 
spectrum. Furthermore, the theory of the wire-length effect shows that, 
if the measured spectrum falls as k~*, the true spectrum should fall as k~*. 
Since this correction has not been applied to the results of figure 5, the 
apparent seventh power law must be interpreted as evidence that the true 
spectrum falls off as k-®. ‘This means that the law proposed by Heisenberg 
has no support from these experiments, and that the results of figure 5 
suggest that the spectrum of the third derivative becomes flat for k > 70. 
There is little doubt that the correction for wire length applies also to 
Laufer’s measurements (curve no. 7, figure 6). 

As a subsidiary verification of the theory for wire-length effect, measure- 
ments F(k) were made with two wires 1 and 6 mmlong. The wave analyser 
was set on a particular frequency, and the two values of F were measured 
by simply switching from one wire to the other; this method minimizes 
the errors introduced by slow changes in the flow or the wire sensitivities. 
The same compensating network was used for both wires, since it was not 
possible to obtain an accurate adjustment of the compensating network 
in the presence of the large turbulence. The preamplifiers were operated 
with capacitive coupling (noise of 300 ohms), and therefore the measure- 
ments could not be extended beyond 15 kc/s. Figure 7 shows the measured 


values of », where 


= | F(R) |); 6mm 
TS | F(R) }! 1mm °* 
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‘Lhe theory of the wire-length effect predicts that 7 should start from 6 
at low frequencies and begin to drop at w = 250 c/s, finally levelling off to 
\6 at high frequencies. ‘lhe measurements support the theory. 
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Figure 7. Effect of wire length on the measurement of a spectrum. 
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5. SKEWNESS AND KURTOSIS OF Ou/dx 

Measurements of the mean-square of the square, cube and fourth 
powers of cu/ex were made, using three chains of 10 double triodes, properly 
biased. ‘Two chains were used for squaring and one for cubing, and the 
frequency responses were satisfactory. The largest errors were due to 
deviations from exact power law responses. 

It was found that S 0:-45+0-05 and y= 45+0-5; these results 
are not different from those obtained in grid flows with 30 < R,; < 60 
or in shear flows (Betchov 1956). 

The largest contribution to S originated from pulses of odu/dt lasting 
approximately 0-5 millisecond. ‘This corresponds to a length of about 
2mm and a correction for wire length may be necessary. 

For the dissipation length parameter, the result was A = 0-7 cm which 
corresponds to a wave-number 1-5 cm! and to R,; = 250. This result 
illustrates the fact that, at large values of R,, the length A is intermediate 
between the energy scale k,;! (about 10 cm) and the scale of dissipation k,! 
(about 0-025 cm). As is now generally appreciated, the name ‘ microscale’ 
which is sometimes given to A is not appropriate. 


6. MOLECULAR AGITATION IN A PERFECT GAS 

In these experiments, the spectrum could not be measured beyond the 
frequency at which the noise of the hot-wire becomes dominant. It can 
be shown theoretically that a hollow wire having the same outside diameter 
as the wire used here and heated to the same temperature would give a 
better signal-to-noise ratio. ‘Thus, by removing metal from within the 
wire, one could in theory extend the measurements of F(R), with no obvious 
limit until the atomic structure of the wire becomes significant. 
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This raises the question of other forms of noise and in particular of the 
noise produced by the gas itself. Somewhere, the spectrum of turbulence 
must merge with the spectrum of molecular agitation, and we shall now 
attempt to locate this point. 

Let us consider a volume of gas of linear dimension ZL containing 7 
molecules. ‘lhe molecular velocities are of the order of (k’ 7'/m)!*, where 
k’ is Boltzmann’s constant, 7’ the absolute temperature and m the mass of 
a molecule. ‘The velocity of the gas can be defined as that of the centre-of- 
mass of the 2 molecules and it will fluctuate with a mean-square velocity 
of the order of k’T/mn. Consider now two such elements of volume, with 
centres distance D apart. ‘The correlation between the velocities of the 
two centres of mass will vary from 1 if D = 0 to zeroif D > 2L. From the 
correlation we can compute the three-dimensional spectrum for RL < 1; 
in this way we find 


Evoise = (3k T/427p)k?, 
where p = nm/L* is the gas density. ‘his spectrum increases with k? and 
should reach a maximum near kL = 1. For the present purpose, we 
can assume that L is small (say L = 1 micron, or 50 mean-free-paths in 
air). 

If the one-dimensional spectrum of turbulence F falls as k~®, the three 
dimensional spectrum (Batchelor 1953) 1s given by £,,,), 48F. The 
spectrum of turbulence can therefore be expected to cross that of molecular 
agitation (so that F,,,,,, = Evie) at the wave-number where 

F(k) = (k'T/647*p)k?. (1) 


‘This expression should be regarded as a crude estimate. It is based ona 
simplified picture of molecular agitation which does not include fluctuations 
of pressure, density or temperature. It is also based on the application of 
the 'l'aylor hypothesis to the finest components of the turbulence. ‘This 
certainly ceases to be valid of k > Uv, since such eddies are dissipated before 
they are displaced by the mean flow. In the present case this occurs at 
k > 3000 cm™!. 

‘The expression on the right-hand side of (1) is shown as dotted lines 
on figures 4 and 6 and it suggests that, unless F(k) has an abrupt cut-off, the 
spectrum of turbulence disappears in molecular noise with a logarithmic 
slope of —6. ‘This would leave quantities such as ¢(0%u/éx*)?) without 
significance in the sense of fluid mechanics. Stating it a different way, it 
appears that the gap between turbulence and molecular agitation may not 
be as wide as is generally assumed. 
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SUMMARY 

In this work a method is given for finding the position of ring 
singularities in a three-dimensional potential field having axial 
symmetry, by consideration of the very much easier case of a 
similar two-dimensional potential function. It is seen that 
the traces of these ring singularities on a plane through the axis 
of symmetry occur at points corresponding to those of the 
singularities existing in the two-dimensional plane when the 
axial velocity potential functions are the same. It is thought 
that this might be of value in the plotting of stream surfaces. 


INTRODUCTION 

When an arbitrary velocity potential function for perfect fluid flow is 
assumed, it is essential that the singularities of this function or its 
derivatives should be excluded from the fluid region that it is hoped may 
be reproduced physically. In the following work a method has been 
developed for obtaining the positions of the singularities in velocities, 
or velocity potential, by consideration of similar velocity or velocity 
potential distributions in the two-dimensional case. It will be supposed 
that the velocity potential function, or one of its derivatives, is known 
along the axis of symmetry, and it can be shown that this defines the 
velocity potential function throughout all space. 

First it will be shown that a relationship exists between the two- 
dimensional case and the three-dimensional case with axial symmetry. 

The procedure for then obtaining the ring singularities when ¢ = f(x) 
along the axis in three dimensions, with symmetry assumed, will be to 
consider the position of the singularities existing when ¢ = f(x) along 
the x-axis in two dimensions. 


SIMILARITY OF CERTAIN TWO-DIMENSIONAL AND THREE-DIMENSIONAL 
VELOCITY POTENTIAL FUNCTIONS 
The axis of symmetry will be taken as the x-axis in both the two- and 
three-dimensional cases. ¢ denotes a velocity potential function in both 
cases, and the value of ¢ along the axis of symmetry is supposed to be 2f(x). 
If y is defined as the coordinate measured perpendicular to the x-axis 
in two dimensions and z = «+7y, it follows that f(z) = f(z), and so 


$ = f(z) +f(2) (1) 
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at any point of the two-dimensional field, with the usual notation for 
conjugate complex quantities. Here, f(z) is the analytic function of z 
that takes the real value f(v) when y = 0. 

For the three-dimensional case with axial symmetry, the definition of 
the axial velocity potential as 2f(x) determines the potential function 
throughout the whole of space except along a number of lines. It may be 
shown (Whittaker & Watson 1927, p. 388) that 

2 px 
= ~ | f(x +ircos @) dé (2) 


—/¥@ 

at every point of the field at which the integral is analytic, where r is defined 
as the coordinate measured perpendicular to the x-axis in three dimensions. 
Putting A = cos@ in (2), one obtains 
2 (1 f(x+irda) 27 sf) f(x +ird) "1 f(x—ird) : 
eee | aaa Os = =| | O57; eat sal | Tm |. (3) 
2 1 V(1—A*) TL Jo V(L—-A?*) Jo VQ —4*) 

There is clearly great similarity between equations (1) and (3). 

Suppose that f(z) has a singularity at s = «+78; then A = 1 gives the 
smallest value of r for which Ar = f, and when A= 1, f(x+zrd) has a 
singularity at x = a, r = f. 

When no finite limit f(«+78) exists, the integral (3) may not converge 
and then ¢ will be infinite at the singularity v = «, r = 8; note that the 
traces of this singularity are at the same points on a plane through the axis 
of symmetry as in the two-dimensional case. If however the integral (3) 
converges for the range A? <1 the velocity potential still exists. ‘The 
condition that the integral (3) converges at A = 1 is that 

san (LAM e+ iB) 

it V(IL—) 
is finite. Even in this case, however, it is easily seen that higher derivatives 
of ¢ must become infinite, so that the point is still a singularity. 

‘Taking the velocity in the direction of the x-axis to be u (= —0¢/éx), 
the preceding work may be repeated with f(x) replaced by —f’(x). Hence the 
singularities in the velocity distributions in the field follow a similar pattern. 





Some examples given below illustrate this. 


Example | 


Suppose the axial velocity is given as (1+.«?)"" = —2f"(x), say. 
Then, in two dimensions, 
od rdf 
1 F we Vt | Pm ‘ 
w= —- 2 = yseyn +t 2y". 
OX 
This becomes infinite for 7 0 when 27+1=0, ory = +1,x=0. 


In three dimensions 
f'(%+7BA) S[1+a«+7BA)?]-". 


In this case a: (1—Ayt 
; 5, 2[1+(a+7BA)?]" 
0 when « = 0, B = 1. 


1 


for i 
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Hence a ring singularity exists symmetrically about the x-axis, of 


radius r = 1 in the plane x = 0, and the value of w is not finite at points on 
this ring provided n—} 0. For the special case when n= 1, the 
integral (3) becomes 


| [’ dn [’ ar | 
m1 Jo ((x+irdA?+1]/(1—A?) © J 9 [(x—ard)? + 1] (1 — A’) |’ 
and this can be found explicitly as 


| 
st. ¥) = PF [{7? + (v—2)?h U2 — {72 + (4 +72)? }-17], 


which clearly has a singularity at x = 0, 7? = 1. It is, in addition, dis- 
continuous across the whole ‘branch line’ x = 0, r > 1. 


Example 2 

Suppose the axial velocity is given as (1+e”’)" = —2f'(x), say. 

Then, in two dimensions, uv = }{(1+e)"+(1+e)“"}; this becomes 
infinite for n 0 when e*°+1=0, or x =0, y =(2K+1)z, where K is 
an integer. 

For three dimensions, f’(a+7BA) = — $(1+e**'’*)-” and 

, 1(—A)!2 
we Tit area = L 
for n— 4 > 0 when « = 0 and B = (2K + 1)z. 

Hence ring singularities exist symmetrically about the x-axis, of radit 
r=, 3a, 52, ... (2K +1), in the plane x = 0. 

A variation of this axial velocity has been used by Szczeniowski (1943) 
for designing wind tunnel contractions; but he does not appear to know 
where these singularities occur. 


Example 3 

Suppose the axial velocity is given as 2(1 +.«')"' = —2f’(x), say. 

In two dimensions, u = (1 +24)!+(1+4+ 2"), which ceases to be analytic 
when 3'+ 1 = 0, that is, when zs = exp{}7(2K + 1)z}, where K is an integer 
and x+7y = (+1+72)/v2. 

For the three-dimensional case, 








f'(a+iBa) = —{14+(a+2Bayt 4, 
and fora = +1/v2, B= 4 
a) —-O-A | | et 
;.,|1+(%+z7BA) 5.4 4+(1400)4] 
Hence tor this axial velocity there exist two ring singularities, given by 
= +1/v2, °° = 4, at which the value of wv is not bounded. 


CONCLUSIONS 

The procedure of defining an axial velocity distribution in three 
dimensions has been used by several authors in the design of nozzles and 
wind tunnel contractions. Prior to this work, it has been necessary to 
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calculate the stream functions (or potential functions) and then deduce 
the regions of the fluid in which the velocity components tend to infinity. 
Since these functions have often been in series form this has not always 
been easy and has usually been tedious. It is suggested here that, following 
the adoption of the axial velocity potential function, the position of the 
singularities in the two-dimensional (and hence of the three-dimensional) 
field should be found, as this will give a good indication of the regions where 
sharp rises in velocity along stream surfaces may occur. If, then, it is of 
importance to eliminate such rises, the plotting of a few stream surfaces in 
the region of the singularities alone might prevent a considerable wastage 
of time. 

A point of interest is the reduction in degree of the singularities occurring 
in the three-dimensional case. For example, when the axial velocity is 
(1 + x”), then on the plane x = 0 the value of the u(0, y) in two dimensions 
is (1—y?)* whilst it can be shown that the value of u(0,r) in three 
dimensions is 

te! —r?) 3/24 (1 —r?) ait 2 
The ‘ order of infinity’ has therefore dropped from 2 to 3; and a reduction 
in order of } (compare example 1) seems general when going from three 
to two dimensions. 

Whilst the preceding work was done mainly for the consideration of 
perfect fluid flow, it is quite general and could be applied to other problems 
with axial symmetry. 
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REVIEWS 


Combustion Processes, edited by B. Lewis, R. N. Pease and H. S. 
‘TAYLOR (Volume II of High Speed Aerodynamics and Jet Propulsion). 
Princeton University Press, 1956. 662 pp. 512.50 or 84s. 


Almost every major advance in civilization has involved an increased 
understanding by man of the use of fire. ‘The use of the camp fire for 
preparation of food was probably one of man’s first steps away from animal 
ways. Fire was used as an aid in fashioning stone tools, and in the smelting 
of ores in the making of metal tools. ‘The discovery of gunpowder made 
untenable the types of civilization existent at the time of its introduction. 
Although the industrial revolution was made possible by the steam engine 
which used fire only in a rather simple furnace, modern society has 
developed as the industrial furnace developed. ‘The automobile and 
airplane are practical devices only because of the compact internal combus- 
tion engines. 

It might be supposed that the scientific, quantitative knowledge of 
combustion processes would have reached an advanced stage in view of 
the great importance of fire. Much indeed is known, and great effort is 
being expended to increase this knowledge, but in spite of all effort the 
available information is in large part empirical. 

The empirical nature of much of our knowledge and the amount of 
effort that has been put into the study of combustion processes is illustrated 
by the volume under review. ‘This book would have been impossible 50 years 
ago. We have 662 pages devoted to a brief account of the state of knowledge 
of combustion processes in gases, liquids, and solids for premixed fuel and 
oxidizer and in non-premixed systems. And yet in spite of the great 
progress made since the last century, this book makes amply clear, in 
nearly every chapter, the unsatisfactory state of our present knowledge and 
frequently points the way to the obviously needed further research. 

During the nineteenth century, while one group of physical chemists 
were struggling with the study of the kinetics of chemical reactions, another 
group were laying the foundations of thermodynamics and its application 
to chemical equilibrium. ‘The present volume very appropriately opens 
with a discussion of the essential relations from thermodynamics that 
permit the calculation of the equilibrium composition of reaction products, 
the heat release and the flame temperature. Only the briefest treatment is 
given since the entire first volume of this series is entitled Thermodynamics 
and Physics of Matter. 

In spite of the logical completeness of the understanding of chemical 
equilibrium afforded by thermodynamics, its practical usefulness is limited, 
as is made clear in the next two sections. From section C it is clear that, 
between the specification of a problem, involving a fuel and oxidizer 
composition, and the desired values of pressure, temperature, entropy 
and other properties of the combustion gases, lies a long path of 
computation. ‘This path, even for simple systems for which all thermal 
and equilibrium constants are available, is so involved that it is seldom 
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pursued unless large scale automatic computers are available. ‘The 
immediate future should see great progress in this area since the computing 
machines make practical many calculations which to date have been 
merely possible. 

It would seem that section B on “ Expansion Processes”? was placed 
between A on “ High Temperature Equilibrium ”’ and C on “* Computation 
Methods ’’ so that the reader would see immediately the limitations as 
well as the importance of the thermodynamic methods. Section B, in 
a sense, is the real meat of jet propulsion ; it is the analysis of the production 
of the jet from high temperature combustion gases. Besides the usual 
compressible flow discussion of the flow of a gas through a nozzle, there 
is treated the appropriate modifications required by the fact that chemical 
equilibrium may or may not be preserved during the expansion. ‘The 
additional problems presented by the fact that sometimes a condensed 
phase forms during the expansion are also considered. 

‘he enormous progress that has been made during the past century in 
the understanding of the kinetic rate mechanism for various chemical 
processes is attested by the 60 pages devoted to the fundamentals of this 
subject. A short text on chemical kinetics would be a good name for 
these chapters. ‘They start with an empirical introduction to reaction rate 
theory, introducing order, half life, temperature, rate constants, activation 
energy. ‘This is followed by the theoretical understanding of the empirical 
facts through kinetic collision theory and absolute rate theory. By 
considering a series of specific cases, the ideas of free radicals, chain 
mechanisms, explosions, wall reactions, solid and gas reactions and 
experimental techniques are presented. Finally the bearing of all these 
effects on non-equilibrium in systems involving rapid changes is explained. 

With the above as introduction, several reactions of importance in 
propulsion are studied in detail. ‘The reader is brought up to date in our 
understanding of the hydrogen-oxygen, carbon monoxide-oxygen and 
hydrocarbon-oxygen systems. ‘lhe uninitiated reader learns of the great 
complexity of even the simplest of these reactions and is thus brought to 
appreciate the present-day potentialities and limitations in the analysis of 
combustion systems. A brief treatment of several systems of interest in 
modern propulsion is included: boron hydrides, hydrogen—fluorine and 
nitromethane decomposition. 

One more section on fundamental ideas is included, namely, ‘‘ The 
Mechanics of Reacting Continua’’. In this 11 page section the reader is 
taken through the ideas and derivations of the continuity, momentum and 
energy equations for flowing, reacting, viscous, heat conducting 
compressible fluids in general vector form. Only 7 pages are required for 
this, the last 4 pages of the chapter being devoted to a discussion of 
irreversible processes, including entropy transport, forces, currents, and 
the kinetic coefficients, and the introduction of the Onsager reciprocal 
relations. ‘This section is an excellent very concise review of fundamentals 
from a general viewpoint for readers already well prepared in advanced 
theoretical fluid mechanics. As a text it’s next to useless. 


—EEE 
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In present day combustion, a careful observer notes a discontinuity 
between the discussions of fundamentals and of applications. ‘The 
fundamentals are rather precisely known physical principles, worked out 
in large part in precise mathematical form by carefully reasoned mathe- 
matical arguments. ‘The applications on the other hand consist of many 
empirical facts ‘ explained’ by a series of ad hoc assumptions assisted by 
the fundamentals as by a crutch. ‘his discontinuity was very sharp half 
a century ago. ‘l’oday—as if by a diffusion process—the discontinuity 
has become dulled. ‘The detailed nature of the blurred discontinuity in 
ideas and methods occupies the last two thirds of this volume. 

Almost all of the middle third of the book is devoted to the combustion 
of gases. First the phenomena encountered in non-turbulent gases are 
discussed. A rather brief account of ‘ the flame’ is given which indicates 
the qualitative ideas of flame thickness, ignition temperature, etc. Only 
the simplest and least adequate relation for the flame speed, that of Mallard 
and LeChatelier, is derived, although a very adequate bibliography of the 
literature on more adequate flame speed relations is included. ‘There 
follows a discussion of ideas and assumptions and corresponding derivations 
on a host of important topics; quenching distance, flame stabilization, 
blowoff and flashback, flame shape and fluid velocity field, propagation 
in channels, ignition, limits of inflammability, minimum ignition energy, 
etc. For those readers who approach combustion from the more exact 
field of aerodynamics, this chapter will leave them reeling. However, 
in most cases combustion is complicated, our understanding is not complete, 
and those who desire a more exact understanding should regard these 
phenomena as challenging and the ad hoc assumptions as the best treatment 
yet possible, a zeroth order approximation, so to speak. Some 20 pages 
devoted to quantitative experimental results on laminar flames serve to 
make these available for use and to show the need for more analytical studies. 

The problems associated with the turbulent flame are, as would be 
expected, more complex and less well established than for laminar flames. 
In the turbulent case, even the simplest ideas, such as the flame ‘ speed ’, 
are hard to define and measure. ‘The fact that a conflict rages between the 
proponents of the wrinkled laminar flame sheet and the proponents of the 
extended reaction zone is not surprising. ‘The author of this section (H) 
is a wrinkled flame sheet man, and most of the 50 pages of the section is 
devoted to an examination of available data in the light of the most modern 
form of this theory. It should be stated that so far the flame sheet wrinklers 
have succeeded in deriving quantitative consequences of their ideas in 
a reasonably basic way, while the extended reaction zone proponents have 
not. In the reviewer’s opinion, the turbulent flame contains significant 
proportions of both wrinkled flame sheets and extended reaction zones. 
The success of one group in deriving formulas should not be interpreted 
as showing that turbulent flames do not involve extended reaction zones. 

The diffusion flame, i.e. the combustion of a jet of fuel in air (say), 
is controlled more by the aerodynamic mixing processes than by the 
chemical rates. Hence there is a clear separation into laminar, transition 
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and turbulent diffusion flames. ‘The experimental results are compared 
with present theoretical ideas of the parts played by molecular and turbulent 
diffusion. Again the understanding, as for example of flame height, is 
adequate to zero order but leaves much to be learned in the future. 

Liquid fuels are usually burned by first atomizing and evaporating 
them. ‘The study of such processes has yielded more precise knowledge 
than some other areas of combustion. However the added complexity 
still leaves much for future study. ‘lhe same is true of solid fuels both 
in lumps and as dust. ‘The treatment given both of these subjects is such 
as to bring the uninitiated reader a long way toward the understanding he 
will need for further work in these areas. 

The application of the fundamental and empirical knowledge of 
combustion processes to jet propulsion is discussed in two long sections, 
one on liquid propellants and one on solid propellants. As is so often the 
case, the knowledge of fundamental processes serves more as a guide than 
as a quantitative calculation procedure on complete engineering systems. 
The complexity is of too high an order. Monopropellants, bipropellants 
and common solid propellants are considered. Equipment, ignition and 
ignition delay problems, chemical and physical mechanisms, pressure 
dependence, thermal decomposition, and performance characteristics 
are all treated. ‘The performance attainable and the problems remaining 
are clearly presented. 

The chemical part of this book on combustion closes with a qualitative 
discussion of detonation processes. Because of the small thickness of the 
detonation front, the overall gas dynamics of the detonation process is 
perhaps simpler than the more extended combustion zones of ordinary 
fires but the chemical and physical processes are at least as complex. Often 
the chemical and physical nature of explosives is complex and the processes 
confined to narrow zones are hard to study experimentally. ‘The parts 
played by composition, reaction time, lateral expansion, detonation 
characteristics, including initiation and decay, are presented from a physical 
chemist’s (rather than aerodynamicist’s) point of view. 

The last section gets into a book on combustion as a most modern 
touch, “ Energy Production by Nuclear Reactions”. ‘This section in 
a book on combustion processes in a series on high speed flight and jet 
propulsion suggests the inference that the future development of rocket 
vehicles both for terrestrial and space flight could be influenced in a major 
way by nuclear developments. ‘There are some close relationships between 
combustion and nuclear energy production from both the hydrodynamic 
and transport mechanism points of view. ‘The present chapter of 40 pages 
can, however, do no more than review the elementary physical principles, 
and this it does in a simple and clear manner. 

This book is the best presentation of the whole of the knowledge of 
combustion currently available. ‘The general excellence of the selection 
of material and its presentation will leave this book in first place until the 
growth of the field makes the present understanding inadequate. 


H, W. EMMons 








~ 


ab 








emuan ee ee AIS ian eialarls is , 








JOURNAL OF FLUID MECHANICS 


Volume 3, Part-2 November 1957 


CONTENTS 


The fundamental solution for small steady three-dimensional 
disturbances to a two-dimensional parallel shear flow 


by M. J. LIGHTHILL 


Skin friction on a strip of finite width moving parallel to its 
length 
by HAROLD LEVINE - % 7 “s 


The asymptotic expansions at large Reynolds numbers for 
steady motion between non-coaxial rotating cylinders 


by W. W. WOOD 


The supersonic motion of an aerofoil through a temperature 
front 


by J. W. CRAGGS 


On the generation of surface waves by shear flows 
by JOHN W. MILES 


On the fine structure of turbulent flows 
by R. BETCHOV 


Note on the position of ring singularities in an axisymmetric 
potential field 
by N. T. BLOOMER 


Reviews . 


Printed in Great Britain by Taylor & Francis Ltd. 





113 


145 


Ny 
bad 
“I 


221 











